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Ephemeral-Stream  Channel  and  Basin-Floor  Infiltration 
and  Recharge  in  the  Sierra  Vista  Subwatershed  of  the 
Upper  San  Pedro  Basin,  Southeastern  Arizona 


By  A.L.  Coes  and  D.R.  Pool 

Abstract 

The  timing  and  location  of  streamflow  in  the  San  Pedro 
River  are  partially  dependent  on  the  aerial  distribution 
of  recharge  in  the  Sierra  Vista  subwatershed.  Previous 
investigators  have  assumed  that  recharge  in  the  sub  watershed 
occurs  only  along  the  mountain  fronts  by  way  of  stream- 
channel  infiltration  near  the  contact  between  the  low- 
permeability  rocks  of  the  mountains  and  the  basin  fill.  Recent 
studies  in  other  alluvial  basins  of  the  Southwestern  United 
States,  however,  have  shown  that  significant  recharge  can 
occur  through  the  sediments  of  ephemeral- stream  channels 
at  locations  several  kilometers  distant  from  the  mountains. 

The  purpose  of  this  study  was  to  characterize  the  spatial 
distribution  of  infiltration  and  subsequent  recharge  through  the 
ephemeral-stream  channels  in  the  Sierra  Vista  subwatershed. 

Infiltration  fluxes  in  ephemeral- stream  channels  and 
through  the  basin  floor  of  the  subwatershed  were  estimated 
using  several  methods.  Data  collected  during  the  drilling 
and  coring  of  16  boreholes  included  physical,  thermal,  and 
hydraulic  properties  of  sediments;  chloride  concentrations  of 
sediments;  and  pore- water  stable-isotope  values  and  tritium 
activity.  Surface  and  subsurface  sediment  temperatures  were 
continuously  measured  at  each  borehole. 

Twelve  boreholes  were  drilled  in  five  ephemeral- 
stream  channels  to  estimate  infiltration  within  ephemeral- 
stream  channels.  Active  infiltration  was  verified  to  at  least 
20  meters  at  1 1  of  the  12  borehole  sites  on  the  basis  of  low 
sediment-chloride  concentrations,  high  soil- water  contents, 
and  pore- water  tritium  activity  similar  to  that  of  present- 
day  precipitation.  Consolidated  sediments  at  the  twelfth 
site  prevented  core  recovery  and  estimation  of  infiltration. 
Analytical  and  numerical  methods  were  applied  to  determine 
the  surface  infiltration  flux  required  to  produce  the  observed 
sediment- temperature  fluctuations  at  six  sites.  Infiltration 
fluxes  were  determined  for  summer  ephemeral  flow  events 
only  because  no  winter  flows  were  recorded  at  the  sites  during 
the  monitoring  period. 

Four  boreholes  were  drilled  in  the  basin  floor  to 
estimate  infiltration  in  areas  between  ephemeral- stream 
channels.  Infiltration  fluxes  through  the  basin  floor  ranged 


from  less  than  1  centimeter  to  6  centimeters  per  year. 

At  a  site  in  semiconsolidated  to  consolidated  basin-fill 
conglomerate,  the  long-term  infiltration  fluxes  were  very 
low  (less  than  1  centimeter  per  year).  Chloride,  tritium, 
and  stable-isotope  data  indicate  long  periods  of  no  net  deep 
downward  percolation  flux  beneath  the  basin  floor.  At  a  site 
in  unconsolidated  to  semiconsolidated  basin-fill  sand  and 
gravel,  infiltration  fluxes  were  high  (2  to  6  centimeters  per 
year).  Chloride,  tritium,  and  stable-isotope  data  indicate  active 
infiltration  to  8  meters  and  a  decrease  in  infiltration  below 
8  meters.  The  change  in  the  infiltration  rate  below  8  meters 
is  controlled  by  an  increase  in  the  silt  and  clay  content  of  the 
sediment. 

Ephemeral- stream  channel  recharge  for  the  entire 
subwatershed  was  estimated  by  upscaling  the  calculated 
infiltration  fluxes  and  weighting  the  fluxes  by  streamflow 
duration,  evaporation,  and  transpiration.  In  contrast  to 
previous  assumptions,  recharge  from  ephemeral- streamflow 
infiltration  occurs  not  only  near  the  mountain  fronts,  but 
also  along  significant  lengths  of  ephemeral- stream  channels. 
Although  most  of  the  ephemeral  streams  in  the  subwatershed 
flow  less  than  a  few  days  per  year,  the  available  streamflow 
quickly  infiltrates  past  depths  where  it  is  available  for 
evapotranspiration.  This  water  likely  stays  in  the  unsaturated 
zone  until  it  is  vertically  displaced  by  infiltrated  water  from 
subsequent  streamflows  and  eventually  recharges  the  regional 
aquifer.  Ephemeral-stream  channel  infiltration  during  2001 
and  2002  was  estimated  to  account  for  about  12  to  19  percent 
of  the  estimated  average  annual  recharge  in  the  Sierra  Vista 
subwatershed. 


Introduction 

An  improved  understanding  of  the  recharge  distribution 
in  the  Sierra  Vista  sub  watershed  of  the  Upper  San  Pedro 
ground-water  basin  was  needed  to  better  estimate  the  effects 
of  ground- water  withdrawals  on  natural  discharges  from 
the  system  through  base  flow  and  riparian  vegetation.  The 
interception  of  ground  water  that  would  flow  naturally  from 
recharge  areas  to  the  San  Pedro  River  by  pumping  could  result 


2  Ephemeral-Stream  Channel  and  Basin-Floor  Infiltration  and  Recharge  in  the  Sierra  Vista  Subwatershed,  Arizona 


in  decreased  streamflows  in  the  river.  The  timing  and  location 
of  the  reduced  streamflow  are  dependant  on  the  location  of 
the  withdrawal  wells,  the  hydraulic  properties  of  the  aquifer, 
and  the  aerial  distribution  of  recharge.  Recharge  occurring 
downgradient  of  withdrawal  wells  is  less  likely  to  be  captured 
by  withdrawal  wells  and  more  likely  to  discharge  to  the 
river.  Improved  understanding  of  recharge  distributions  will 
enable  local,  State,  and  Federal  organizations  and  agencies  to 
effectively  manage  available  water  resources  to  meet  water 
needs  in  the  basin.  Toward  this  goal,  the  U.S.  Geological 
Survey  (USGS),  in  cooperation  with  the  Bureau  of  Land 
Management  (BLM)  and  the  Upper  San  Pedro  Partnership, 
has  collected  and  analyzed  hydrologic  information  describing 
infiltration  and  recharge,  ground-  and  surface-water 
interactions,  water-level  and  ground- water  storage  trends,  and 
geometry  and  lithology  of  the  aquifer.  This  report  documents 
the  results  of  investigations  of  recharge  by  stream-channel 
and  basin-floor  infiltration,  potentially  important  recharge 
mechanisms  in  the  basin. 

The  Sierra  Vista  subwatershed  encompasses  about 
2,460  square  kilometers  (km2)  in  the  southern  part  of  the 
San  Pedro  River  Basin  in  southeastern  Arizona  (fig.  1)  and 
extends  about  43  km  north  from  the  international  boundary 
with  Mexico  to  near  Fairbank.  The  subwatershed  is  part  of 
a  broad  alluvial  basin  and  is  bounded  on  the  west  by  the 
Huachuca  Mountains  and  on  the  east  by  the  Mule  Mountains 
and  the  Tombstone  Hills.  Elevations  of  the  basin  floor  vary 
from  about  1,150  to  1,310  m;  mountain  elevations  range  from 
about  1,520  to  2,900  m. 

The  San  Pedro  River  flows  northward  through  the 
subwatershed  (fig.  1).  The  river  and  its  associated  ground- 
water  system  support  a  corridor  of  riparian  vegetation  that 
hosts  several  endangered  species  and  is  an  important  habitat 
for  American  migratory  birds.  About  64  km  of  the  river 
and  riparian  area  has  been  protected  by  designation  as  the 
San  Pedro  Riparian  National  Conservation  Area  (SPRNCA; 
fig.  1),  and  is  managed  by  the  BLM. 

Sierra  Vista,  with  a  population  of  37,775  in  2000,  is 
the  largest  city  in  the  sub  watershed  (U.S.  Census  Bureau, 
2000).  By  2020,  the  population  of  Sierra  Vista  is  projected 
to  be  52,571  (Arizona  Department  of  Economic  Security, 
1997),  an  increase  of  39  percent.  Populations  of  other  cities 
in  the  sub  watershed  in  2000  were  6,090  in  Bisbee,  1,751  in 
Huachuca  City,  and  1,504  in  Tombstone  (U.S.  Census  Bureau, 
2000). 

Ground  water  is  the  primary  source  for  municipal, 
domestic,  livestock,  and  irrigation  water  needs  in  the 
subwatershed,  and  ground- water  withdrawals  have  increased 
as  the  population  has  increased.  The  magnitude  of  the  effect 
that  these  increasing  withdrawals  have  on  the  river  and  its 
associated  riparian  vegetation  is  proportional  to  the  amount 
and  distribution  of  recharge  to  the  ground- water  system. 
Ground- water  pumping  intercepts  part  of  the  natural  flow  of 
water  through  the  aquifer  before  the  water  reaches  natural 
discharge  areas.  To  better  understand  and  predict  the  effects 
of  historic,  present,  and  future  ground- water  withdrawals 
on  natural  discharge,  the  spatial  and  temporal  variability  of 
recharge  rates  must  be  adequately  understood. 


Previous  investigators  have  assumed  that  recharge  in 
the  Sierra  Vista  subwatershed  occurs  only  near  the  mountain 
fronts  (Freethey,  1982;  Corell  and  others,  1996;  Goode 
and  Maddock,  2000).  They  envisioned  recharge  occurring 
as  stream-channel  infiltration  near  the  contact  of  the  low- 
permeability  rocks  of  the  mountains  and  the  basin  fill.  Recent 
studies  in  other  alluvial  basins  of  the  Southwestern  United 
States,  however,  have  shown  that  significant  recharge  can 
occur  through  the  sediments  of  stream  channels  far  from  the 
mountains  (Blasch  and  others,  2000;  Izbicki  and  others,  2000; 
Constantz  and  others,  2003).  Direct  infiltration  of  precipitation 
through  the  floors  of  alluvial  basins  in  the  Southwest  has 
been  considered  negligible  because  of  high  evapotranspiration 
and  low  precipitation.  The  possibility  of  significant  recharge 
along  ephemeral-stream  channels  far  from  the  mountains  or 
through  direct  infiltration  on  the  basin  floor  places  uncertainty 
on  the  adequacy  of  previous  recharge  concepts  for  use  in  the 
development  of  water-management  plans. 

Purpose  and  Scope 

This  report  presents  estimates  of  infiltration  and  recharge 
through  the  ephemeral- stream  channels  and  the  basin  floor 
of  the  Sierra  Vista  subwatershed,  and  documents  the  data 
collection  and  analysis.  Data  resulted  from  drilling  and  coring 
of  16  boreholes  and  included  physical,  hydraulic,  and  thermal 
properties  of  sediment;  chloride  concentrations  of  sediment; 
pore- water  stable-isotope  values  and  tritium  activity;  neutron- 
soil  moisture  and  electromagnetic  induction  (EMI)  logs; 
and  continuous  sediment-temperature  profiles  and  water 
levels  during  2001  and  2002.  Unsaturated-zone  infiltration 
processes  at  six  ephemeral- stream  channel  sites  were  modeled 
to  estimate  infiltration  and  percolation  fluxes  during  and 
following  streamflows. 
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Figure  1.  Sierra  Vista  subwatershed,  locations  of  boreholes,  temperature  sensors.  Agricultural  Research  Service  flumes, 

U.S.  Geological  Survey  streamflow-gaging  stations,  and  precipitation  gages,  southeastern  Arizona  and  northeastern  Sonora,  Mexico. 
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Land  Cover 

Major  types  of  land  cover  include  grassland 
(32.6  percent),  desert  scrub  (28.8  percent),  mesquite 
woodland  (13.5  percent),  and  urban  (5.3  percent;  fig.  2; 

U.S.  Environmental  Protection  Agency,  2003).  The  remaining 
19.8  percent  includes  oak  woodlands  (15.2  percent)  and 
forest  (2.2  percent)  in  and  near  the  mountains  and  the  riparian 
(0.9  percent)  and  agricultural  (0.8  percent)  areas  near  the 
San  Pedro  and  Babocomari  Rivers  (U.S.  Environmental 
Protection  Agency,  2003).  Although  the  urban  percentage  is 
small,  most  of  the  urban  land  is  concentrated  in  the  north- 
central  part  of  the  subwatershed. 

The  hydrologic  influence  of  each  major  land-cover 
type  on  the  portion  of  precipitation  that  is  transported  to 
ephemeral  stream  channels  has  not  been  qualified.  Some  basic 
assumptions  can,  however,  be  inferred.  Runoff  produced  in 
mesquite  woodlands  and  desert  scrub  areas  most  likely  is  a 
small  percentage  of  precipitation  because  leaves  and  branches 
can  intercept  precipitation  and  the  high  soil  temperatures 
and  low  soil-moisture  contents  and  potentials  result  in 
high  soil  permeability.  Runoff  from  grasslands  likely  is  a 
greater  percentage  of  precipitation  because  the  lower  soil 
temperatures  and  higher  soil-moisture  contents  result  in  lower 
soil  permeability  than  that  in  mesquite  woodlands  and  desert 
scrub  areas.  Urban  areas  produce  the  highest  percentage 
of  runoff  into  stream  channels  because  of  large  expanses 
of  impermeable  surfaces  that  prevent  soil  infiltration  and 
engineered  drainage  infrastructures  that  concentrate  and 
focus  runoff. 

Soils 

The  distribution  of  soil  types  in  the  Sierra  Vista 
sub  watershed  is  dominated  by  lithic,  loamy,  and  fine 
paleosoils  near  the  mountains,  and  clayey  soils  and  calcium 
carbonate  soils  along  the  lower  reaches  of  the  watersheds 
(fig.  3;  U.S.  Environmental  Protection  Agency,  2003).  Narrow 
ribbons  of  sandy  loam  typically  lie  along  ephemeral- stream 
channels;  sandy,  clayey,  and  calcium  carbonate  soils  typically 
lie  between  the  channels.  Calcium  carbonate  soils  dominate 
areas  below  the  mountains  across  much  of  the  eastern  part  of 
the  subwatershed.  A  wedge  of  calcium  carbonate  soils  also 
extends  nearly  to  the  mountains  in  the  central  part  of  the  study 
area  west  of  the  San  Pedro  River.  The  southwestern  part  of  the 
subwatershed  is  dominated  by  clayey  and  sandy  soils. 

Soil  types  can  also  have  a  strong  influence  on  the  fraction 
of  precipitation  that  is  transported  to  ephemeral-stream 
channels.  Sandy  loam  soils  likely  are  the  most  permeable, 
enabling  direct  infiltration  of  precipitation  and  limiting  runoff 
into  stream  channels.  Clayey  soils  can  limit  deep  infiltration 
where  they  underlie  channels.  Soils  having  secondary  calcium 
carbonate  are  less  permeable  and  enhance  runoff.  Shallow 
lithic  soils  underlain  by  bedrock  or  rocky  layers  promote 
runoff  into  stream  channels  and  limit  infiltration. 


Precipitation 

The  subwatershed  climate  is  arid  to  semiarid. 

Precipitation  data  are  available  from  the  Western  Regional 
Climate  Center  (2003)  for  a  gage  in  Tombstone  at  an 
elevation  of  1,384  m,  from  1905  to  2002;  and  for  a  gage  in  the 
Huachuca  Mountains  at  an  elevation  of  1,692  m,  from  1956  to 
2002  (fig.s  1  and  4).  Historically,  July  through  September  are 
the  wettest  months;  precipitation  during  this  period  averages 
21  cm  at  the  Tombstone  gage  and  27  cm  at  the  Huachuca 
Mountains  gage.  Precipitation  from  October  through  March 
averaged  12  and  23  cm  at  the  Tombstone  and  Huachuca 
Mountains  gages,  respectively;  and  precipitation  from  April 
through  June  averaged  less  than  3  cm  at  both  gages. 

Long-term  monitoring  indicates  that  annual  precipitation 
rates  were  above  average  before  about  1940  and  during  the 
early  to  mid-1980s,  below  average  from  about  1940  to  1980, 
and  about  average  after  the  mid-1980s  (Pool  and  Coes,  1999). 
Trends  in  seasonal  precipitation,  however,  are  different  from 
trends  in  annual  precipitation.  From  1956  to  1997,  winter 
(November-February)  precipitation  rates  increased  while  wet- 
season  (June-October)  precipitation  rates  decreased  (Pool  and 
Coes,  1999). 

Estimates  of  recharge  in  this  report  are  biased  by 
precipitation  immediately  prior  to  and  during  the  monitoring 
period.  Precipitation  during  the  year  that  preceded  the 
monitoring  period,  2000,  was  average  (36  cm)  at  the 
Tombstone  gage,  but  wetter  than  average  (78  cm)  at  the 
Huachuca  Mountains  gage  (fig.  4).  During  2000,  the  Huachuca 
Mountains  gage  recorded  an  abnormally  wet  June  and  August, 
and  October  was  the  wettest  month  on  record  (32  cm). 
Precipitation  during  the  monitoring  period,  2001  and  2002, 
was  dominated  by  dry  winters  and  slightly  drier  than  average 
summers  (fig.  4).  During  July  through  September  of  2001 
and  2002,  precipitation  was  below  average  at  the  Tombstone 
gage  (16  and  17  cm,  respectively)  and  slightly  above  average 
at  the  Huachuca  Mountains  gage  (29  cm  during  both  years). 
Precipitation  in  October  through  June  of  2001  and  2002  was 
below  average  at  both  gages. 

Surface-Water  Drainage  System 

The  San  Pedro  River,  an  intermittent  stream,  drains  the 
Sierra  Vista  sub  watershed  (fig.  1).  Annual  and  wet-season 
runoff  in  the  San  Pedro  River  have  declined  since  the  mid- 
1910s,  probably  as  a  result  of  reduced  precipitation  duration 
and  intensity,  increased  vegetation,  and  increased  streamflow 
infiltration  along  ephemeral  reaches  of  the  river  and  its 
tributary  streams  (Pool  and  Coes,  1999).  Base  flow  in  the 
river  is  lowest  during  the  summer  when  rates  of  near- stream 
ground-water  withdrawals  by  wells  and  phreatophytes  are 
highest.  Long-term  records  indicate  that  summer  base  flows 
have  been  declining  since  1937  (Pool  and  Coes,  1999). 
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Figure  2.  Land  cover  and  locations  of  borehole  sites.  Sierra  Vista  subwatershed,  Arizona. 
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Figure  3.  Soil  type  and  locations  of  borehole  sites,  Sierra  Vista  subwatershed,  Arizona. 
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Figure  4.  Monthly  and  annual  precipitation.  Tombstone  and  Huachuca  Mountains  gages.  Sierra  Vista  subwatershed,  Arizona. 


Tributaries  of  the  San  Pedro  River  are  ephemeral  except 
for  the  Babocomari  River,  which  is  perennial  in  some  areas 
(fig.  5).  Flow  occurs  only  in  direct  response  to  rainfall  or 
snowmelt.  The  amount  of  flow  is  generally  dependant  on  the 
spatial  distribution  and  intensity  of  precipitation,  the  runoff 
characteristics  of  the  channel,  and  the  associated  drainage 
area.  The  temporal  distribution  of  flow  in  ephemeral- stream 
channels  in  the  arid  Southwest  generally  is  the  result  of 
monthly,  annual,  and  decadal  oscillations  in  precipitation 
frequency  (Redmond  and  Koch,  1991;  Webb  and  Betancourt, 
1992). 

The  majority  of  tributaries  to  the  San  Pedro  River 
originate  in  the  mountains  bounding  the  subwatershed  and 
coalesce  downstream  to  form  higher-ordered  ephemeral- 
stream  channels  (fig.  5).  Tributaries  generally  have  narrow, 
shallow  channels,  close  to  the  mountain  fronts,  that  become 
wide,  deep  channels  below  a  knickpoint.  Channels  draining 
the  east  side  of  the  Huachuca  Mountains  typically  have  long 
and  narrow  drainage  areas  and  few  channels  coalescing 
downstream;  the  southern  tributaries  become  wider  and  less 
channelized  towards  the  San  Pedro  River.  Channels  draining 


the  west  side  of  the  Mule  Mountains  typically  have  much 
wider  drainage  areas  and  multiple  channels  coalescing 
downstream. 

Ground-Water  System 

The  ground-water  system  of  the  Sierra  Vista 
subwatershed  has  been  described  by  previous  investigators 
(Brown  and  others,  1966;  Pool  and  Coes,  1999);  an  overview 
is  given  here.  The  primary  aquifer  in  the  subwatershed  is 
within  basin  fill  that  overlies  crystalline  and  sedimentary 
rocks.  Secondary  aquifers  comprise  stream  alluvium 
beneath  the  flood  plains  of  the  San  Pedro  and  Babocomari 
Rivers,  and  pre-basin  fill  and  limestones  that  crop  out  in  the 
mountains  and  hills  surrounding  the  subwatershed.  Ground 
water  generally  flows  from  recharge  areas  to  discharge 
areas  along  the  San  Pedro  and  Babocomari  Rivers.  Ground 
water  also  discharges  from  the  aquifer  through  ground- water 
withdrawals,  evapotranspiration,  springs,  and  ground- water 
underflow  to  the  north. 
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drainage  basins  refers  to  site  listed  in  table  1 

Figure  5.  Tributary  drainage  basins  and  locations  of  borehole  sites,  Sierra  Vista  subwatershed,  Arizona  and  northeastern 
Sonora,  Mexico. 
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Aquifers 

The  Pantano  Formation,  a  consolidated  to 
semiconsolidated  conglomerate,  overlies  the  bedrock  in  the 
sub  watershed  and  is  the  oldest  of  the  basin  sediments.  The 
formation  probably  is  greater  than  900  m  thick  in  some  areas 
(Halverson,  1984;  Gettings  and  Houser,  1995).  The  Pantano 
Formation  is  a  locally  important  water-bearing  unit,  as  it  yields 
water  through  fractures  to  many  wells  in  the  Sierra  Vista  area 
(Pool  and  Coes,  1999). 

The  basin  fill  that  overlies  the  Pantano  Formation  has 
been  divided  into  upper  and  lower  units  on  the  basis  of 
compositional,  textural,  and  depositional  differences.  The 
lower  unit  consists  of  partially  cemented  interbedded  gravel, 
sand,  silt,  and  clay  that  is  often  described  as  conglomerate, 
and  ranges  in  thickness  from  about  45  to  90  m  (Pool  and 
Coes,  1999).  Permeability  of  the  lower  unit  varies  because  of 
variation  in  sediment  size,  sorting,  and  cementation  (Brown 
and  others,  1966).  The  upper  unit  consists  of  poorly  cemented 
gravel,  sand,  silt,  and  clay  (Brown  and  others,  1966),  and  is 
generally  less  than  120  m  thick  (Pool  and  Coes,  1999).  The 
upper  unit  includes  a  permeable  fan-gravel  facies  near  the 
mountain  fronts  that  grades  to  a  poorly  permeable  silt  and  clay 
facies  with  interspersed  sand  and  caliche  beds  near  the  basin 
center  (Pool  and  Coes,  1999).  Between  Sierra  Vista  and  the 
San  Pedro  River,  and  between  Hereford  and  Highway  90,  the 
upper  unit  is  primarily  a  saturated  confining  bed  of  silt  and 
clay  (Pool  and  Coes,  1999).  The  upper  and  lower  basin-fill 
units  are  the  major  water-yielding  units  in  the  Sierra  Vista 
subwatershed  and  constitute  the  basin-fill  aquifer.  Water  levels 
in  the  basin-fill  aquifer  range  from  less  than  1  m  to  as  much  as 
165  m  (USGS  unpublished  data). 

Terrace  deposits  of  clay,  silt,  sand,  and  gravel  overlie  the 
upper  basin  fill  in  most  areas.  These  deposits  form  a  veneer 
of  stream  alluvium  that  is  thin  near  the  mountain  fronts  but 
is  as  much  as  15  to  30  m  thick  in  erosional  channels  that 
parallel  the  San  Pedro  River  (Pool  and  Coes,  1999).  Older 
terrace  deposits  include  lakebed  deposits  of  silt,  clay,  and 
marl  (Haynes,  1968)  that  are  local  confining  layers  where 
they  occur  near  the  San  Pedro  River;  erosion  has  removed 
most  of  the  lakebed  deposits  along  tributaries  (Pool  and  Coes, 
1999).  Younger  terrace  deposits  include  stream  alluvium  along 
the  San  Pedro  River  and  some  of  its  tributaries  (Pool  and 
Coes,  1999).  The  older  deposits  are  rarely  saturated  outside 
of  the  flood  plains  of  the  San  Pedro  and  Babocomari  Rivers. 
The  younger  deposits  are  saturated  in  some  areas  near  the 
San  Pedro  and  Babocomari  Rivers  and  along  some  tributaries 
near  the  base  of  the  mountains  and  are  locally  important 
water-bearing  units  (Pool  and  Coes,  1999). 

Recent  stream  alluvium  lies  along  the  San  Pedro  and 
Babocomari  Rivers.  Stream  alluvium  deposited  prior  to  the 
entrenchment  of  the  San  Pedro  River  consists  of  clay,  silt,  and 
fine  sand  interbedded  with  coarse  sand  and  gravel,  and  is  as 
much  as  6  m  deep  and  1.5  km  wide  (Pool  and  Coes,  1999). 
Stream  alluvium  deposited  after  the  entrenchment  of  the 
San  Pedro  River  consists  of  sand  and  gravel  and  is  less  than 


6  m  deep  and  only  about  1  km  wide  (Pool  and  Coes,  1999). 
The  post-entrenchment  stream  alluvium  is  highly  permeable; 
the  pre-entrenchment  stream  alluvium  is  less  permeable  (Pool 
and  Coes,  1999).  The  pre-  and  post-entrenchment  stream 
alluvium  forms  a  long,  narrow,  shallow  aquifer  along  the 
San  Pedro  River  that  is  referred  to  as  the  flood-plain  aquifer. 
Water  levels  in  the  flood-plain  aquifer  range  from  less  than 
1  m  to  as  much  as  6  m  below  land  surface  (USGS  unpublished 
data).  Water  levels  in  this  aquifer  are  a  major  influence  on 
river  discharges  and  show  seasonal  trends  reflecting  cycles  of 
surface-water  runoff  and  evapotranspiration. 

Recharge  Mechanisms 

Historically,  recharge  to  the  basin-fill  aquifer  in  the 
Sierra  Vista  subwatershed  was  thought  to  primarily  occur 
near  the  mountain  fronts  through  ephemeral- stream  channel 
infiltration  of  runoff  from  the  relatively  impermeable  rocks 
of  the  mountains.  Previous  conceptual  models  of  the  ground- 
water  flow  system  included  annual  recharge  rates  of  15  to 
23  hm3/yr,  all  of  which  was  distributed  within  about  1.5  km 
of  the  mountains  (Freethey,  1982;  Corell  and  others,  1996; 
Goode  and  Maddock,  2000).  Infiltration  and  recharge  could, 
however,  be  distributed  along  much  greater  lengths  of  the 
channels.  Recent  studies  (Blasch  and  others,  2000;  Izbicki 
and  others,  2000;  Constantz  and  others,  2003)  have  shown 
that  significant  recharge  can  occur  through  the  beds  of 
ephemeral  streams  far  from  mountain  fronts.  In  the  Sierra 
Vista  sub  watershed,  Corell  and  others  (1996)  estimated  that 
infiltration  along  the  length  of  Greenbush  Draw  contributed 
about  1.20  hm3/yr  of  recharge.  Goodrich  and  others  (2004) 
estimated  that  a  6.8-km  reach  of  Walnut  Gulch  contributed 
about  0.14  to  0.47  hm3  of  recharge  in  1999  and  about  0.13  to 
0.37  hm3  of  recharge  in  2000.  The  occurrence  of  similar 
rates  of  recharge  along  other  major  ephemeral  streams  would 
result  in  significant  alterations  to  the  pre-existing  conceptual 
model  of  recharge  distributions  in  the  subwatershed.  Recharge 
through  direct  infiltration  of  precipitation  through  the  basin 
floor  is  thought  to  occur  rarely  because  of  high  evaporation, 
low  precipitation,  and  water  use  by  desert  vegetation  in  the 
semiarid  Southwest  (Anderson  and  others,  1992;  Scott  and 
others,  2000).  Rare  basin-floor  recharge  could,  however,  be  a 
significant  source  of  recharge  during  extended  periods  of  high 
precipitation  and  low  evapotranspiration. 

The  occurrence  of  significant  recharge  along  the  lengths 
of  ephemeral  stream  channels  or  through  the  basin  floor, 
however,  would  not  result  in  increased  estimates  of  average 
annual  basin- wide  recharge  rates.  Water-budget  analysis 
requires  that  estimates  of  recharge  rates  be  maintained  to 
balance  inflows  and  outflows.  The  change  would  simply 
be  in  the  conceptual  model  of  recharge  distribution:  less 
concentrated  recharge  near  the  mountains  and  more  extensive 
recharge  downstream  in  tributary  channels. 

Previous  recharge  estimates  for  the  Sierra  Vista 
sub  watershed  have  generally  assumed  that  recharge  rates  are 
constant  with  time.  Trends  in  seasonal  and  annual  precipitation 
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(Pool  and  Coes,  1999),  however,  likely  result  in  long-term 
variations  in  recharge  rates.  Water-chemistry  data  for  the 
Sierra  Vista  sub  watershed  and  the  adjacent  Tucson  Basin 
indicate  that  a  greater  percentage  of  precipitation  infiltrates 
during  winter  streamflow  than  during  summer  streamflow 
(Gallaher,  1979;  Keith,  1981;  Pool  and  Coes,  1999).  Rates  of 
summer  streamflow  have  varied  randomly  from  year  to  year, 
but  rates  of  winter  streamflow  vary  over  decadal  periods  with 
Pacific  Ocean  climate  indicators  (Cayan  and  Webb,  1992; 
Webb  and  Betancourt,  1992;  Cayan  and  others,  1999;  McCabe 
and  Dettinger,  1999). 

Ground-Water  Flow  System 

Ground  water  in  the  basin-fill  aquifer  is  transmitted  from 
recharge  areas  to  discharge  areas  primarily  through  the  sands 
and  gravels  of  the  upper  and  lower  basin-fill  units.  Near  the 
mountain  fronts,  shallow  bedrock  and  relatively  impermeable 
sediments  support  perched  aquifers.  Silt  and  clay  layers 
within  the  upper  and  lower  basin-fill  units  on  the  west  side 
of  the  subwatershed  split  ground-water  flow  in  the  regional 
aquifer  into  deep  and  shallow  flow  systems  (Pool  and  Coes, 
1999).  Exceptions  to  the  regional  directions  of  ground- water 
flow  exist  near  the  Sierra  Vista-Fort  Huachuca  area,  where 
withdrawals  have  removed  ground  water  from  aquifer  storage 
and  altered  pre-development  ground-water  flow  paths  (Pool 
and  Coes,  1999). 

Post-development  discharge  from  the  basin-fill 
aquifer  in  the  Sierra  Vista  subwatershed  occurs  primarily 
as  ground-water  withdrawals,  discharge  to  the  San  Pedro 
and  Babocomari  Rivers  (base  flow),  and  evapotranspiration. 
Ground- water  withdrawals  in  the  subwatershed  began  in 
about  1940  and  peaked  in  the  early  1980s  at  about 
18.5  hm3/yr  (Freethey,  1982;  Corell  and  others,  1996). 
Withdrawals  declined  to  about  13.5  hm3/yr  by  1991  (Corell 
and  others,  1996).  Post-development  discharge  to  streams  is 
about  7.4  hm3/yr,  and  post-development  evapotranspiration 
is  estimated  to  be  about  7.0  hm3/yr  (Freethey,  1982).  About 
3.7  to  4.9  hm3/yr  discharges  from  the  sub  watershed  as 
underflow  through  basin  fill  and  stream  alluvium  in  the 
Fairbank  area  (Freethey,  1982;  Corell  and  others,  1996). 

In  addition,  a  small  amount  of  ground  water  discharges  as 
springs  near  the  San  Pedro  River. 

Pool  and  Coes  (1999)  identified  four  types  of  long-term 
water-level  changes  in  the  Sierra  Vista  subwatershed.  The  first 
type  is  periodic  decline  and  recovery  of  water  levels  near  the 
mountains  related  to  variations  in  precipitation  and  rates  of 
recharge.  The  second  is  a  decline  of  0.2  to  more  than  0.3  m/yr 
in  the  Sierra  Vista-Fort  Huachuca  area  caused  by  ground- water 
withdrawals.  The  third  is  a  regional  decline  of  0.1  to  0.2  m/yr 
during  1940  through  the  mid-1960s  or  early  1980s  followed  by 
a  period  of  no  decline  or  slight  recovery.  This  regional  decline 
could  be  caused  by  variations  in  recharge  rates  or  regional 
response  to  incision  of  the  San  Pedro  River.  The  fourth  type  of 
long-term  water-level  change  is  recovery  near  the  San  Pedro 
River  after  the  mid-1980s  that  is  probably  due  to  a  decrease  in 
agricultural  ground- water  withdrawals. 


Methods  of  Investigation 

Multiple  methods  were  used  to  estimate  infiltration 
through  ephemeral-stream  channels  and  the  basin  floor  in 
the  Sierra  Vista  subwatershed.  Information  gathered  through 
unsaturated-zone  borehole  drilling  and  monitoring  was  used 
to  better  describe  infiltration  through  the  unsaturated  zone. 
Tritium  data  were  used  to  describe  the  extent  to  which  water 
has  moved  through  the  unsaturated  zone  since  the  1950s. 
Ratios  of  the  stable  isotopes  of  oxygen  and  hydrogen  were 
used  to  identify  evaporative  and  seasonal  signatures  of  pore 
water.  The  distribution  of  chloride  in  the  subsurface  was 
used  to  identify  active  or  inactive  recharge  areas,  and  where 
possible,  was  balanced  to  determine  infiltration  fluxes  and 
water  residence  times.  The  transport  of  heat  with  infiltrating 
water  was  used  to  quantify  infiltration  rates  from  individual 
ephemeral  streamflow  events. 

Field  and  Laboratory  Methods 

To  characterize  ephemeral- stream  channel  and  basin- 
floor  recharge  in  the  Sierra  Vista  sub  watershed,  12  tributary 
sites  and  4  intertributary  sites  were  chosen  for  intensive 
unsaturated-zone  data  collection.  At  each  site,  a  borehole  was 
drilled  into  the  subsurface  and  sediment  cuttings  and  cores 
were  collected.  Of  the  16  boreholes,  15  were  completed  with 
casing  for  subsurface  monitoring.  Sediment  textures  were 
described  in  the  field,  and  sediment-chloride  concentrations 
were  measured.  Sediment  cores  were  measured  for  hydraulic, 
physical,  and  thermal  properties;  pore- water  stable-isotopes 
values;  and  tritium  activity.  Twelve  of  the  completed  boreholes 
were  instrumented  for  continuous  monitoring  of  sediment 
temperature.  Ground- water  levels  of  saturated  intervals  were 
continuously  monitored  in  three  of  the  completed  boreholes. 

Borehole  Site  Selection 

Tributary  and  intertributary  borehole  locations  were 
chosen  by  using  a  multistep  approach.  First,  the  Automated 
Geospatial  Watershed  Assessment  (AGWA)  tool  (Miller 
and  others,  2002),  an  extension  to  the  Environmental 
Systems  Research  Institute’s  Arc  View  version  3.0,  was 
used  to  delineate  the  drainage  basins  of  the  major  tributaries 
of  the  San  Pedro  River  (fig.  5).  For  each  basin,  the  main 
channel  length,  main  channel  slope,  total  basin  area,  mean 
basin  elevation,  basin  perimeter,  basin  length,  major  land 
use,  and  major  soil  type  were  quantified  using  the  Xtools 
and  Spatial  Analyst  extensions  to  Arc  View  (table  1).  The 
basins  were  then  divided  into  seven  groups  on  the  basis  of 
similar  morphometric,  land-use,  and  soil-type  characteristics 
(fig.  5).  Twelve  tributary  sites  were  chosen  in  five  drainage 
basins  that  represent  four  of  the  seven  groups:  four  sites  in 
Woodcutters  Wash,  two  sites  in  Miller  Canyon  Wash,  two 
sites  in  Greenbush  Draw,  one  site  in  Banning  Creek,  and  three 
sites  in  Walnut  Gulch.  In  addition,  four  intertributary  (basin 
floor)  sites  were  chosen  in  three  drainage  basins:  two  sites  in 
Walnut  Gulch,  one  site  in  Woodcutters  Wash,  and  one  site  in 
Carr  Canyon. 


Methods  of  Investigation  1 1 


Table  1.  Tributary  drainage  basin  geomorphic,  land-cover,  and  soil-type  characteristics.  Sierra  Vista  subwatershed,  Arizona 

[Group  numbers  correspond  with  groups  shown  on  figure  5;  km,  kilometers;  m/km,  meters  per  kilometer;  km2,  square  kilometers;  ND,  not  determined] 


Iden¬ 
tifier 
(see 
fig.  5) 

Tributary  name 

Main 

channel 

length1 

(km) 

Main 

channel 

slope 

(m/km) 

Total 

basin 

area 

(km2) 

Mean 

basin 

eleva¬ 

tion 

(m) 

Basin 

peri¬ 

meter 

(km) 

Basin 

length 

(km) 

Major  land  cover2 
(percentage) 

Major  soil  type2  (percentage) 

Group  1 

l 

Babocomari  River 

47.8 

7.9 

528.0 

1,512 

187.5 

48.2 

Grassland  (48.3) 

ND 

Group  2 

2 

Unnamed  Wash  1 

15.4 

10.9 

41.8 

1,288 

39.2 

14.4 

Desert  scrub  (44.1) 

Loam  (27.9) 

3 

Graveyard  Gulch 

14.4 

12.3 

30.9 

1,310 

35.1 

12.7 

Desert  scrub  (38.3) 

Loam  (49.7) 

4 

Coyote  Wash 

20.7 

13.8 

50.3 

1,370 

44.8 

19.4 

Urban  (37.9) 

Loam  (54.8) 

Group  3 

5 

Woodcutters  Wash 

23.0 

13.9 

37.5 

1,434 

54.7 

32.6 

Desert  scrub  (26.7) 

Loam  (34.4) 

6 

Garden  Canyon  Wash 

21.7 

13.5 

86.4 

1,651 

63.8 

33.5 

Oak  woodland  (37.0) 

Loam  (17.3) 

7 

Ramsey  Canyon  Wash 

17.4 

18.1 

32.5 

1,701 

51.5 

25.2 

Forest  (31.0) 

Clay  (22.6) 

8 

Carr  Canyon  Wash 

16.8 

18.1 

57.2 

1,590 

42.5 

23.0 

Urban  (23.1) 

Loam  (26.9) 

9 

Miller  Canyon  Wash 

16.2 

20.5 

32.8 

1,693 

48.1 

27.9 

Forest  (25.2) 

Loam  (23.3) 

10 

Hunter  Canyon  Wash 

15.5 

19.1 

24.1 

1,478 

43.0 

23.9 

Grassland  (38.5) 

Loam  (32.9) 

11 

Stump  Canyon  Wash 

14.6 

19.5 

18.5 

1,624 

36.3 

16.4 

Grassland  (35.3) 

Clay  (41.5) 

12 

Ash  Canyon  Wash 

12.9 

23.7 

17.5 

1,818 

37.5 

15.6 

Oak  woodland  (67.2) 

Paleosoil  (13.3) 

Brown-Bob 

13 

9.5 

23.3 

21.0 

1,536 

25.7 

11.0 

Grassland  (38.4) 

Loam  (31.9) 

Thompson  Wash 

Group  4 

14 

Lewis  Springs  Wash 

12.1 

12.0 

13.0 

1,300 

22.0 

12.5 

Desert  scrub  (46.7) 

Secondary  calcium  carbonate  (71.8) 

15 

Bakarich-McCool  Wash 

11.8 

12.6 

18.5 

1,320 

28.5 

13.2 

Desert  scrub  (52.0) 

Secondary  calcium  carbonate  (58.8) 

16 

Three  Canyons  Wash 

6.6 

12.7 

21.2 

1,354 

29.1 

12.7 

Grassland  (41.7) 

Clay  (85.0) 

17 

Palominas  Wash 

11.1 

19.6 

15.4 

1,385 

23.6 

17.2 

Desert  scrub  (47.3) 

Sandy  stream  channel  (52.3) 

Group  5 

18 

Greenbush  Draw 

19.5 

8.8 

306.3 

1,476 

81.6 

23.9 

Grassland  (38.6) 

ND 

Group  6 

19 

Banning  Creek 

8.6 

15.6 

41.3 

1,655 

47.8 

21.2 

Desert  scrub  (46.4) 

Secondary  calcium  carbonate  (1 1.0) 

20 

Unnamed  Wash  2 

10.1 

24.1 

19.0 

1,436 

24.7 

13.3 

Desert  scrub  (5 1 .7) 

Secondary  calcium  carbonate  (59.4) 

21 

Stagg  Ranch  Wash 

10.5 

25.7 

21.9 

1,456 

25.0 

10.2 

Desert  scrub  (63.3) 

Secondary  calcium  carbonate  (53.8) 

22 

Little  Dry  Creek 

10.6 

23.4 

19.2 

1,466 

27.6 

13.4 

Desert  scrub  (40.5) 

Secondary  calcium  carbonate  (61.1) 

23 

Spring  Creek 

11.4 

22.9 

22.6 

1,527 

33.3 

14.6 

Desert  scrub  (43.3) 

Loam  (29.0) 

Group  7 

24 

Willow  Wash 

20.9 

16.0 

130.8 

1,531 

76.6 

23.1 

Grassland  (40.1) 

Rocky  calcium  carbonate  (29.0) 

25 

Walnut  Gulch 

26.5 

12.1 

152.1 

1,411 

86.1 

34.1 

Desert  scrub  (57.6) 

Rocky  calcium  carbonate  (35.2) 

26 

Government  Draw 

17.8 

7.6 

148.4 

1,399 

80.8 

20.5 

Desert  scrub  (65.0) 

Lithic  (26.4) 

High  Knolls  Canyon 

27 

Wash 

18.7 

12.5 

98.9 

1,443 

56.5 

29.4 

Desert  scrub  (45.2) 

Lithic  (43.6) 

Trom  the  mountain  front. 

2U.S.  Environmental  Protection  Agency,  2003. 
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Borehole  Drilling  and  Completion 

Unsaturated-zone  boreholes  were  drilled  and  completed 
in  January  2001  and  November-December  2001  (table  2). 

All  boreholes  were  drilled  using  the  ODEX  air-hammer 
method,  also  known  as  the  under-reamer  method  (Driscoll, 
1986;  Hammermeister  and  others,  1986).  This  drilling  method 
minimized  disturbance  of  the  formation  near  the  borehole, 
and  enabled  collection  of  representative  cuttings  and  cores  at 
predetermined  intervals.  Boreholes  were  drilled  to  depths  of 
less  than  30.0  m  and  were  either  15.2  or  20.3  cm  in  diameter. 
At  each  borehole,  drill  cuttings  were  collected  every  0.3  m, 
and  0.6-m-long  (24  inches)  cores  were  collected  at  differing 
intervals  using  a  10.2-cm-diameter  piston  core  barrel.  The 
0.6-m-long  core  barrel  was  lined  with  three  15.2-cm-long 
sleeve  liners  and  two  7.6-cm-long  sleeve  liners.  Upon  core 
collection,  each  of  the  liners  was  immediately  capped,  taped, 
plastic  wrapped,  and  placed  in  an  aluminum  pouch,  which 
was  then  heat  sealed,  following  the  procedure  described  by 
Hammermeister  and  others  (1986),  to  preserve  the  integrity  of 
each  core  sample. 

Of  the  16  boreholes,  15  were  completed  to  enable  future 
subsurface  monitoring  (table  2).  Permission  was  not  granted 
from  the  land  owner  to  complete  the  borehole  at  WC2.  Each 
borehole  was  completed  with  5.1-cm-diameter  schedule-40 
polyvinyl  chloride  (PVC)  pipe.  Twelve  of  the  boreholes  were 
completed  with  solid- walled  PVC  and  capped  at  the  bottom 
and  top.  During  drilling,  three  of  the  boreholes  intercepted 
saturated  intervals  and  were  completed  with  solid- walled  PVC 
above  the  saturated  interval  and  a  1.5-m  section  of  slotted 
PVC  spanning  the  unsaturated/saturated  interface.  The  annular 
space  of  the  12  boreholes  completed  with  solid- walled  PVC 
was  filled  with  a  mixture  of  cuttings  from  the  boreholes  and 
sand.  The  annular  space  of  the  three  boreholes  completed 
with  slotted  PVC  was  filled  with  gravel  in  the  slotted  interval, 
bentonite  above  the  slotted  interval,  and  a  mixture  of  cuttings 
and  sand  above  the  bentonite.  Surface  casing  consists  of 
10.2-cm-diameter  schedule-80  PVC  pipe  from  above  land 
surface  to  0.9  to  1.5  m  depth  around  the  5.1-cm-diameter 
PVC  casing.  The  uppermost  0.9  to  1.5  m  of  annular  space 
between  the  surface  casing  and  the  borehole  was  filled  with 
cement.  The  surface  casing  was  capped  with  a  watertight, 
locking,  aluminum  cap.  Completion  of  the  15  boreholes  with 
PVC  casing  resulted  in  12  unsaturated-zone  monitoring  access 
points  and  3  wells.  Use  of  the  term  “borehole”  is  retained, 
however,  for  convenience  of  discussion  in  this  report. 

Sediment  Analysis 

Cuttings 

Cuttings  from  each  borehole  were  visually  described 
for  grain-size  distribution  at  0.3-m  intervals  (fig.  6).  Cuttings 
from  three  boreholes  were  also  analyzed  in  the  laboratory  for 
particle- size  distribution  (fig.  6). 


Sediment-chloride  concentrations  of  borehole  cuttings 
(except  for  borehole  WG5)  collected  at  0.3-m  intervals  were 
determined  in  the  field  (table  3).  A  subsample  of  cuttings 
representative  of  each  0.3-m  interval  was  sieved  to  obtain 
particles  less  than  2  mm  in  diameter.  Fifty  grams  of  the 
sieved  subsample  was  mixed  well  with  50  mL  of  deionized 
water  and  allowed  to  settle.  The  chloride  concentration  of  the 
leachate  was  then  measured  using  an  Orion  chloride-specific 
electrode,  Model  96-17  (Thermo  Electron  Corporation, 

2000).  This  electrode  measures  chloride  concentrations  of 
1.8  to  35,500  mg/L  with  a  reproducibility  of  ±2  percent. 

The  electrode  was  calibrated  in  the  field  at  the  beginning  of 
each  set  of  measurements;  calibration  was  checked  every  1  to 
2  hours.  Chloride  is  reported  in  milligrams  of  chloride  per 
kilogram  of  soil  (mg/kg;  MCI,  in  M/M): 

Mcl  =  Clsamp^-  ,  (1) 

where 

Cl  =  chloride  concentration  of  leachate  [M/L3], 

samp  L  J  7 

Vw  =  volume  of  deionized  water  in  leachate  [L3], 

Ms  =  mass  of  sediment  in  leachate  [M]. 

Cores 

Cores  were  analyzed  to  determine  pore- water  stable- 
isotope  values  and  tritium  activity;  thermal,  physical,  and 
hydraulic  properties;  and  particle-size  distributions.  Pore 
water  was  extracted  and  analyzed  for  the  stable  isotopes  of 
oxygen  and  hydrogen  at  the  USGS  Isotope  Fractionation 
Project  Laboratory  in  Reston,  Virginia  (table  4).  Pore 
water  was  extracted  and  analyzed  for  tritium  activity  at  the 
USGS  Water-Quality  Laboratory  in  Menlo  Park,  California 
(table  4).  The  precision  of  individual  tritium  measurements 
is  a  function  of  the  volume  of  water  extracted  from  the 
core  material  and  ranged  from  ±0.2  tritium  units  (TU)  for 
moist  cores  to  ±5.0  TU  for  an  extremely  dry  core.  Thermal 
properties  (specific  heat  and  thermal  conductivity,  at  field 
soil-water  content;  table  5)  were  measured  at  the  USGS 
Hydrologic  Research  Laboratory  in  Sacramento,  California, 
using  a  dual-needle,  line-heat  source  probe  and  a  ThermoLink 
measurement  system  (Campbell  and  others,  1991;  Decagon 
Devices,  Inc.,  1999). 


Table  2.  Borehole  construction  and  instrumentation.  Sierra  Vista  subwatershed,  Arizona 

[km,  kilometers;  m,  meters;  cm,  centimeters;  do.,  ditto;  NA,  not  applicable;  NC,  not  completed.  Surface  TidbiT:  Y,  yes;  N,  no] 


Distance 

Construction  information 

Instrumentation  information 

down¬ 

stream 

Approx- 

Bore- 

from 

imate 

hole 

Bore- 

Casing 

Depth  to 

Depth  to 

Trans¬ 

Borehole 

mountain 

channel 

diam- 

hole 

diam- 

Casing 

top  of 

bottom  of 

ducer 

Sur¬ 

Borehole 

identification 

front 

width 

End  drill 

eter 

depth 

eter 

depth 

opening 

opening 

depth 

face 

Down-hole  TidbiT  depths 

identifier 

number 

(km) 

(m) 

date 

(cm) 

(m) 

(cm) 

(m) 

Opening  type 

(m) 

(m) 

(m) 

TidbiT 

(m) 

Ephemeral-stream  channel  boreholes 

BCl 

(D-22-22)23bbd 

3.9 

4 

01-17-2001 

20.3 

26.3 

5.1 

26.3 

None 

NA 

NA 

NA 

Y 

1.5,  3.0,  4.6,  6.1,  7.6,  9.1,  12.2 

GDI 

(D-24-23)06dbd 

12.4 

14 

11-30-2001 

15.2 

19.5 

do. 

19.5 

None 

NA 

NA 

NA 

N 

1.5,  3.0,  6.1,  9.1 

GD2 

(D-24-23)14aba 

5.6 

2 

11-29-2001 

do. 

28.2 

do. 

28.2 

None 

NA 

NA 

NA 

Y 

1.5,  3.0,  4.6,  6.1,  7.6,  9.1,  12.2 

MCI 

(D-22-21)36bad 

7.1 

4 

11-18-2001 

do. 

23.0 

do. 

23.0 

None 

NA 

NA 

NA 

Y 

1.5,  3.0,  4.6,  6.1,  9.1 

MC2 

(D-23-21)bca2 

0 

2 

11-20-2001 

do. 

16.0 

do. 

16.0 

None 

NA 

NA 

NA 

N 

1.5,  3.0,  4.6,  6.1,  7.6,  9.1,  12.2 

WC1 

(D-21-21)14dbcl 

16.4 

10 

01-14-2001 

20.3 

9.4 

do. 

9.4 

None 

NA 

NA 

NA 

Y 

1.5,  3.0,  4.6,  6.1,  7.6,  9.1 

WC1 

(D-21-21)14dbc2 

do. 

do. 

do. 

do. 

do. 

do. 

8.5 

Slotted  PVC 

7.3 

8.8 

8.8 

N 

NA 

WC2 

(D-21-21)28bbb 

12.3 

22 

01-12-2001 

do. 

12.5 

NC 

NC 

NC 

NC 

NC 

NA 

N 

NA 

WC3 

(D-21-20)36caa 

7.2 

8 

11-09-2001 

15.2 

22.6 

5.1 

22.6 

None 

NA 

NA 

NA 

N 

1.5,  3.0,  4.6,  6.1,  7.6,  9.1,  12.2 

WC4a 

(D-22-20)16abcl 

.4 

0.5 

01-10-2001 

20.3 

12.5 

do. 

12.7 

None 

NA 

NA 

NA 

Y 

1.5,  3.0,  4.6,  6.1,  7.6,  9.1,  12.2 

WC4a 

(D-22-20)16abc2 

do. 

do. 

do. 

do. 

do. 

do. 

5.3 

Slotted  PVC 

3.8 

5.3 

4.6 

N 

NA 

WG1 

(D-20-2 1)01  aba 

25.6 

31 

11-11-2001 

15.2 

23.0 

do. 

23.0 

None 

NA 

NA 

NA 

N 

1.5,  3.0,  4.6,  6.1,  9.1 

WG2 

(D-20-22)01bdc 

13.2 

13 

11-16-2001 

do. 

12.2 

do. 

11.4 

None 

NA 

NA 

NA 

N 

1.5,  3.0,  4.6 

WG3 

D-19-23)33aad 

7.2 

2 

11-15-2001 

do. 

12.3 

do. 

12.3 

None 

NA 

NA 

NA 

N 

NA 

Basin-floor  boreholes 

BF1 

(D-22-21)32dad 

NA 

NA 

12-01-2001 

15.2 

15.8 

5.1 

15.1 

None 

NA 

NA 

NA 

Y 

1.5,  3.0,  4.6,  6.1,  7.6,  9.1,  12.2 

WC4c 

(D-22-20)16acb 

NA 

NA 

01-08-2001 

20.3 

12.5 

do. 

12.5 

Slotted  PVC 

4.3 

5.8 

9.1 

Y 

1.5,  3.0,  4.6,  6.1,  7.6,  9.1,  12.2 

WG4 

(D-19-22)36abb 

NA 

NA 

11-12-2001 

15.2 

13.4 

do. 

13.4 

None 

NA 

NA 

NA 

N 

NA 

WG5 

(D-19-23)36dbc 

NA 

NA 

11-13-2001 

do. 

12.3 

do. 

12.3 

None 

NA 

NA 

NA 

N 

NA 

Methods  of  Investigation 


DEPTH  BELOW  LAND  SURFACE,  IN  METERS  DEPTH  BELOW  LAND  SURFACE,  IN  METERS 


14  Ephemeral-Stream  Channel  and  Basin-Floor  Infiltration  and  Recharge  in  the  Sierra  Vista  Subwatershed,  Arizona 


A  WC4a 


FIELD 

ELECTRICAL  LITHOLOGY, 
CONDUCTIVITY,  IN 

IN  PERCENTAGE 
MILLISIEMENS  OF: 

PER  METER  GO 


NEUTRON,  ELECTRICAL  LITHOLOGY, 
IN  1,000  CONDUCTIVITY,  IN 
COUNTS  IN  PERCENTAGE 


CO 


LU 

O 

< 


oo 

□ 


D.  WC1 

NEUTRON, 
IN  1,000 
COUNTS 
PER 

SECOND 


FIELD 

LAB 

ELECTRICAL 

LITHOLOGY, 

LITHOLOGY, 

CONDUCTIVITY, 

IN 

IN 

IN 

PERCENTAGE 

PERCENTAGE 

MILLISIEMENS 

OF: 

OF: 

PER  METER 

OO  Q  > 

yd  ^  < 

OO  Q  > 

yd  ^  < 

0  5  10  15  0  50  100  □=  w  S  □=  w  § 


□ 


Figure  6.  Neutron  soil-moisture,  electrical  conductivity,  field  grain-size,  and  laboratory  grain-size  data  for  boreholes.  Sierra  Vista 
subwatershed,  Arizona.  A,  WC4a;  B,  WC3;  C,  WC2;  0,  WC1;  Er  MC2;  F,  MCI;  G,  WG3;  //,  WG2;  /,  WG1;  J,  BC1;  K  GD2;  L,  GDI; 

M,  WG4;  /V,  WG5;  0,  BF1;  P,  WC4c. 
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Figure  6.  Continued. 
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Figure  6.  Continued. 


BASIN  FILL  STREAM  ALLUVIUM 

(clay,  silt,  sand,  and  gravel)  (clay,  silt,  sand,  and  gravel) 
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Figure  6.  Continued 


BASINFUL  CALICHE  STREAM  ALLUVIUM 

(clay,  silt,  sand,  and  gravel)  II  (clay,  silt,  sand,  and  gravel) 
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Figure  6.  Continued 
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Table  3.  Sediment-chloride  concentrations  of  borehole  cuttings,  Sierra  Vista  subwatershed,  Arizona 

[cm,  centimeters;  m,  meters;  mg/kg,  milligrams  chloride  per  kilograms  sediment;  — ,  no  data;  <,  less  than] 


Depth 

to  top  of _ 

15.2  cm 

interval  - 

(m)  BC1 


Chloride 

(mg/kg) 


Ephemeral-stream  channel  boreholes 


Basin-floor  boreholes 


GDI  GD2  MCI  MC2  WC1  WC2  WC3  WC4a  WG1  WG2  WG3 


BF1  WC4c  WG4 


0.0 

.3 

.6 

— 

— 

— 

<1.8 

_ 

<1.8 

.9 

<1.8 

<1.8 

<1.8 

1.2 

<1.8 

<1.8 

<1.8 

1.5 

<1.8 

<1.8 

<1.8 

1.8 

<1.8 

<1.8 

<1.8 

2.1 

<1.8 

<1.8 

<1.8 

2.4 

1.8 

<1.8 

<1.8 

2.7 

2.4 

<1.8 

<1.8 

3.0 

<1.8 

<1.8 

<1.8 

3.4 

<1.8 

<1.8 

<1.8 

3.7 

<1.8 

<1.8 

<1.8 

4.0 

<1.8 

<1.8 

<1.8 

4.3 

1.9 

<1.8 

<1.8 

4.6 

2.7 

<1.8 

<1.8 

4.9 

2.8 

<1.8 

<1.8 

5.2 

3.4 

<1.8 

<1.8 

5.5 

3.0 

<1.8 

___ 

5.8 

3.2 

<1.8 

2.6 

6.1 

2.9 

<1.8 

<1.8 

6.4 

2.5 

<1.8 

<1.8 

6.7 

2.0 

<1.8 

<1.8 

7.0 

— 

<1.8 

<1.8 

7.3 

___ 

<1.8 

<1.8 

7.7 

<1.8 

<1.8 

<1.8 

7.9 

<1.8 

<1.8 

<1.8 

8.2 

<1.8 

<1.8 

<1.8 

8.5 

<1.8 

<1.8 

___ 

8.8 

<1.8 

<1.8 

<1.8 

9.1 

<1.8 

<1.8 

___ 

9.4 

<1.8 

<1.8 

<1.8 

9.8 

<1.8 

<1.8 

<1.8 

10.1 

<1.8 

<1.8 

<1.8 

10.4 

<1.8 

<1.8 

<1.8 

10.7 

3.2 

<1.8 

<1.8 

11.0 

<1.8 

<1.8 

<1.8 

11.3 

2.8 

<1.8 

<1.8 

11.6 

2.5 

<1.8 

<1.8 

11.9 

5.2 

<1.8 

<1.8 

12.2 

8.4 

<1.8 

<1.8 

12.5 

5.1 

<1.8 

<1.8 

12.8 

5.9 

<1.8 

<1.8 

13.1 

3.5 

<1.8 

— 

13.4 

3.2 

<1.8 

___ 

13.7 

5.1 

<1.8 

___ 

14.0 

4.1 

<1.8 

___ 

14.3 

3.2 

<1.8 

<1.8 

14.6 

2.6 

<1.8 

<1.8 

14.9 

2.3 

<1.8 

<1.8 

15.2 

5.2 

<1.8 

<1.8 

15.5 

3.6 

<1.8 

<1.8 

15.8 

2.7 

<1.8 

<1.8 

16.2 

3.6 

<1.8 

<1.8 

16.5 

2.4 

<1.8 

<1.8 

16.8 

2.5 

<1.8 

<1.8 

17.1 

___ 

<1.8 

<1.8 

— 

1.9 

— 

— 

<1.8 

<1.8 

— 

___ 

<1.8 

<1.8 

<1.8 

3.3 

<1.8 

<1.8 

<1.8 

6.5 

2.1 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

5.0 

2.2 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

3.6 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

3.7 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

2.6 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

2.9 

<1.8 

1.8 

___ 

<1.8 

<1.8 

1.9 

— 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

3.9 

— 

<1.8 

<1.8 

2.3 

___ 

<1.8 

<1.8 

1.9 

___ 

<1.8 

<1.8 

4.2 

___ 

___ 

<1.8 

3.4 

___ 

___ 

<1.8 

3.2 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

<1.8 

1.9 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

<1.8 

1.9 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

<1.8 

___ 

___ 

___ 

<1.8 

___ 

___ 

___ 

<1.8 

___ 

___ 

___ 

<1.8 

___ 

— 

___ 

9.9 

___ 

— 

___ 

7.9 

3.5 

___ 

___ 

8.5 

3.6 

___ 

<1.8 

27.3 

4.6 

___ 

<1.8 

9.1 

7.5 

___ 

<1.8 

8.9 

5.1 

— 

___ 

9.3 

7.3 

___ 

<1.8 

36.9 

3.9 

<1.8 

<1.8 

32.6 

<1.8 

___ 

<1.8 

73.1 

<1.8 

2.6 

___ 

7.8 

<1.8 

<1.8 

— 

3.0 

<1.8 

<1.8 

___ 

2.7 

<1.8 

<1.8 

___ 

24.2 

<1.8 

<1.8 

___ 

<1.8 

2.0 

<1.8 

___ 

<1.8 

1.9 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

— 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

3.3 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

___ 

___ 

<1.8 

<1.8 

<1.8 

— 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

___ 

___ 

<1.8 

___ 

<1.8 

___ 

2.0 

2.1 

___ 

___ 

2.9 

<1.8 

<1.8 

___ 

2.2 

2.3 

<1.8 

___ 

2.7 

2.1 

<1.8 

— 

4.1 

<1.8 

<1.8 

___ 

<1.8 

<1.8 

___ 

___ 

<1.8 

<1.8 

<1.8 

— 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

— 

<1.8 

___ 

<1.8 

___ 

<1.8 

— 

<1.8 

___ 

<1.8 

— 

<1.8 

— 

<1.8 

— 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

— 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

<1.8 

___ 

___ 

3.1 

3.0 

— 

___ 

<1.8 

2.8 

___ 

___ 

3.3 

6.1 

— 

<1.8 

2.4 

5.2 

42.2 

<1.8 

3.3 

17.1 

130.0 

3.8 

<1.8 

___ 

222.0 

___ 

2.0 

11.1 

226.0 

___ 

<1.8 

8.1 

168.0 

___ 

<1.8 

8.2 

73.6 

___ 

1.9 

3.8 

54.6 

___ 

2.5 

2.4 

34.2 

___ 

2.1 

2.5 

90.3 

___ 

<1.8 

2.4 

77.5 

___ 

<1.8 

3.3 

66.1 

___ 

<1.8 

2.8 

54.8 

— 

1.8 

3.2 

50.3 

___ 

2.3 

___ 

61.7 

___ 

3.9 

2.4 

33.2 

___ 

4.8 

3.0 

15.0 

___ 

10.1 

6.0 

20.0 

___ 

16.4 

4.4 

9.1 

___ 

15.8 

2.5 

7.3 

___ 

21.2 

4.8 

6.7 

___ 

22.1 

4.4 

7.7 

___ 

27.7 

<1.8 

20.1 

___ 

22.8 

<1.8 

21.1 

___ 

29.3 

<1.8 

22.6 

___ 

24.2 

2.7 

23.4 

___ 

25.2 

3.2 

34.3 

___ 

20.2 

<1.8 

40.7 

___ 

15.0 

<1.8 

32.5 

___ 

21.7 

<1.8 

16.1 

___ 

30.4 

<1.8 

18.8 

___ 

21.2 

<1.8 

22.7 

___ 

19.3 

<1.8 

3.2 

___ 

23.8 

<1.8 

4.3 

— 

28.1 

<1.8 

5.8 

___ 

32.1 

<1.8 

5.6 

___ 

26.4 

<1.8 

4.3 

___ 

22.6 

<1.8 

4.3 

___ 

23.5 

<1.8 

4.1 

___ 

20.4 

___ 

3.7 

___ 

17.1 

___ 

3.6 

___ 

11.8 

___ 

___ 

___ 

10.7 

— 

— 

___ 

14.3 

___ 

___ 

___ 

15.5 

— 

— 

___ 

9.9 

___ 

___ 

___ 

7.1 

___ 

___ 

___ 

7.1 

— 

— 

___ 

6.0 

___ 

___ 

___ 

5.5 

— 

___ 
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Table  3.  Sediment-chloride  concentrations  of  borehole  cuttings,  Sierra  Vista  subwatershed,  Arizona — Continued 


Depth 
to  top  of 

Chloride 

(mg/kg) 

15.2  cm 
interval 
(m) 

Ephemeral-stream  channel  boreholes 

Basin-floor  boreholes 

BC1 

GDI 

GD2 

MCI 

MC2 

WC1 

WC2 

WC3  WC4a 

WG1 

WG2 

WG3 

BF1 

WC4c 

WG4 

17.4 

4.9 

<1.8 

<1.8 

<1.8 

___ 

— 

___ 

<1.8 

<1.8 

___ 

— 

___ 

___ 

___ 

17.7 

6.6 

<1.8 

<1.8 

<1.8 

___ 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

___ 

18.0 

4.5 

<1.8 

<1.8 

<1.8 

— 

— 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

___ 

18.3 

8.1 

<1.8 

<1.8 

<1.8 

___ 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

___ 

18.6 

10.5 

<1.8 

<1.8 

<1.8 

— 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

— 

18.9 

11.1 

<1.8 

<1.8 

<1.8 

___ 

— 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

___ 

19.2 

11.9 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

___ 

19.5 

9.6 

— 

<1.8 

<1.8 

___ 

___ 

— 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

___ 

19.8 

8.2 

___ 

<1.8 

<1.8 

___ 

___ 

— 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

— 

20.1 

2.4 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

<1.8 

<1.8 

— 

___ 

___ 

___ 

___ 

20.4 

<1.8 

— 

<1.8 

<1.8 

___ 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

— 

20.7 

<1.8 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

___ 

21.0 

2.5 

— 

<1.8 

<1.8 

___ 

___ 

___ 

<1.8 

<1.8 

___ 

— 

___ 

___ 

___ 

21.3 

2.0 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

— 

21.6 

3.1 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

<1.8 

<1.8 

— 

___ 

___ 

___ 

___ 

21.9 

6.8 

___ 

<1.8 

<1.8 

___ 

— 

— 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

___ 

22.3 

9.4 

___ 

<1.8 

___ 

— 

___ 

___ 

<1.8 

<1.8 

___ 

___ 

___ 

___ 

___ 

22.6 

4.8 

— 

<1.8 

___ 

___ 

___ 

___ 

<1.8 

___ 

___ 

___ 

___ 

___ 

___ 

22.9 

5.0 

— 

<1.8 

___ 

___ 

___ 

___ 

___ 

___ 

___ 

___ 

___ 

___ 

23.2 

7.2 

___ 

<1.8 

___ 

___ 

___ 

___ 

— 

___ 

___ 

___ 

— 

___ 

___ 

23.5 

6.2 

___ 

<1.8 

___ 

___ 

— 

___ 

___ 

___ 

— 

___ 

___ 

— 

___ 

23.8 

3.5 

___ 

<1.8 

___ 

___ 

___ 

___ 

— 

___ 

___ 

___ 

___ 

___ 

___ 

24.1 

4.5 

— 

<1.8 

— 

— 

— 

___ 

— 

___ 

___ 

___ 

___ 

___ 

___ 

24.4 

3.1 

___ 

<1.8 

___ 

___ 

___ 

___ 

___ 

— 

___ 

___ 

___ 

___ 

___ 

24.7 

5.1 

___ 

<1.8 

___ 

— 

___ 

___ 

___ 

___ 

___ 

___ 

___ 

___ 

— 

25.0 

25.3 

3.5 

4.5 

_ 

<1.8 

_ 

— 

— 

— 

— 

_ 

_ 

_ 

— 

_ 

_ 

25.6 

25.9 

3.7 

3.3 
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Table  4.  Stable-isotope  and  tritium  data  for  pore  water  samples  from  cores.  Sierra  Vista  subwatershed,  Arizona 

[cm,  centimeters;  m,  meters;  %c,  per  mil;  TU,  tritium  units;  — ,  no  data] 


Depth  to  top  of 
15.2-cm  interval 
(m) 

Delta  deuterium 

(%.) 

Delta  oxygen-18 

(%.) 

Tritium 

(TU) 

Tritium  error 
(±TU) 

Ephemeral-stream  channel  boreholes 

BC1 

2.0 

-63.5 

-9.50 

3.8 

0.2 

3.5 

-64.0 

-9.22 

3.3 

.2 

5.0 

-60.8 

-8.80 

3.2 

.6 

6.6 

-61.8 

-8.96 

3.8 

.6 

8.1 

-60.8 

-9.23 

4.8 

.3 

9.6 

-60.0 

-8.52 

4.4 

.3 

11.1 

-63.4 

-9.25 

4.1 

.3 

12.6 

-60.0 

-8.81 

4.5 

.7 

14.2 

-63.8 

-9.24 

3.4 

.3 

15.7 

-63.6 

-9.00 

3.9 

.3 

17.2 

-59.5 

-8.51 

6.7 

.9 

20.3 

-57.0 

-8.30 

3.6 

.6 

23.3 

-60.0 

-8.65 

4.5 

.9 

GDI 

.5 

-60.0 

-.41 

4.0 

5.0 

.8 

-65.6 

-9.64 

3.7 

.4 

5.3 

-64.7 

-9.67 

4.0 

.3 

6.9 

-65.1 

-9.41 

3.3 

.3 

8.4 

-64.4 

-9.18 

4.6 

.6 

9.9 

-62.0 

-9.09 

4.3 

.3 

19.1 

-65.6 

-9.68 

4.1 

.3 

GD2 

2.3 

-55.5 

-8.32 

3.5 

.3 

3.8 

-58.4 

-8.57 

3.7 

.2 

5.5 

-59.3 

-8.24 

4.1 

.3 

6.9 

-62.0 

-8.89 

4.3 

.3 

8.4 

-64.5 

-9.02 

3.9 

.3 

9.9 

-67.9 

-9.34 

3.8 

.3 

13.0 

-63.5 

-8.83 

4.0 

.4 

16.0 

-63.7 

-9.02 

4.3 

.3 

19.1 

-67.0 

-9.03 

3.8 

.4 

MCI 

.2 

-39.4 

-1.25 

7.8 

.6 

1.1 

-59.8 

-8.59 

4.2 

.3 

2.3 

-59.5 

-8.44 

4.7 

.3 

3.8 

-65.8 

-9.07 

4.2 

.3 

5.3 

-69.3 

-9.23 

4.1 

.3 

6.9 

-63.1 

-8.96 

4.0 

.3 

8.4 

-64.9 

-8.94 

4.4 

.3 

9.9 

-64.9 

-9.10 

4.0 

.3 

13.0 

-72.5 

-10.09 

3.7 

.2 

16.0 

-66.0 

-9.00 

3.6 

.3 

19.1 

-64.9 

-8.81 

4.1 

.4 

22.1 

-84.6 

-11.77 

3.3 

.3 

MC2 

.2 

-65.1 

-5.28 

___ 

— 

1.7 

-64.7 

-9.14 

1.9 

.5 

6.9 

-71.3 

-9.93 

3.8 

.3 

9.9 

-69.1 

-9.71 

5.9 

.4 
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Table  4.  Stable-isotope  and  tritium  data  for  pore  water  samples  from  cores.  Sierra  Vista  subwatershed,  Arizona — Continued 


Depth  to  top  of 

Delta 

15.2-cm  interval 

deuterium 

Delta  oxygen-18 

Tritium 

Tritium  error 

(m) 

(%o) 

(%o) 

(TU) 

(TU) 

WC1 

0.2 

-63.5 

-7.70 

7.2 

0.8 

2.0 

-78.7 

-11.60 

5.9 

.5 

3.5 

-57.1 

-8.02 

4.3 

.3 

5.0 

-61.7 

-8.80 

4.4 

.2 

6.6 

-57.8 

-8.40 

4.6 

.2 

8.1 

-57.8 

-8.30 

3.9 

.2 

WC2 

2.0 

-59.2 

-6.29 

3.9 

.7 

3.5 

-60.0 

-7.70 

5.4 

.6 

5.0 

-59.5 

-7.70 

5.0 

.6 

6.6 

-46.8 

-6.70 

4.6 

.6 

8.1 

-50.5 

-6.90 

5.7 

.6 

9.6 

-50.0 

-6.80 

4.3 

.5 

11.1 

-63.8 

-8.43 

3.9 

.3 

WC3 

.2 

-36.1 

-1.36 

6.3 

.5 

.8 

-48.8 

-5.83 

4.5 

.4 

1.4 

-48.6 

-6.10 

5.0 

.4 

2.9 

-55.9 

-7.05 

3.8 

.3 

8.1 

-52.2 

-6.90 

5.5 

.3 

10.5 

-49.4 

-6.75 

5.4 

.5 

13.6 

-59.9 

-7.84 

2.6 

.3 

16.6 

-62.5 

-8.26 

1.0 

.6 

19.7 

-62.1 

-8.04 

0.4 

.3 

WC4a 

.8 

-63.0 

-7.70 

3.7 

.7 

2.1 

-59.0 

-7.38 

3.4 

1.1 

3.6 

-64.6 

-8.82 

5.0 

.4 

5.1 

-60.0 

-8.20 

10.1 

.7 

6.6 

-66.0 

-8.90 

1.6 

.5 

8.2 

-63.0 

-8.30 

1.0 

.4 

9.6 

-73.0 

-9.90 

1.6 

.5 

11.1 

-68.8 

-9.55 

1.0 

.4 

WG1 

2.9 

-50.9 

-6.94 

2.9 

1.2 

3.5 

-67.7 

-9.60 

3.6 

.7 

7.5 

-64.7 

-9.14 

3.9 

.4 

8.7 

-62.0 

-9.06 

4.5 

.4 

10.2 

-63.3 

-9.24 

4.4 

.4 

13.3 

-52.3 

-7.60 

3.9 

.8 

16.3 

-55.5 

-8.10 

3.7 

.4 

19.4 

-54.6 

-8.14 

4.2 

.3 

22.4 

-59.6 

-8.66 

5.7 

.5 

WG2 

.5 

-29.9 

-2.80 

5.6 

1.1 

.8 

-28.5 

-3.72 

6.7 

.6 

2.3 

-55.3 

-8.38 

4.7 

.3 

WG3 

.2 

-32.0 

3.27 

8.5 

2.4 

.8 

-29.1 

3.32 

6.5 

1.5 
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Table  4.  Stable-isotope  and  tritium  data  for  pore  water  samples  from  cores.  Sierra  Vista  subwatershed,  Arizona — Continued 


Depth  to  top  of  15.2-cm 

Delta 

interval 

deuterium 

Delta  oxygen-18 

Tritium 

Tritium  error 

(m) 

(%o) 

(%o) 

(TU) 

(TU) 

Basin-floor  boreholes 

BF1 

0.5 

-45.8 

-0.60 

6.4 

0.9 

3.8 

-62.2 

-7.82 

6.4 

.9 

5.3 

-59.5 

-7.34 

2.5 

.6 

6.9 

-57.2 

-7.03 

.8 

.6 

8.4 

-59.8 

-7.11 

.8 

.5 

10.2 

-58.0 

-7.18 

.6 

.4 

13.0 

-56.5 

-7.04 

.3 

.2 

WC4c 

.2 

-67.0 

-9.20 

4.1 

.7 

.8 

-75.5 

-10.37 

2.6 

.8 

1.7 

-68.0 

-9.40 

4.6 

1.0 

2.3 

-68.4 

-8.99 

7.0 

2.0 

3.7 

-77.0 

-9.60 

4.9 

1.4 

5.2 

-64.4 

-8.07 

3.7 

.4 

6.7 

-59.8 

-7.74 

5.6 

.3 

8.2 

-63.0 

-8.56 

11.4 

.4 

9.8 

-63.0 

-8.68 

10.4 

.4 

11.3 

-71.0 

-9.40 

.3 

.3 

WG4 

.5 

-44.1 

-0.32 

6.3 

.7 

1.1 

-42.5 

-0.94 

9.1 

.9 

2.0 

-52.2 

-2.82 

10.2 

1.2 

4.1 

-53.8 

-3.91 

8.2 

.4 

7.5 

-43.1 

-3.37 

4.0 

2.4 

13.3 

-55.9 

-5.23 

.4 

.6 

WG5 

.2 

-60.0 

-6.12 

.4 

.4 
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Table  5.  Thermal  properties  of  core  material,  Sierra  Vista  subwatershed,  Arizona 

[cm,  centimeters;  m,  meters;  cm3/cm3,  cubic  centimeter  per  cubic  centimeter;  J/m3-°C,  joules  per  cubic  meter  degrees  Celsius;  J/s-m-°C,  joules  per  second 
meter  degrees  Celsius;  — ,  no  data] 


Volumetric  Volumetric 


Depth  to  top 
of  15.2-cm 
interval 
(m) 

Field  water 

content 

(cm3/cm3) 

specific 
heat,  field 
water  content 
(J/m3-°C-106) 

Thermal  con¬ 
ductivity,  field 
water  content 
(J/s-m-C) 

Depth  to  top 
of  15.2-cm 
interval 
(m) 

Field  water 

content 

(cm3/cm3) 

specific 
heat,  field 
water  content 
(J/m3-°C-106) 

Thermal  con¬ 
ductivity,  field 
water  content 
(J/sm-C) 

Ephemeral-stream  channel  boreholes 

BC1 

WC4a 

2.3 

0.20 

2.06 

2.25 

1.1 

___ 

1.71 

1.38 

3.8 

___ 

2.24 

2.03 

2.4 

.25 

1.36 

.79 

5.3 

.09 

1.94 

2.87 

3.9 

— 

1.63 

1.64 

6.9 

— 

1.68 

2.32 

5.4 

— 

2.98 

1.66 

8.4 

___ 

3.93 

2.30 

6.9 

— 

3.25 

2.94 

9.9 

— 

2.78 

4.11 

8.4 

— 

3.33 

2.32 

11.4 

— 

3.46 

2.11 

9.9 

___ 

1.94 

.90 

14.5 

.27 

2.15 

2.59 

11.4 

.20 

1.43 

1.91 

16.0 

.29 

1.07 

1.76 

WG1 

17.5 

— 

2.54 

1.50 

2.7 

0.25 

1.63 

0.86 

20.6 

.08 

2.32 

2.61 

10.1 

— 

1.41 

.93 

23.6 

___ 

1.34 

1.74 

19.2 

___ 

2.40 

1.55 

GDI 

WG2 

6.7 

.11 

1.74 

1.56 

2.1 

___ 

1.23 

.80 

18.9 

___ 

2.42 

1.63 

WG3 

GD2 

.6 

___ 

1.15 

.27 

2.1 

.44 

2.54 

1.66 

Basin-floor  boreholes 

9.8 

___ 

2.30 

1.63 

BF1 

18.9 

___ 

2.58 

.77 

6.7 

.11 

2.47 

1.73 

MCI 

12.8 

___ 

2.42 

1.45 

0 

— 

1.34 

1.09 

WC4c 

8.2 

.28 

3.01 

1.99 

.5 

___ 

1.22 

.56 

18.9 

___ 

2.76 

1.54 

1.1 

.09 

1.25 

.53 

MC2 

2.0 

___ 

1.48 

.69 

6.7 

.30 

1.90 

2.08 

2.6 

.07 

9.70 

.42 

WC1 

4.0 

.15 

1.55 

.67 

.5 

.13 

1.37 

.95 

5.5 

___ 

2.67 

1.49 

2.3 

.07 

1.11 

.60 

7.0 

— 

3.20 

2.31 

3.8 

___ 

5.08 

1.32 

8.5 

— 

2.42 

2.05 

5.3 

___ 

4.41 

1.22 

10.1 

___ 

2.00 

1.89 

6.9 

.37 

3.41 

2.00 

11.6 

.24 

1.44 

1.12 

8.4 

___ 

2.60 

1.57 

WG4 

WC2 

0 

___ 

1.49 

.31 

2.3 

.15 

1.34 

1.20 

13.1 

.10 

1.33 

.89 

3.8 

.33 

1.75 

1.00 

WG5 

5.3 

___ 

2.71 

1.93 

0 

.10 

1.24 

.34 

6.9 

.47 

2.79 

1.33 

8.4 

.40 

1.27 

1.50 

9.9 

.18 

2.22 

1.50 

11.4 

___ 

3.02 

1.63 

WC3 

0 

___ 

2.72 

1.35 

10.4 

— 

1.23 

1.03 
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Core  analyses  for  hydraulic  properties  (saturated 
hydraulic  conductivity  and  van  Genuchten  fitting  parameters) 
and  physical  properties  (bulk  density,  volumetric  soil-water 
content,  porosity,  and  particle  density)  were  completed  in  a 
laboratory  at  the  University  of  Arizona  in  Tucson  (table  6). 
Initially,  each  core’s  weight  was  measured  and  its  volume 
calculated.  Cores  were  then  fully  saturated  and  weighed  again. 
Saturated  hydraulic  conductivity  (Ks,  in  L/t)  was  measured  by 
the  falling-head  method  (Klute  and  Dirksen,  1986).  Reported 
Ks  values  are  the  average  of  two  to  five  falling-head  tests  run 
sequentially.  Van  Genuchten  fitting  parameters  alpha  (a),  in 
1/L;  n ,  dimensionless;  residual  soil- water  content  (0R),  in 
L3/L3;  saturated  soil- water  content  (0S),  in  L3/L3)  were 
determined  from  moisture-release  data  measured  by  the 
pressure-step  outflow  method  (Gardner,  1956)  from  about 
15  to  15,000  cm  water  pressure.  Cores  were  then  oven  dried 
at  105°C  for  24  hours,  and  weighed  a  third  time.  Dry  bulk 
density  (p  ,  in  M/L3)  was  determined  by  dividing  the  weight 
of  the  core  after  oven  drying  by  the  calculated  volume  of  the 
core.  Gravimetric  field  soil-water  content  (0  ,  in  M/M)  was 
calculated  as  the  difference  between  the  initial  core  weight 
and  the  weight  of  the  core  after  oven  drying,  divided  by  the 
weight  of  the  core  after  oven  drying.  Volumetric  field  soil- 
water  content  (0 ,  in  L3/L3)  was  calculated  as  the  product  of 
the  gravimetric  field  soil- water  content  and  the  bulk  density. 
Porosity  (d>,  dimensionless)  was  assumed  equal  to  the 
saturated  soil-water  content  and  was  calculated  as  the  product 
of  the  bulk  density  and  the  difference  between  the  saturated 
core  weight  and  the  dry  core  weight,  divided  by  the  dry  core 
weight.  Particle  density  (pp,  in  M/L3)  was  calculated  by: 

_  pb 

Pp  i  —  (j)  (Danielson  and  Sutherland,  1986).  (2) 

The  Soil,  Water,  and  Plant  Analysis  Laboratory  at  the 
University  of  Arizona  measured  particle-size  distributions  for 
the  cores  (table  6).  Gravel-size  fractions  (greater  than  2  mm) 
were  determined  by  sieve  analysis.  Sand-,  silt-,  and  clay- sized 
fractions  were  determined  by  hydrometer  analysis  (Gee  and 
Bauder,  1986). 

Borehole  Logging 

Electromagnetic-induction  (EMI)  and  neutron  soil- 
moisture  geophysical  borehole  tools  were  used  to  augment 
subsurface  sampling  of  borehole  cuttings  and  cores  (fig.  6). 
EMI  is  used  to  measure  the  electrical  conductivity  of  the 
sediments  near  a  borehole  by  measuring  the  strength  of  a 
secondary  magnetic  field  induced  in  the  sediments  by  a 
primary  electromagnetic  field.  The  EMI  tool  has  a  variable 
radius  of  investigation  dependent  on  the  frequency  of  the 
primary  field,  the  spacing  between  the  primary  and  secondary 
magnetic  coils,  and  the  subsurface  conductivity.  A  frequency 
of  about  40  kilohertz  and  a  coil  spacing  of  0.5  m  was  used 
for  this  investigation.  At  this  EMI  setting,  about  90  percent 
of  the  secondary  magnetic  field  results  from  materials 


within  a  radius  of  about  1.3  m  from  the  center  of  the  coil 
spacing  (Peterson,  1999).  Variations  in  subsurface  electrical 
conductivity  are  primarily  caused  by  variations  in  amounts  of 
electrically  conductive  silt,  clay,  and  dissolved  ions  in  water. 
Electrical  conductivity  is  proportional  to  silt,  clay,  and  soil- 
water  contents.  Dry  zones  with  little  silt  and  clay  can  have 
electrical  conductivity  values  well  below  10  millisiemens 
per  meter  (mS/m),  but  saturated  clay  can  have  values  greater 
than  100  mS/m.  Variations  in  electrical  conductivity  can 
also  be  related  to  variations  in  electrically  conductive  salts; 
however,  the  salt  contents  in  most  of  the  boreholes  drilled  for 
this  investigation  were  very  low  and  did  not  likely  contribute 
significantly  to  variations  in  electrical  conductivity. 

The  neutron  tool  used  in  this  study  contained  a 
50-millicurie  americium-241  /beryllium  neutron  source 
and  detector  to  measure  variations  in  the  amount  of  water 
molecules  near  the  borehole.  The  tool  had  a  radius  of 
investigation  of  about  60  to  75  cm.  Response  of  the  neutron 
tool  is  reported  as  measured  neutron  counts  at  the  detector; 
higher  neutron  counts  indicate  greater  amounts  of  water  near 
the  borehole.  Neutron  counts  can  be  converted  to  soil- water 
content  by  using  a  soil-dependent  calibration;  however, 
calibration  was  not  attempted  for  this  investigation,  and  the 
resultant  neutron  counts  from  different  boreholes  can  only  be 
compared  in  a  relative  sense. 

Borehole  Monitoring 

Onset  Stow  Away  TidbiT  temperature  sensors  were 
installed  in  boreholes  WC1,  WC4a,  WC4c,  and  BC1  in 
February  2001,  and  in  boreholes  WC3,  MCI,  MC2,  BF1, 

GDI,  GD2,  WG1,  and  WG2  from  January  to  May  2002. 
Sensors  were  installed  from  land  surface  to  a  maximum  depth 
of  12.2  m  at  intervals  of  1.5  to  3.0  m  (table  2).  Additionally,  at 
most  boreholes,  a  sensor  was  buried  about  15  cm  below  land 
surface  next  to  the  borehole  surface  casing.  Temperature  was 
measured  and  recorded  every  30  minutes.  The  TidbiT  sensors 
measure  temperatures  of  -5°C  to  +37°C  with  an  accuracy 
of  ±0.2°C. 

In-Situ  Inc.  miniTROLL  pressure  transducers  were 
installed  in  the  three  boreholes  that  intercepted  saturated 
intervals  (WC1,  WC4a,  and  WC4c)  during  February  2001. 
Depth  to  water  was  measured  and  recorded  every  30  minutes. 

Precipitation  Collection  and  Analysis 

Composite  atmospheric  deposition  samples  were 
collected  twice  a  year  after  the  summer  (April  to  October) 
and  winter  (October  to  April)  rains  at  the  Upper  Babocomari 
River  and  Greenbush  Draw  streamflow-gaging  stations 
(fig.  1).  Sample  collection  began  in  July  2000  and  continued 
through  April  2003.  Samples  were  analyzed  for  bulk 
chloride  concentrations,  stable-isotope  values,  and  tritium 
activity  (table  7).  Bulk  chloride  was  determined  by  ion 
chromatography  at  the  USGS  laboratory  in  San  Diego, 
California.  The  stable  isotopes  of  oxygen  and  hydrogen 
were  analyzed  at  the  USGS  Isotope  Fractionation  Project 
Laboratory  in  Reston,  Virginia.  Tritium  activity  was 
determined  at  the  Laboratory  of  Isotope  Geochemistry  at  the 
University  of  Arizona  in  Tucson. 


Table  6.  Physical  and  hydraulic  properties  of  core  material,  Sierra  Vista  subwatershed,  Arizona  c? 

[cm,  centimeters;  m,  meters;  %,  percent;  g/cm3,  grams  per  cubic  centimeter;  cm3/cm3,  cubic  centimeter  per  cubic  centimeter;  cm/s,  centimeters  per  second;  1/cm,  inverse  centimeters;  — ,  no  data. 

Cores  of  basin  fill  material  could  not  be  collected  at  WGZ  and  WG3,  and  properties  were  not  determined  for  the  stream  alluvium  at  these  boreholes] 


van  Genuchten  fitting  parameters 


Depth 
to  top  of 
7.6-cm 
interval 
(m) 

Gravel 

(%) 

Sand 

<%) 

Silt 

(%) 

Clay 

(%) 

Bulk 

density 

(g/cm3) 

Particle 

density 

(g/cm3) 

Porosity 

(cm3/cm3) 

Volumetric 
water  content, 
field 

(cm3/cm3) 

Saturated 

hydraulic 

conductivity 

(cm/s) 

a 

(1/cm) 

n 

(dimen¬ 

sionless) 

m 

(1-1/n) 

Residual 

water 

content 

(cm3/cm3) 

Saturated 

water 

content 

(cm3/cm3) 

R2 

Ephemeral-stream  channel  boreholes 

BC1 

2.1 

63.7 

22.6 

8.9 

4.8 

1.85 

2.81 

0.34 

0.20 

9.6- 103 

0.145 

1.86 

0.46 

0.19 

0.34 

0.99 

5.2 

79.7 

16.2 

2.6 

1.5 

1.82 

2.65 

.31 

.09 

1.2-10 1 

.122 

2.01 

.50 

.11 

.31 

1.00 

8.2 

— 

— 

— 

— 

— 

— 

— 

— 

4.  MO'7 

— 

— 

— 

— 

— 

— 

14.3 

70.8 

19.9 

7.4 

1.9 

1.82 

3.00 

.39 

.27 

1.7  106 

.138 

1.45 

.31 

.07 

.39 

1.00 

15.8 

35.3 

37.8 

23.5 

3.4 

1.63 

2.88 

.43 

.29 

3.3- 10'5 

.009 

1.62 

.38 

.16 

.42 

1.00 

20.4 

77.7 

19.5 

1.9 

.9 

2.02 

2.70 

.25 

.08 

8.9- 10-3 

.145 

1.96 

.49 

.08 

.25 

1.00 

GDI 

4.0 

47.6 

48.3 

2.2 

1.9 

1.82 

2.49 

.27 

.15 

1.4-10'2 

.145 

2.19 

.54 

.11 

.26 

.95 

7.0 

46.1 

47.6 

2.8 

3.5 

1.89 

— 

— 

.11 

4.6- 10-3 

.145 

1.29 

.23 

.06 

.36 

1.00 

13.1 

19.4 

38.9 

27.3 

14.4 

— 

— 

— 

— 

6.8- 10  s 

— 

— 

— 

— 

— 

— 

GD2 

2.4 

2.7 

66.1 

21.9 

9.3 

1.47 

3.11 

.53 

.44 

5.9- 10-3 

.067 

1.11 

.10 

.20 

.50 

.99 

5.3 

50.7 

33.0 

4.2 

12.1 

1.44 

2.60 

.45 

.14 

5.3- 10-2 

.145 

2.54 

.61 

.18 

.44 

.97 

13.1 

70.9 

12.5 

4.2 

12.4 

1.57 

2.64 

.41 

.27 

1.6-102 

.145 

2.40 

.58 

.25 

.40 

.97 

MCI 

0.9 

4.2 

60.1 

22.9 

12.8 

1.84 

2.87 

.36 

— 

9.3- 10-3 

.145 

1.20 

.17 

.11 

.36 

.99 

4.0 

18.9 

44.6 

19.1 

17.4 

1.84 

2.91 

.37 

.31 

3.5- 10-6 

.005 

1.33 

.25 

.12 

.35 

.95 

8.5 

4.2 

57.6 

25.3 

12.9 

1.78 

2.92 

.39 

.28 

1.3-10 3 

.005 

1.21 

.17 

.12 

.39 

.95 

19.2 

4.5 

31.5 

32.1 

31.9 

— 

— 

— 

— 

1.8-106 

— 

— 

— 

— 

— 

— 

MC2 

1.8 

36.4 

42.7 

12.2 

8.7 

1.37 

2.35 

.42 

.18 

3.3-10-3 

.066 

1.22 

.19 

.08 

.41 

.99 

7.0 

3.7 

63.7 

21.0 

11.6 

1.87 

2.92 

.36 

.30 

3.2-106 

.013 

1.37 

.27 

.16 

.36 

.95 

10.1 

73.0 

19.0 

4.9 

3.1 

1.60 

2.76 

.42 

.22 

1.7  102 

.136 

1.33 

.25 

.14 

.42 

.97 

Ephemeral-Stream  Channel  and  Basin-Floor  Infiltration  and  Recharge  in  the  Sierra  Vista  Subwatershed,  Arizona 


Table  6.  Physical  and  hydraulic  properties  of  core  material.  Sierra  Vista  subwatershed,  Arizona — Continued 


van  Genuchten  fitting  parameters 

Depth 


to  top  of 
7.6-cm 
interval 
(m) 

Gravel 

(%) 

Sand 

(%) 

Silt 

(%) 

Clay 

(%) 

Bulk 

density 

(g/cm3) 

Particle 

density 

(g/cm3) 

Porosity 

(cm3/cm3) 

Volumetric 
water  content, 
field 

(cm3/cm3) 

Saturated 

hydraulic 

conductivity 

(cm/s) 

a 

(1/cm) 

n 

(dimen¬ 

sionless) 

m 

(1-1/n) 

Residual 

water 

content 

(cm3/cm3) 

Saturated 

water 

content 

(cm3/cm3) 

R2 

Ephemeral-stream  channel  boreholes 

WC1 

0.3 

25.6 

72.5 

<0.7 

1.9 

1.49 

2.74 

0.46 

0.13 

7.2-10'2 

0.145 

2.44 

0.59 

0.05 

0.45 

0.98 

2.1 

53.0 

43.5 

1.0 

2.5 

1.84 

2.90 

.37 

.07 

3.6-102 

.145 

2.32 

.57 

.07 

.36 

.97 

3.7 

45.1 

36.1 

13.0 

5.8 

— 

— 

— 

— 

3.2-10'3 

.006 

1.23 

.19 

.28 

.54 

1.00 

6.7 

61.1 

20.8 

12.1 

6.0 

1.69 

— 

— 

.37 

2.1-1 0*6 

.008 

1.59 

.37 

.22 

.52 

.98 

WC2 

2.1 

34.7 

61.8 

.2 

3.3 

1.71 

2.64 

.35 

.15 

6.3-103 

.145 

2.56 

.61 

.07 

.35 

.97 

3.7 

22.0 

39.5 

22.7 

15.8 

1.64 

3.18 

.48 

.33 

2.2-10'2 

.106 

1.25 

.20 

.22 

.48 

.96 

6.7 

59.0 

14.5 

12.7 

13.8 

1.43 

2.95 

.52 

.47 

5.0-10-4 

.028 

1.15 

.13 

.20 

.52 

.95 

8.2 

27.4 

30.1 

30.7 

11.8 

1.52 

2.92 

.48 

.40 

9.8-10-6 

.014 

1.41 

.29 

.21 

.47 

.95 

9.8 

15.9 

76.1 

3.7 

4.3 

1.65 

2.48 

.34 

.18 

3.0-103 

.068 

1.51 

.34 

.10 

.32 

.97 

WC3 

0.9 

38.7 

45.3 

4.8 

11.2 

1.64 

2.48 

.34 

.11 

8.5-10'3 

.145 

2.03 

.51 

.13 

.34 

1.00 

1.5 

— 

— 

— 

— 

— 

— 

— 

— 

1.7-104 

— 

— 

— 

— 

— 

— 

8.2 

92.4 

5.0 

1.7 

.9 

1.75 

2.60 

.33 

.21 

2.3-103 

.043 

1.36 

.27 

.08 

.32 

.97 

13.7 

20.7 

62.4 

9.8 

7.1 

1.64 

2.71 

.39 

.11 

9.9-10'3 

.162 

1.53 

.35 

.10 

.39 

1.00 

19.8 

11.3 

58.9 

22.6 

7.2 

1.72 

2.88 

.40 

.19 

5.8-103 

.068 

1.37 

.27 

.13 

.39 

.98 

WC4a 

2.1 

78.7 

15.9 

1.9 

3.5 

1.81 

3.16 

.43 

.25 

6.4- Id2 

.145 

2.03 

.51 

.15 

.42 

.99 

11.3 

48.5 

32.7 

11.7 

7.1 

1.80 

2.69 

.33 

.20 

4.2-10-4 

.049 

1.48 

.32 

.22 

.33 

1.00 

WG1 

3.0 

68.5 

27.9 

1.4 

2.2 

1.68 

2.97 

.43 

.25 

2.0-10-2 

.146 

1.39 

.28 

.05 

.43 

1.00 

3.7 

72.8 

24.0 

1.6 

1.6 

1.84 

2.81 

.35 

.13 

4.6-102 

.145 

1.51 

.34 

.04 

.34 

.98 

13.4 

43.9 

52.9 

.7 

2.5 

1.37 

2.58 

.47 

.10 

4.0- Id2 

.145 

1.98 

.50 

.07 

.47 

.99 

22.6 

66.0 

27.8 

5.7 

.5 

1.59 

2.66 

.40 

.19 

2.0-102 

.145 

1.29 

.23 

.05 

.39 

.97 
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Table  6.  Physical  and  hydraulic  properties  of  core  material ,  Sierra  Vista  subwatershed,  Arizona — Continued 


van  Genuchten  fitting  parameters 


Depth 
to  top  of 
7.6-cm 
interval 
(m) 

Gravel 

(%) 

Sand 

(%) 

Silt 

(%) 

Clay 

(%) 

Bulk 

density 

(g/cm3) 

Particle 

density 

(g/cm3) 

Porosity 

(cm3/cm3) 

Volumetric 
water  content, 
field 

(cm3/cm3) 

Saturated 

hydraulic 

conductivity 

(cm/s) 

a 

(1/cm) 

n 

(dimen¬ 

sionless) 

m 

(1-1/n) 

Residual 

water 

content 

(cm3/cm3) 

Saturated 

water 

content 

(cm3/cm3) 

R2 

Basin-floor  boreholes 

BF1 

0.3 

41.7 

43.3 

7.0 

8.0 

1.52 

2.51 

0.39 

<0.02 

2.  MO2 

0.076 

2.11 

0.53 

0.09 

0.40 

1.00 

4.0 

68.0 

25.8 

2.5 

3.7 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

5.5 

76.5 

16.5 

1.7 

5.3 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

7.0 

64.1 

28.5 

3.1 

4.3 

1.77 

2.72 

.35 

.11 

1 .4*  10 2 

.137 

1.72 

.42 

.14 

.35 

.99 

8.5 

34.8 

34.5 

21.0 

9.7 

1.60 

2.92 

.45 

.26 

5.1-1 0-3 

.073 

1.28 

.22 

.18 

.46 

.95 

10.1 

77.0 

18.5 

2.4 

2.1 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

13.1 

60.7 

24.6 

7.7 

7.0 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

WC4c 

1.1 

41.1 

54.4 

1.0 

3.5 

1.70 

2.54 

.33 

.09 

1.5  102 

.166 

1.58 

.37 

.03 

.33 

.98 

2.4 

51.6 

43.9 

1.7 

2.8 

1.69 

2.57 

.34 

.07 

5.  MO'3 

.145 

1.81 

.45 

.05 

.34 

.99 

3.8 

23.9 

68.5 

2.1 

5.5 

1.66 

2.76 

.40 

.15 

1.5103 

.071 

1.89 

.47 

.06 

.40 

1.00 

9.9 

68.3 

19.1 

9.3 

3.3 

— 

— 

— 

— 

1.9  105 

.025 

3.89 

.74 

.27 

.40 

.99 

11.3 

56.6 

31.0 

9.6 

2.8 

1.83 

2.87 

.36 

.24 

7.2- 10'6 

.021 

2.09 

.52 

.27 

.36 

.97 

WG4 

4.3 

47.6 

26.0 

17.2 

9.2 

1.24 

2.38 

.48 

.19 

1.3  102 

.122 

1.21 

.17 

.15 

.48 

1.00 

13.4 

36.3 

37.2 

18.6 

7.9 

1.30 

2.61 

.50 

.10 

1.7  102 

.077 

1.63 

.39 

.17 

.47 

.96 

WG5 

.3 

34.0 

47.2 

15.6 

3.2 

1.41 

2.71 

.48 

.10 

1.5  102 

.126 

1.40 

.29 

.10 

.48 

.98 
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Table  7.  Isotope  and  chloride  data  for  composite  precipitation  samples,  July  2000— April  2003,  Sierra  Vista  subwatershed,  Arizona 


[cm,  centimeters;  mg/cm2,  milligrams  per  square  centimeter;  %c,  per  mil;  TU,  tritium  units;  — ,  no  data] 


Begin 

collection 

End 

collection 

Precipitation 

amount 

Chloride 
deposition 
rate  for  time 
period 

Delta 

deuterium 

Delta 

oxygen-18 

Tritium 

Tritium  error 

date 

date 

(cm) 

(mg/cm2) 

(%o) 

(%o) 

(TU) 

(TU) 

Upper  Babocomari  River  streamflow-gaging  station 


07-27-2000 

10-17-2000 

25.7 

0.021 

-34.0 

-5.72 

— 

— 

10-17-2000 

04-11-2001 

19.1 

.011 

-73.6 

-10.97 

3.3 

0.4 

04-11-2001 

10-11-2001 

26.4 

.013 

-38.7 

-5.31 

4.7 

.3 

10-11-2001 

04-25-2002 

5.3 

.004 

-90.7 

-12.80 

3.3 

.3 

04-25-2002 

10-17-2002 

24.1 

.010 

-54.7 

-8.17 

4.8 

.3 

10-17-2002 

04-23-2003 

___ 

___ 

-80.1 

-11.47 

4.2 

.3 

Greenbush  Draw  streamflow-gaging  station 


07-26-2000 

10-18-2000 

27.2 

.011 

-38.5 

-6.05 

___ 

___ 

10-18-2000 

04-11-2001 

15.5 

.017 

-65.6 

-10.41 

3.2 

.3 

04-11-2001 

10-11-2001 

14.5 

.010 

-46.5 

-6.12 

4.6 

.3 

10-11-2001 

04-25-2002 

4.3 

.001 

-133.4 

-17.71 

2.6 

.3 

04-25-2002 

10-17-2002 

15.0 

.012 

-34.2 

-5.23 

4.8 

.3 

10-17-2002 

04-21-2003 

___ 

___ 

-68.8 

-10.37 

1.7 

.3 

Streamflow-Data  Collection  and  Analysis 

Knowledge  of  ephemeral  streamflow  timing  and  duration 
in  the  tributaries  was  needed  to  determine  the  availability  of 
water  for  ephemeral-stream  channel  infiltration.  In  this  study, 
streamflow  availability  was  estimated  using  frequency  and 
duration  of  flow  determined  by  USGS  streamflow-gaging 
stations,  ARS  flumes,  and  temperature  methods.  Onset 
Stow  Away  TidbiT  temperature  sensors  were  buried  about 
15  cm  below  the  surface  of  streambed  sediments  at  7  of  the 
12  tributary  borehole  sites  (table  2)  and  at  selected  channel 
locations  between  boreholes  (fig.  1).  Sediment  temperatures 
were  logged  every  15  to  30  minutes.  A  rapid  drop  in 
temperature  of  more  than  0.2°C  indicated  onset  of  streamflow; 
a  low-temperature  inflection  point  in  the  thermograph 
was  interpreted  as  the  cessation  of  streamflow  (Constantz 
and  others,  2001;  Gungle,  2003).  Durations  of  individual 
ephemeral-streamflow  events  ranged  from  15  minutes  to 
48  days  (table  8). 

Data-Analysis  Methods 

Infiltration  is  defined  as  the  flow  of  water  from  above 
ground  into  the  unsaturated  zone;  recharge  is  defined  as  the 
flow  of  water  across  the  regional  water  table.  The  rate  of 
infiltration,  or  infiltration  flux,  is  strongly  dependent  on  the 
spatial  and  temporal  distribution  of  source  water  and  the 
hydraulic  and  physical  properties  of  the  subsurface  sediments. 
At  the  onset  of  infiltration,  the  initial  infiltration  flux  is  at  its 
highest  and  is  extremely  transient.  As  infiltration  continues 
over  time,  the  infiltration  flux  decreases  to  a  nearly  constant 
rate  approaching  the  effective  mean  saturated  hydraulic 
conductivity  of  the  subsurface  sediments. 


Infiltration  can  lead  directly  to  recharge  when  the  wetting 
front  reaches  the  regional  water  table  while  infiltration  is  still 
occurring.  In  the  arid  and  semiarid  Southwest,  however,  deep 
water  tables  and  short  infiltration  times  lead  to  infiltration 
ending  long  before  the  wetting  front  crosses  the  water 
table.  In  these  instances,  water  that  has  infiltrated  into  the 
unsaturated  zone  can  be  removed  by  evapotranspiration  or, 
if  the  water  has  moved  past  the  root  zone,  remain  in  storage 
until  it  is  displaced  by  water  from  a  subsequent  infiltration 
event.  Unsaturated-zone  water  in  storage  below  the  root 
zone  can  continue  to  redistribute  over  time.  This  long-term 
movement  of  water,  or  percolation  flux,  occurs  both  vertically, 
by  drainage,  and  horizontally  in  response  to  heterogeneity  or 
anisotropy  of  the  sediments.  Long-term  percolation  can  result 
in  recharge. 

Unsaturated  Hydraulic  Properties 

Moisture-release  curves  describe  the  nonlinear  decrease 
in  soil- water  content  at  decreasing  matric  potentials.  The  van 
Genuchten  (1980)  algebraic  expression  relates  soil- water 
content  to  matric  potential  by: 

ev=(l  +  (a\hp\nm  ,  (3) 

where 

6v  =  volumetric  soil- water  content  [L3/L3], 
a  =  van  Genuchten  fitting  parameter  [1/L], 
h  =  matric  potential  [L], 

n  =  van  Genuchten  fitting  parameter  (dimensionless), 
and 

m  =  van  Genuchten  fitting  parameter  equal  to  l-l/n 
(dimensionless). 
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Table  8.  Streamflow  durations  at  streambed-temperature  sensors,  Agricultural  Research  Service  flumes,  and  U.S.  Geological 
Survey  streamflow-gaging  stations;  and  distance  of  duration-measurement  locations  from  the  mountain  front,  Sierra  Vista 
subwatershed,  Arizona 

[hr,  hours;  km,  kilometers;  — ,  no  data;  values  in  parentheses  are  distance  from  the  mountain  front  in  kilometers] 


Streamflow  duration 
(hr) 

Banning  Creek 

Streamflow 
begin  date 

Gaging  station  094707001  Temperature  sensor  3778221 

(0.0  km)  (0.0  km) 

Temperature  sensor  at  BC1 
(3.9  km) 

07-15-2001 

1,152 

— 

0 

07-20-2001 

0 

— 

1.5 

07-22-2001 

0 

— 

11 

07-25-2001 

0 

— 

1.0 

08-01-2001 

0 

— 

1.5 

10-03-2001 

0 

___ 

2.3 

08-03-2002 

— 

909 

13 

08-04-2002 

— 

0 

1.0 

08-05-2002 

— 

0 

5.8 

Total 

1,152 

909 

37.1 

Greenbush  Draw 

Streamflow 
begin  date 

Temperature  sensor 
377815 
(1.8  km) 

Temperature  sensor 
at  GD2 
(5.6  km) 

Temperature  sensor 
374866 
(11.3  km) 

Temperature  sensor  Temperature  sensor 
at  GDI2  377820 

(12.4  km)  (18.5  km) 

01-09-2001 

1.0 

0 

0 

___ 

0 

04-06-2001 

7.5 

0 

0 

— 

0 

07-09-2001 

0 

2.0 

0 

___ 

0 

07-10-2001 

0 

3.0 

0 

___ 

0 

07-13-2001 

0 

.5 

12 

— 

0 

07-17-2001 

0 

0 

15 

— 

0 

07-18-2001 

0 

0 

.8 

___ 

2.5 

07-21-2001 

0 

13 

0 

— 

5.3 

07-24-2001 

0 

2.3 

0 

— 

0 

07-25-2001 

0 

6.0 

0 

— 

0 

07-28-2001 

0 

0 

1.8 

— 

16 

08-05-2001 

0 

16 

1.0 

— 

0 

08-11-2001 

0 

3.3 

0 

— 

0 

08-12-2001 

0 

5.3 

0 

___ 

0 

08-13-2001 

0 

7.0 

.8 

— 

0 

09-12-2001 

0 

0 

0 

— 

16 

07-14-2002 

0 

0 

0 

0 

1.5 

07-21-2002 

17 

0 

0 

0 

0 

08-01-2002 

0 

2.0 

.5 

0 

0 

08-02-2002 

0 

0 

0 

1.0 

0 

08-03-2002 

0 

.8 

0 

0 

0 

08-18-2002 

0 

0 

0 

0 

19 

08-29-2002 

0 

0 

0 

0 

0 

09-09-2002 

7.0 

0 

0 

0 

0 

09-10-2002 

0 

1.3 

0 

0 

0 

Total 

37.5 

62.5 

94.4 

1 

60.3 

See  footnotes  at  end  of  table. 
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Table  8.  Streamflow  durations  at  streambed-temperature  sensors,  Agricultural  Research  Service  flumes,  and  U.S.  Geological 
Survey  streamflow-gaging  stations;  and  distance  of  duration-measurement  locations  from  the  mountain  front,  Sierra  Vista 
subwatershed,  Arizona — Continued 


Streamflow  duration 
(hr) 

Miller  Canyon  Wash 

Temperature  sensor 
Streamflow  377806 

be9,ndate  (0.0  km) 

Temperature  sensor 
at  MC22 
(0.0  km) 

Temperature  sensor 
377797 
(3.9  km) 

Temperature  sensor 
at  MCI2 
(7.1  km) 

Temperature  sensor 
374919 
(10.9  km) 

04-06-2001 

0 

— 

6.0 

— 

0 

06-25-2001 

0 

— 

2.0 

— 

0 

07-18-2001 

0 

— 

2.3 

— 

0 

07-28-2001 

.5 

— 

3.5 

— 

0 

08-01-2001 

0 

— 

1.8 

— 

0 

08-05-2001 

15 

— 

3.0 

— 

1.5 

08-16-2001 

.8 

— 

0 

— 

3.5 

08-17-2001 

0 

— 

0 

— 

2.5 

10-03-2001 

0 

— 

0 

— 

2.5 

07-16-2002 

20 

18 

0 

0 

0 

07-25-2002 

0 

0 

0 

0 

2.3 

08-19-2002 

0 

0 

1.0 

0 

0 

09-03-2002 

0 

0 

2.0 

0 

0 

09-09-2002 

0 

0 

1.3 

0 

4.5 

09-10-2002 

0 

0 

3.5 

0 

0 

Total 

36.3 

18 

26.4 

0 

16.8 

Walnut  Gulch3 

Streamflow 

Flume  11 

Flume  6 

Flume  1 

begin  date 

(7.2  km) 

(13.2  km) 

(25.6  km) 

01-06-2001 

0 

0 

0.7 

01-13-2001 

38 

0 

0 

01-16-2001 

82 

0 

0 

07-07-2001 

1.2 

0 

0 

07-13-2001 

.5 

0 

0 

08-11-2001 

1.7 

0 

0 

08-13-2001 

.8 

0 

2.6 

08-14-2001 

0 

0 

2.9 

08-17-2001 

0 

0 

.4 

09-14-2001 

0 

0 

3.4 

07-08-2002 

1.2 

0 

0 

07-16-2002 

0 

.5 

0 

07-19-2002 

1.2 

2.8 

0 

07-26-2002 

0 

1.3 

10 

08-01-2002 

.6 

5.0 

0 

08-04-2002 

3.0 

6.0 

5.0 

08-06-2002 

1.0 

0 

0 

09-08-2002 

0 

0 

.2 

09-09-2002 

0 

3.5 

0 

09-11-2002 

1.3 

1.0 

0 

Total 

132.5 

20.1 

25.2 

See  footnotes  at  end  of  table. 
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Table  8.  Streamflow  durations  at  streambed-temperature  sensors,  Agricultural  Research  Service  flumes,  and  U.S.  Geological  Survey 
streamflow-gaging  stations;  and  distance  of  duration-measurement  locations  from  the  mountain  front,  Sierra  Vista  subwatershed, 
Arizona — Continued 


Streamflow  duration 
(hr) 

Woodcutters  Wash 

0jt  Temperature  sensor 

Streamflow  atWC4a2 

begMlda,e  (0.4  km) 

Temperature  sensor 
375010 
(3.1  km) 

Temperature  sensor 
at  WC32 
(7.2  km) 

Temperature  sensor 
377793 
(8.0  km) 

Temperature  sensor 
at  WC1 
(16.4  km) 

04-05-2001 

0 

0 

— 

13 

0 

04-06-2001 

0 

.5 

— 

0 

3.0 

06-19-2001 

0 

16 

— 

12.5 

0 

06-25-2001 

0 

0 

— 

2.0 

.3 

07-05-2001 

0 

13 

— 

0 

0 

07-07-2001 

0 

22 

— 

1.5 

.5 

07-09-2001 

0 

5.0 

— 

1.3 

0 

07-16-2001 

0 

0 

— 

0 

.8 

07-24-2001 

0 

5.0 

— 

0 

0 

07-25-2001 

0 

7.5 

— 

0 

0 

07-28-2001 

0 

13 

— 

2.8 

.3 

08-01-2001 

0 

1.0 

— 

0 

0 

08-05-2001 

0 

2.0 

— 

0 

0 

08-11-2001 

0 

7.5 

— 

0 

0 

08-13-2001 

.5 

6.5 

— 

0 

14 

08-14-2001 

0 

1.5 

— 

0 

0 

08-16-2001 

0 

2.0 

— 

0 

0 

08-17-2001 

0 

8.5 

— 

0 

4.5 

08-29-2001 

0 

5.5 

— 

0 

12 

09-12-2001 

0 

6.5 

— 

0 

0 

07-08-2002 

0 

1.0 

2.5 

.8 

12 

07-16-2002 

0 

.5 

1.5 

1.0 

.5 

07-18-2002 

0 

0 

0 

0 

.8 

07-26-2002 

0 

0 

0 

0 

2.0 

08-03-2002 

0 

0 

2.0 

0 

0 

08-08-2002 

0 

0 

0 

0 

1.8 

08-16-2002 

0 

0 

0 

0 

.8 

08-19-2002 

0 

0 

1.0 

1.0 

.8 

09-08-2002 

0 

0 

0 

0 

15 

09-09-2002 

0 

0 

0 

0 

3.3 

10-07-2002 

0 

0 

1.5 

2.5 

1.0 

Total 

0.5 

124.5 

8.5 

38.4 

73.4 

Streamflow-gaging  station  09470700  and  temperature  sensor  377822  are  at  the  same  location. 

2Data  not  used  in  analysis  of  streamflow  duration  versus  distance  from  the  mountain  front. 

Provisional  data  provided  by  the  U.S.  Department  of  Agriculture,  Agricultural  Research  Service,  Southwest  Watershed  Research  Center. 
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For  this  study,  the  pressures  (in  cm  water)  at  which  soil- water 
contents  were  determined  are  considered  to  be  equivalent  to 
matric  potentials. 

The  RETC  (RETention  Curve)  computer  code  (van 
Genuchten  and  others,  1991)  was  used  to  fit  moisture-release 
curves  to  the  induced-pressure  and  measured  water-content 
data  to  determine  the  van  Genuchten  fitting  parameters  of  a 
and  n  (table  6).  All  curves  were  fit  so  that  their  coefficients  of 
determination  (. R 2)  for  the  regression  of  the  measured  values 
versus  the  fitted  values  were  greater  than  0.95.  Volumetric 
residual  and  saturated  soil- water  contents  (0R  and  0S,  in  L3/L3) 
were  also  determined  using  the  RETC  code. 

Unsaturated-Zone  Chloride  Mass  Balance 

The  mass  of  chloride  (Cl)  in  the  unsaturated  zone  can  be 
used  to  estimate  long-term  vertical  water  movement.  Chloride 
is  continuously  deposited  on  the  land  surface  by  precipitation 
and  dry  fallout.  The  high  solubility  of  Cl  enables  its  transport 
into  the  subsurface  by  infiltrating  water.  Because  Cl  is 
essentially  nonvolatile  and  Cl  uptake  by  plants  is  minimal, 

Cl  is  retained  in  the  sediment  when  water  is  removed  by 
evaporation  and  transpiration.  An  increase  in  Cl  within  the 
root  zone  of  the  shallow  subsurface,  therefore,  is  proportional 
to  the  amount  of  water  lost  by  evapotranspiration  (Allison  and 
Hughes,  1978).  In  areas  where  active  infiltration  is  occurring, 
an  increase  in  Cl  in  the  shallow  subsurface  will  generally 
be  absent,  and  concentrations  will  be  very  low  through  the 
unsaturated  zone.  In  areas  where  little  to  no  active  infiltration 
is  occurring,  an  increase  in  Cl  in  the  shallow  subsurface  will 
be  present,  and  concentrations  will  be  relatively  high  through 
the  unsaturated  zone. 

If  the  Cl  deposition  rate  on  the  land  surface  is  known,  the 
average  travel  time  of  Cl  (ta\  in  t)  to  a  depth  in  the  unsaturated 
zone  (z,  in  L)  can  be  calculated  by: 


where 

ClsoiI  =  chloride  mass  in  the  sample  interval  [M/L3],  and 
Cldep  =  chloride  deposition  rate  [M/L2/t] . 

The  above  calculation  entails  several  assumptions:  (1)  flow 
in  the  unsaturated  zone  is  downward- vertical  and  piston  type, 
(2)  bulk  precipitation  (precipitation  plus  dry  fallout)  is  the  only 
source  of  Cl  and  there  are  no  mineral  sources  of  Cl,  (3)  the  Cl 
deposition  rate  has  stayed  constant  over  time,  and  (4)  there  is 
no  recycling  of  Cl  within  the  unsaturated  zone. 

This  model  of  Cl  movement  predicts  that  in  interdrainage 
areas  with  little  to  no  infiltration,  Cl  concentrations  should 
increase  through  the  root  zone  and  then  remain  constant  down 
to  the  water  table.  Numerous  published  chloride  profiles, 
however,  show  that  Cl  concentrations  commonly  decrease 
below  the  peak  concentration  in  the  root  zone.  Some  of 


the  above  assumptions,  therefore,  may  not  be  valid  for  all 
unsaturated-zone  systems.  The  decrease  in  Cl  concentrations 
below  the  peak  has  been  attributed  to  nonpiston  flow  (Sharma 
and  Hughes,  1985),  ground- water  dilution  (Allison  and  others, 
1985),  or  nonsteady  flow  as  a  result  of  paleoclimatic  variations 
(Phillips  and  Stone,  1985;  Scanlon,  1991). 

For  some  geomorphic  settings,  such  as  ephemeral 
streams,  assumption  2  is  not  valid  because  streamflow 
provides  an  additional  source  of  Cl  to  the  subsurface.  To 
correctly  determine  travel  time  below  ephemeral  streams  with 
the  Cl  mass  balance  method,  the  mass  of  Cl  in  streamflow 
needs  to  be  included  in  equation  4.  In  this  study,  ephemeral 
streamflow  discharges  and  Cl  concentrations  were  not 
quantified;  therefore,  the  Cl  mass  method  was  not  used  to 
determine  travel  times  below  ephemeral  streams. 

Chloride  deposition  rates  in  the  Sierra  Vista  subwatershed 
were  determined  from  the  composite  precipitation  samples 
collected  for  this  study  (table  7).  Annual  rates  ranged  from 
0.011  to  0.032  mg/cm2  yr.  The  average  annual  Cl  deposition 
rate  for  the  subwatershed  was  determined  to  be 
0.022  mg/cm2  yr. 

Meteoric  Tritium 

Tritium  (3H),  a  naturally  occurring  radioactive  isotope  of 
hydrogen,  can  be  used  to  estimate  the  age  of  subsurface  water 
and  infiltration  fluxes  into  the  subsurface.  Large  quantities  of 
tritium  were  released  to  the  atmosphere  during  thermonuclear- 
weapons  testing  from  1952  until  the  late  1960s;  maximum 
releases  occurred  in  the  early  1960s.  As  a  result,  the  amount 
of  tritium  in  precipitation  sharply  increased  during  testing  as 
tritium  was  introduced  into  the  water  cycle,  and  decreased 
after  testing  ended.  The  amount  of  tritium  in  subsurface  water 
at  a  given  time  is  a  function  of  the  amount  of  tritium  in  the 
atmosphere  when  infiltration  occurred  and  the  radioactive- 
decay  rate  of  tritium.  If  flow  in  the  unsaturated  zone  is 
assumed  to  be  downward- vertical  and  piston  type,  the  average 
infiltration  flux  (q.,  in  L/T)  can  be  estimated  by: 

_  A2 

q'  At 9"  ’  (5) 

where 

A z  =  depth  of  maximum  tritium  activity  [L], 

At  =  elapsed  time  between  sampling  and  maximum 
atmospheric  tritium  activity  [t]  (38  years  for  this 
study),  and 

0v  =  volumetric  soil- water  content  [L3/L3] . 

Tritium  activity  in  the  composite  precipitation  samples 
collected  during  this  study  ranged  from  1.7  to  4.2  TU  for 
winter  precipitation  and  from  4.6  to  4.8  TU  for  summer 
precipitation  (table  7).  Pore  water  having  a  tritium  activity 
within  this  range  was,  therefore,  considered  to  be  recent 
infiltration.  Pore  water  having  an  activity  of  less  than  0.5  TU 
was  considered  to  have  infiltrated  prior  to  1952. 
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Stable  Isotopes 


The  stable-isotope  ratios  of  oxygen  (180/160)  and 
hydrogen  (^H^H)  can  be  used  to  estimate  seasonal  cycles 
of  infiltration  and  to  determine  if  pore  water  has  undergone 
evaporation.  Ratios  determined  in  this  study  are  reported 
relative  to  Vienna  Standard  Mean  Ocean  Water  (VSMOW). 
Differences  from  the  standard  are  expressed  by  delta  notation 
(S)  in  per  mil  (%o).  The  delta  notation  is  computed  by: 


6  = 


Rx  ~Ri 

RStd 


•1,000, 


(6) 


where 

Rx  =  ratio  of  isotopes  measured  in  the  sample,  and 
RStd  =  ratio  of  same  isotopes  in  standard. 


The  relation  between  <5180  and  S2 H  in  precipitation  can  be 
represented  by  the  global  meteoric  water  line: 


S2 H  =  8<5180  +  10  (Craig,  1961).  (7) 


Mass  differences  between  180  and  160  and  between  2H  and  lH 
enable  isotopic  fractionation  during  phase  changes,  such  as 
evaporation.  Water  molecules  consisting  of  160  and  lH  evapo¬ 
rate  preferentially,  resulting  in  water  enriched  in  180  and  2H. 
Lines  describing  the  relation  between  <5180  and  <52H  in  enriched 
water  have  smaller  slopes  than  the  slope  of  the  global  meteoric 
water  line. 

The  amount  of  each  isotope  in  precipitation  is  dependant 
on  the  temperature  and  source  location  of  the  precipitation. 
Because  precipitation  temperatures  and  sources  vary  by 
season  in  the  Sierra  Vista  subwatershed,  stable-isotope  values 
in  precipitation  show  seasonal  variation.  The  <5180  for  the 
composite  precipitation  samples  collected  for  this  study  ranged 
from  -5.23  to  -8.17  %c  for  summer  precipitation  and  from 
-10.37  to  -17.71  %o  for  winter  precipitation  (table  7).  The  S2U 
ranged  from  -34.0  to  -54.7  %c  for  summer  precipitation  and 
from  -65.6  to  -133.4  %c  for  winter  precipitation  (table  7). 


Heat  as  a  Tracer 

Sediment- temperature  profiles  can  be  used  to  estimate 
one-dimensional  infiltration  fluxes  in  the  unsaturated  zone 
below  the  beds  of  ephemeral-stream  channels  (Constantz  and 
Thomas,  1996;  Constantz  and  others,  2003).  Temperatures  of 
streambed  surfaces  are  predominantly  influenced  by  heating 
from  solar  radiation  and  ambient  air  temperatures.  Beneath  dry 
streambeds,  shallow  sediment  temperatures  vary  diurnally  and 
seasonally  in  response  to  downward  conductive  heat  transport. 
Deeper  sediment  temperatures  are  less  affected  by  diurnal  or 
seasonal  air  temperatures,  and,  therefore,  do  not  greatly  vary 
over  time.  During  periods  of  stream-channel  infiltration,  heat 
transport  into  the  subsurface  by  advection  predominates  and 
transport  by  conduction  is  considered  negligible. 


An  analytical  method  (Taniguchi  and  Sharma, 

1990;  Constantz  and  Thomas,  1996)  was  used  as  a  first 
approximation  to  calculate  infiltration  flux  (q  in  L/t)  from 
measured  sediment  temperatures  at  each  ephemeral- stream 
site  (table  9): 


Qi  = 


PjM 

A^ 


(8) 


where 

A z  =  vertical  distance  between  temperature 
sensors  [L], 

At  =  travel  time  of  the  temperature  perturbation 
between  sensors  [t], 

PbSoii  ~  density  °f  sediment  [M/L3], 

Cs  =  heat  capacity  of  the  dry  sediment  [E/L3/T], 
phwater  =  bulk  density  of  water  [M/L3],  and 
Cw  =  heat  capacity  of  water  [E/L3/T]. 


Onset  of  streamflow  was  determined  from  streambed-tem- 
perature  data.  The  distance  between  temperature  sensors  was 
known  for  each  borehole,  and  the  travel  time  of  the  tempera¬ 
ture  signal  was  measured.  The  pb  and  Cs  of  the  sediment  were 
determined  from  cores.  For  depth  intervals  where  more  than 
one  property  value  was  available,  average  values  were  used. 
The  pbwater  is  a  constant  1.0  g/ cm3,  and  the  Cw  is  a  constant 
4.2  10rj/m3oC. 

The  infiltration  flux  determined  by  the  use  of  the 
analytical  method  represents  the  initial,  transient  infiltration 
flux.  This  method  is  only  valid  for  cases  where  pore- water 
velocities  are  high  enough  that  heat  transport  by  conduction 
is  negligible  compared  with  heat  transport  by  advection. 

The  high  K  values  determined  for  stream- alluvium  cores 
indicate  that  this  situation  is  typical  for  most  of  the  subsurface 
sediments  beneath  ephemeral- stream  channels  in  the  Sierra 
Vista  subwatershed.  The  analytical  method  is  also  only  valid 
for  infiltration  under  nonisothermal,  continuous- streamflow 
loss;  this  condition  is  assumed  in  the  ephemeral- stream 
channels. 

A  more  comprehensive,  rigorous  numerical  approach 
was  used  to  account  for  heat  transport  with  water  flow  in  the 
unsaturated  zone.  Convective  and  advective  heat  transport  can 
be  expressed  by: 


-[evcu+(i-$)c.] 


d_ 

dz 

d_ 

dz 


KM 


dT_ 

dz 
dT 

SCwDh  — 
dz 


(Kipp,  1987),  (9) 


where 

0v  =  volumetric  soil- water  content  [L3/L3], 

Kt  =  thermal  conductivity  [E/t/L/T], 

T  =  temperature  [T],  and 
Dh  =  a  thermomechanical  dispersion  tensor  [L2/t] . 
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Table  9.  Streamflow  characteristics  and  input  parameters  used  for  the  infiltration-flux  analytical  method  using  heat  as  a  tracer. 
Sierra  Vista  subwatershed,  Arizona 

[m,  meters;  hr,  hour;  g/cm3,  grams  per  cubic  centimeter;  J/m3-°C,  joules  per  cubic  meter  degrees  Celsius;  m/hr,  meters  per  hour; 

Flux  type:  Infiltration,  infiltration  flux;  Percolation,  percolation  flux;  NC,  not  calculated] 


Upper 

Lower 

Travel  time  to 

Heat  capacity. 

Calculated 

Date  of 

measuring  point 

measuring  point 

measuring  point 

Bulk  density 

dry  soil 

transient  flux 

flow  event 

(m) 

(m) 

(hr) 

(g/cm3) 

(J/m3oC-106) 

(m/hr) 

Flux  type 

BC1 

07-20-2001 

0.0 

1.5 

1.0 

1.84 

2.08 

1.4 

Infiltration 

1.5 

3.0 

0 

do. 

do. 

NC 

NC 

3.0 

4.6 

.5 

do. 

do. 

2.8 

Percolation 

4.6 

6.1 

2.0 

do. 

do. 

.7 

Percolation 

07-22-2001 

0 

1.5 

.5 

do. 

do. 

2.8 

Infiltration 

1.5 

3.0 

.5 

do. 

do. 

2.8 

Infiltration 

3.0 

4.6 

1.5 

do. 

do. 

.9 

Percolation 

4.6 

6.1 

1.0 

do. 

do. 

1.4 

Percolation 

07-25-2001 

0 

1.5 

1.0 

do. 

do. 

1.4 

Infiltration 

1.5 

3.0 

0 

do. 

do. 

NC 

NC 

3.0 

4.6 

1.0 

do. 

do. 

1.4 

Percolation 

4.6 

6.1 

1.5 

do. 

do. 

.9 

Percolation 

08-01-2001 

0 

1.5 

1.0 

do. 

do. 

1.4 

Infiltration 

1.5 

3.0 

0 

do. 

do. 

NC 

NC 

3.0 

4.6 

1.0 

do. 

do. 

1.4 

Percolation 

4.6 

6.1 

1.0 

do. 

do. 

1.4 

Percolation 

10-03-2001 

0 

1.5 

1.0 

do. 

do. 

1.4 

Percolation 

1.5 

3.0 

0 

do. 

do. 

NC 

NC 

3.0 

4.6 

.5 

do. 

do. 

2.8 

Percolation 

4.6 

6.1 

1.5 

do. 

do. 

.9 

Percolation 

08-03-2002 

0 

1.5 

1.0 

do. 

do. 

1.4 

Percolation 

1.5 

3.0 

0 

do. 

do. 

NC 

NC 

3.0 

4.6 

1.0 

do. 

do. 

1.4 

Percolation 

4.6 

6.1 

1.5 

do. 

do. 

.9 

Percolation 

08-04-2002 

0 

1.5 

1.0 

do. 

do. 

1.4 

Percolation 

1.5 

3.0 

0 

do. 

do. 

NC 

NC 

3.0 

4.6 

1.0 

do. 

do. 

1.4 

Percolation 

4.6 

6.1 

.5 

do. 

do. 

2.8 

Percolation 

GDI 

08-02-2002 

0 

1.5 

1.5 

1.86 

1.74 

.8 

Percolation 

1.5 

3.0 

0 

do. 

do. 

NC 

NC 

3.0 

6.1 

3.5 

do. 

do. 

.7 

Percolation 

6.1 

9.1 

17.5 

do. 

do. 

.1 

Percolation 

GD2 

08-01-2002 

0 

1.5 

1.5 

1.46 

2.54 

.9 

Percolation 

1.5 

3.0 

1.5 

do. 

do. 

.9 

Percolation 

3.0 

6.1 

5.0 

do. 

do. 

.5 

Percolation 

09-10-2002 

0 

1.5 

2.0 

do. 

do. 

.7 

Percolation 

1.5 

3.0 

2.5 

do. 

do. 

.5 

Percolation 

3.0 

6.1 

1.5 

do. 

do. 

1.8 

Percolation 

6.1 

9.1 

1.5 

do. 

do. 

1.8 

Percolation 

9.1 

12.2 

10 

do. 

do. 

.3 

Percolation 

WC1 

07-07-2001 

0 

1.5 

.5 

1.67 

1.24 

1.5 

Infiltration 

1.5 

3.0 

.5 

do. 

do. 

1.5 

Infiltration 

07-16-2001 

0 

1.5 

.5 

do. 

do. 

1.5 

Percolation 

1.5 

3.0 

.5 

do.. 

do. 

1.5 

Percolation 

07-28-2001 

0 

1.5 

.5 

do. 

do. 

1.5 

Percolation 

1.5 

3.0 

.5 

do. 

do. 

1.5 

Percolation 

07-16-2002 

0 

1.5 

2.0 

do. 

do. 

.4 

Percolation 

1.5 

3.0 

2.0 

do. 

do. 

.4 

Percolation 

WG1 

07-26-2002 

0 

1.5 

1.0 

1.76 

1.63 

1.0 

Infiltration 

1.5 

3.0 

.5 

do. 

do. 

2.1 

Infiltration 

3.0 

4.6 

.5 

do. 

do. 

2.1 

Infiltration 

4.6 

6.1 

5.5 

do. 

do. 

.2 

Percolation 

08-04-2002 

0 

1.5 

1.0 

do. 

do. 

1.0 

Infiltration 

1.5 

3.0 

0 

do. 

do. 

NC 

NC 

3.0 

4.6 

0 

do.. 

do. 

NC 

NC 

4.6 

6.1 

7.4 

do. 

do. 

.1 

Percolation 

09-08-2002 

0 

1.5 

1.0 

do. 

do. 

1.0 

Percolation 

WG2 

07-19-2002 

0 

3.0 

1.5 

‘1.76 

1.23 

1.0 

Infiltration 

Estimated  on  the  basis  of  the  value  for  borehole  WG1. 
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To  describe  water  flow  in  the  unsaturated  zone,  variably 
saturated  ground- water  flow  can  be  expressed  by: 

dO(zjt) 

dt 
__  d 

dz 

where 

K  =  unsaturated  hydraulic  conductivity  [L/t], 
hp  =  matric  potential  [L],  and 
h  =  total  head  [L]. 


KM 


dh(z ,  t) 


dz 


(Richards,  1931),  (10) 


The  numerical  code  VS2DH  (Healy  and  Ronan,  1996)  and  its 
graphical  user  interface  (Hsieh  and  others,  2000)  were  used  to 
solve  equations  9  and  10  simultaneously. 

One-dimensional  numerical  models  were  developed 
to  represent  conceptual  fluid  flow  and  heat  transfer  at  each 
ephemeral-stream  site  (tables  10  and  11).  Models  were  created 
only  for  sites  and  times  where  temperature  perturbations  were 
detected  at  depth  during  streamflow  events.  Models  were  not 
created  for  sites  and  events  where  temperature  perturbations 
were  not  detected  until  after  streamflow  ended.  Similarly, 
simulations  were  terminated  at  the  end  of  streamflow  because 
at  that  point,  heat  transfer  is  dominated  by  both  water  and 
vapor  flux,  a  condition  VS2DH  cannot  simulate. 


Table  10.  Sediment  properties  used  for  the  VS2DH  numerical  simulations.  Sierra  Vista  subwatershed,  Arizona 

[m,  meters;  cm/s,  centimeters  per  second;  cm3/cm3,  cubic  centimeters  per  cubic  centimeters;  1/cm,  inverse  centimeters;  W/m-°C,  watts  per  meter  degrees 
Celsius;  J/m3-°C,  joules  per  cubic  meter  degrees  Celsius] 


Borehole 

Depth  Depth  to 
to  top  of  bottom  of 
domain  domain 
(m)  (m) 

Initial  Adjusted 

saturated  saturated 

hydraulic  hydraulic 

conductivity  conductivity 
(cm/s)  (cm/s) 

Porosity 

(cm3/cm3) 

Residual 

water 

content 

(cm3/cm3) 

Thermal 
conductivity, 
Thermal  residual 

n  conductivity,  water 

a  (dimen-  saturation1  content 

(1/cm)  sionless)  (W/m°C)  (W/m°C) 

Heat 
capacity, 
dry  soil 
(J/m3oC-106) 

BC1 

0.0  3.0 

9.6- 10'3  1.5  10 2 

0.34 

0.10 

0.145  1.9 

2.5  2.1 

2.08 

WC1 

0  3.0 

5.4- lO2  4.0- 10'2 

.42 

.03 

145  2.4 

2.0  .8 

1.24 

WG1 

0  4.6 

3.3  102  3.5- lO2 

.39 

.05 

145  1.4 

2.5  .9 

1.63 

WG2 

0  3.0 

23.310'2  3.5- lO'2 

2.39 

2.05 

2 

145  21.4 

2.5  .8 

1.23 

Estimated  from  thermal  conductivity  at  field  conditions. 

^Estimated  on  the  basis  of  values  for  borehole  WG1. 

Table  11. 

Initial  and  boundary  conditions  for  the  VS2DH  numerical  simulations,  Sierra  Vista  subwatershed,  Arizona 

[m,  meters; c 

C,  degrees  Celsius;  cm3/cm3,  cubic  centimeter  per  cubic  centimeter;  m/hr,  meters  per  hour;  — ,  no  data] 

Initial  sediment  temperature  condition 

Initial  hydraulic 

(°C) 

condition 

Transport  boundary  condition 

Water 

Temperature 

Flux  into 

Simulation 

content1 

of  inflow 

domain2 

Borehole 

date 

0.0  m  1.5  m  3.0  m 

4.6  m 

i  6.1  m 

(cm3/cm3) 

(°C) 

(m/hr) 

BC1 

07-20-2001 

29.5  23.7  19.6 

___ 

___ 

0.15 

24.0 

0.5 

07-25-2001 

37.3  24.1  20.9 

___ 

— 

.25 

22.0 

.4 

08-01-2001 

25.0  24.6  21.9 

— 

— 

.25 

20.0 

.4 

WC1 

07-07-2001 

26.2  25.0  21.0 

— 

— 

.07 

24.0 

1.3 

WG1 

07-26-2002 

28.7  25.9  23.4 

21.5 

— 

.10 

25.0 

.5 

08-04-2002 

26.8  25.7  24.2 

22.6 

— 

.20 

25.0 

.4 

WG2 

07-19-2002 

28.2  —  18.0 

— 

— 

.15 

21.0 

.4 

Estimated  from  soil- water  content  at  time  of  sampling, 
determined  through  calibration. 
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Sediment  properties  were  determined  from  core 
measurements.  For  depth  intervals  where  more  than  one 
property  value  was  measured,  average  values  were  used  in 
the  calculation.  Thermal  conductivities  at  saturation  were 
estimated  on  the  basis  of  measured  thermal  conductivities  at 
ambient  conditions  (table  5).  All  layers  were  assumed  to  be 
isotropic  and  to  have  a  specific  storage  of  zero.  Streamflow 
water  temperatures  were  determined  from  values  measured 
by  the  surface  temperature  sensors.  Upper  boundaries  were 
defined  as  no  flow  for  the  hour  prior  to  the  onset  of  flow,  and 
then  as  specified  normal-flux  boundaries  during  streamflow. 
Initial  sediment  temperatures  were  determined  from  values 
measured  by  temperature  sensors  at  depth.  Initial  hydraulic 
conditions  were  estimated  as  soil- water  contents.  Lower 
boundaries  were  defined  as  possible  seepage  faces.  Onset 
and  length  of  streamflow  were  determined  from  streambed 
sediment-temperature  data. 

The  model  simulations  yielded  subsurface  sediment 
temperatures  due  to  advective  heat  transport.  In  all  cases, 
the  saturated  hydraulic  conductivity  had  to  be  adjusted 
to  accurately  simulate  the  temperature  inflection  point 
identifying  streamflow  arrival  time  (table  10).  The  models 
were  calibrated  by  adjusting  the  specified  infiltration  flux  at 
the  upper  boundary  to  minimize  the  difference  between  the 
calculated  and  the  measured  sediment  temperatures  at  depth. 

The  infiltration  flux  calibrated  for  the  upper  boundary 
represents  an  average  value  of  the  initial,  transient  infiltration 
flux.  The  model  simulations  represent  infiltration  from 
streamflow  events  ranging  in  duration  from  1  to  5.5  hours; 
long-term  constant  infiltration  rates  were  probably  not 
achieved  during  these  ephemeral  streamflows. 

Rates  of  Ephemeral-Stream  Channel  Recharge 

Recharge  deriving  from  ephemeral  streamflow  along  the 
length  of  a  channel  can  be  estimated  by: 

i 

R  =  J (Swqt  -  EV  -  T)d£  ,  (11) 

0 

where 

R  =  recharge  [L3], 

S  =  streamflow  duration  [t], 
w  =  wetted  width  of  channel  [L], 
q.  =  infiltration  flux  [L/t], 

EV  =  water  lost  to  channel  evaporation  per  channel 
length  [L3/L], 

T  =  water  lost  to  near-channel  transpiration  per 
channel  length  [L3/L],  and 
l  =  length  of  channel  containing  alluvial  sediments 
[L]. 


Average  annual  recharge  was  determined  by  dividing  the 
calculated  recharge  rate  for  the  2001-2002  period  by  the 
number  of  years  in  the  time  period  (2  years  for  this  study). 
The  total  streamflow  duration  was  determined  for  the 
2-year  period  at  each  streambed-temperature  sensor,  ARS 
flume,  and  USGS  streamflow-gaging  station  (table  8),  and 
an  exponential  best-fit  relation  between  streamflow  duration 
and  distance  in  the  channel  downstream  from  the  mountain 
front  was  then  developed  for  each  channel  (fig.  7).  Only  sites 
having  streamflow-duration  data  for  the  entire  2001-2002 
period  were  used.  The  width  of  the  channel  at  each  borehole 
was  known  (table  2),  and  was  linearly  extrapolated  between 
boreholes.  The  averaged  wetted-channel  width  during  flow 
events  was  determined  to  be  between  10  and  20  percent 
of  the  channel  width  on  the  basis  of  the  effective  wetted- 
channel  width  along  a  portion  of  Walnut  Gulch  as  estimated 
by  Goodrich  and  others  (2004).  An  average  wetted-channel 
width  of  15  percent  of  the  channel  width  was  used  in  the 
calculations.  Goodrich  and  others  (2004)  also  estimated 
channel  evaporation  losses  and  near-channel  transpiration 
losses  for  1999  and  2000  along  a  6.8-km  reach  of  Walnut 
Gulch.  The  average  losses  were  scaled  so  that  channel 
evaporation  was  represented  by  3.0-10'4hm3/km/2yr  and  near¬ 
channel  transpiration  was  represented  by  1.4-1 0'2  hm3/km/2yr. 
These  values  were  used  for  all  the  channels  in  the  2001-2002 
calculations 

Ephemeral- stream  channel  recharge  rates  determined 
with  equation  1 1  were  compared  with  rates  derived  using 
other  independent  methods  (channel-reach  water  balance, 
chloride  mass  balance,  ground-water  mounding  model, 
and  microgravity  measurements)  for  a  6.8-km  reach  of 
Walnut  Gulch  (Goodrich  and  others,  2004).  The  infiltration 
flux  used  in  this  comparison  was  the  lowest  one  calculated 
for  Walnut  Gulch  by  sediment-temperature  analysis.  As 
previously  stated,  streamflow  events  for  which  infiltration 
fluxes  could  be  determined  were  all  less  than  5.5  hours 
in  duration;  the  calculated  infiltration  fluxes  represent  an 
average  value  of  the  initial,  transient  infiltration  fluxes. 

The  lowest  calculated  infiltration  fluxes  are  assumed  to  best 
represent  the  steady- state  infiltration  rate.  Calculations  of 
recharge  for  1999  were  0.44  hm3  using  equation  11  and 
0.14  to  0.47  hm3  using  other  methods. Recharge  calculations 
for  2000  were  0.16  hm3  using  equation  11  and  0.13  to 
0.37  hm3  using  other  methods.  These  results  validate  the 
assumption  that  the  lowest  calculated  infiltration  flux  can 
represent  the  steady- state  infiltration  flux  when  determining 
ephemeral-stream  channel  recharge. 
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Figure  7.  Regression  of  streamflow  duration  to  distance  from  the 
mountain  front,  2001  and  2002,  Sierra  Vista  subwatershed,  Arizona. 


Ephemeral-Stream  Channel  Infiltration 

Infiltration  fluxes  during  streamflow  events  and 
percolation  fluxes  immediately  following  those  events  were 
determined  for  sites  in  Woodcutters  Wash,  Miller  Canyon 
Wash,  Walnut  Gulch,  Banning  Creek,  and  Greenbush  Draw. 
Data  used  in  the  determination  of  fluxes  included  physical, 
hydraulic,  and  thermal  properties  of  sediment;  chloride 
concentrations  of  sediment;  pore- water  stable-isotope  values 
and  tritium  activity;  neutron  soil-moisture  and  EMI  logs;  and 
sediment- temperature  and  water-level  monitoring.  Where 
possible,  fluxes  through  the  upper  1.5  to  12.2  m  of  sediments 
were  estimated. 

Woodcutters  Wash 

Woodcutters  Wash  is  a  northeast-trending  channel  with 
headwaters  in  the  Huachuca  Mountains.  The  channel  is 
23.0  km  long  and  has  an  average  slope  of  13.9  m/km.  The 
channel  widens  from  less  than  1  m  along  the  upper  reaches 
to  more  than  20  m  along  the  lower  reaches.  It  is  incised  less 
than  0.5  m  along  the  upper  reaches  and  2  m  or  more  along  the 
lower  reaches.  Several  knickpoints  occur  along  the  lower  one- 
half  of  the  channel  in  an  area  where  the  stream  has  cut  through 
several  layers  of  clay  and  caliche. 


The  channel  drains  a  3 7. 5 -km2  area  that  ranges  in 
elevation  from  1,204  to  2,488  m.  Desert  scrub  and  grassland 
are  the  primary  land-cover  types  in  the  drainage  basin; 
mesquite  is  dominant  along  the  mountain  fronts  and  oak 
woodland  is  dominant  in  the  mountainous  areas  (fig.  2). 

About  1  percent  of  the  drainage  basin,  the  central  part,  is 
urban.  The  banks  of  Woodcutters  Wash  have  been  stabilized 
in  this  area. 

Soils  in  the  drainage  basin  vary  from  loamy  soils 
near  the  mountains  to  calcium  carbonate  and  clayey  soils 
downgradient.  Along  the  stream  channel,  soils  are  sandy 
(fig.  3).  Thickness  of  permeable  stream  alluvium  is  about  1  m 
near  the  mountains,  a  few  centimeters  within  Sierra  Vista, 
and  about  3  m  in  the  lower  reaches  of  Woodcutters  Wash. 

The  stream  alluvium  is  underlain  by  the  compacted  clayey 
gravels  of  the  Pantano  Formation  near  the  mountains;  the 
poorly  sorted  sand,  silt,  and  clay  of  the  basin  fill  in  the  Sierra 
Vista  area;  and  the  low-permeability  silt  and  clay,  and  calcium 
carbonate  of  the  basin  fill  along  the  lower  reaches  of  the  wash. 

Ephemeral  streamflows  recorded  from  2001  to  2002  in 
Woodcutters  Wash  ranged  in  duration  from  about  15  min  to 
22  hr  (table  8).  Flows  in  the  upper  reaches  generally  were 
more  frequent  and  were  of  a  longer  duration  than  flows 
in  reaches  downstream  (fig.  7).  The  farthest  upstream 
temperature  sensor,  375010,  detected  19  flows,  which  had 
an  average  duration  of  6.5  hr.  The  farthest  downstream 
temperature  sensor,  at  site  WC1,  detected  18  flows,  which 
had  an  average  duration  of  4.0  hr. 

Borehole  WC4a  is  the  farthest  upstream  site  in 
Woodcutters  Wash  and  is  about  0.4  km  from  the  base  of  the 
Huachuca  Mountains.  The  site  is  in  a  reach  about  0.5  m  wide 
and  incised  less  than  1  m  into  terrace  deposits,  basin  fill,  and 
the  Pantano  Formation.  The  highly  permeable  sand  and  gravel 
stream  alluvium  extends  from  the  land  surface  to  a  depth  of 
3.4  m,  and  the  less  permeable  semiconsolidated  conglomerate 
of  the  Pantano  Formation  extends  from  3.4  to  12.5  m  (fig.  6A). 
Electrical  conductivities  of  the  stream  alluvium  and  the 
Pantano  Formation  ranged  from  about  10  to  30  mS/m,  and 
from  about  25  to  45  mS/m,  respectively  (fig.  6A).  A  neutron 
soil-moisture  log  was  not  recorded  at  this  site. 

Physical  and  hydraulic  properties  were  determined  for 
two  core  samples  from  WC4a  collected  at  depths  of  2.1  and 
11.3  m  (table  6).  For  the  stream  alluvium  and  the  Pantano 
Formation,  measured  saturated  hydraulic  conductivity  was 
6.4- 10'2  and  4.2- 10'4  cm/s,  respectively;  porosity  was  0.43  and 
0.33  cm3/cm3,  respectively;  and  soil- water  content  was 
0.09  and  0.07  cm3/cm3,  respectively. 

When  borehole  WC4a  was  drilled  in  January  2001, 
a  saturated  interval  was  penetrated  above  the  Pantano 
Formation  at  a  depth  of  2.0  to  3.4  m  (fig.  8A).  Sediment- 
chloride  concentrations  were  lower  than  8  mg/kg  in  the 
stream  alluvium  and  generally  lower  than  3  mg/kg  in  the 
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Pantano  Formation  (fig.  9 A).  Tritium  activity  in  the  stream 
alluvium  pore  water  was  similar  to  activity  in  present-day 
precipitation  (between  3.4  and  3.7  TU;  fig.  9 A).  Tritium 
activity  was  10.1  TU  in  the  upper  1.7  m  of  the  Pantano 
Formation  and  was  between  1  and  2  TU  from  a  depth  of  6.6  to 
11.1m.  This  pattern  indicates  that  recent  precipitation  has 
likely  infiltrated  through  the  stream  alluvium  to  the  Pantano 
Formation.  Water  is  infiltrating  into  the  semiconsolidated, 
unsaturated  Pantano  Formation  slowly,  at  a  rate  of  about 
0.9  cm/yr  (assuming  an  average  volumetric  soil- water  content 
of  0.20  cm3/cm3).  Pore  water  within  the  Pantano  Formation 
was  depleted  in  heavy  isotopes,  similar  to  winter  precipitation 
(fig.  10).  This  could  indicate  that  the  Pantano  Formation 
receives  infiltration  primarily  during  winter  streamflows  or 
that  the  existing  pore  water  was  derived  from  infiltration 
received  when  the  climate  was  cooler  than  at  present. 

Temperature  methods  detected  only  one  summer  flow 
event  (30  minutes  in  duration)  at  WC4a  between  January 
2001  and  December  2002  (table  8).  Water  levels  in  the 
saturated  zone,  however,  rose  at  least  six  times  during  the 
summer  of  2001  (fig.  8A).  The  lack  of  sediment- temperature 
signals  may  be  caused  by  a  retention  basin  about  260  m 
upstream  from  the  borehole.  Water-level  recovery  was  likely 
the  result  of  infiltration  at  the  retention  basin  and  subsequent 
downgradient  subflow  above  the  contact  between  the  Pantano 
Formation  and  the  stream  alluvium.  A  recharge  flux  could 
not  be  calculated  from  the  detected  flow  event  because  the 
water  level  began  to  rise  before  the  onset  of  temperature- 
detected  flow  in  the  channel. 


Borehole  WC3  is  about  7.2  km  downstream  from  the 
mountain  front.  This  site  is  within  the  city  limits  of  Sierra 
Vista,  just  downstream  from  an  area  where  the  banks  and 
streambed  of  Woodcutters  Wash  have  been  stabilized.  The  site 
is  in  a  reach  about  8  m  wide  and  incised  1.8  to  3.0  m  into 
terrace  deposits  and  basin  fill.  Stream  alluvium  in  the  channel 
bed  is  only  a  few  centimeters  thick.  The  poorly  sorted  sand, 
silt,  and  clay  of  the  upper  basin  fill  extends  from  the  land 
surface  to  a  depth  of  22.6  m  (fig.  6 B).  The  silt  and  clay 
content  of  the  upper  basin  fill  is  generally  30  to  60  percent. 
Electrical  conductivity  of  the  sediments  ranged  from  about 
25  to  60  mS/m  (fig.  6 B).  Neutron  soil-moisture  counts  ranged 
from  6,000  to  9,000  in  the  borehole  (fig.  6 B). 

Physical  and  hydraulic  properties  were  determined  for 
five  core  samples  from  WC3  collected  at  depths  of  0.9,  1.5, 

8.2,  13.7,  and  19.8  m  (table  6).  The  measured  saturated 
hydraulic  conductivity  of  the  basin  fill  ranged  from  1.7-1  O'4 
to  9.9-10'3  cm/s.  Porosity  ranged  from  0.33  and  0.40  cm3/cm3, 
and  soil- water  content  ranged  from  0.1 1  to  0.19  cm3/cm3. 

Sediment-chloride  concentrations  ranged  from  3  to 
70  mg/kg  from  the  surface  to  4.9  m,  and  were  below  4  mg/kg 
from  5.2  to  22.6  m  (fig.  9 B).  Stable-isotope  data  indicate  that 
the  high  chloride  concentrations  near  the  surface  are  not  the 
result  of  chloride  accumulation  due  to  evapotranspiration; 
stable-isotope  values  in  pore  water  at  0.2  m  had  evaporative 
signatures,  but  the  values  at  0.8  m  and  below  did  not  (fig.  10). 
The  high  chloride  concentrations  in  the  shallow  part  of 
the  unsaturated  zone  are  probably  indicative  of  chloride 
concentrations  in  runoff  from  Sierra  Vista.  Pore- water  tritium 
activity  values  from  the  surface  to  16.6  m  were  similar  to  those 
of  present-day  precipitation  (between  1.0  and  6.3  TU;  fig.  9 B). 


Figure  8.  Water  levels  and  streamflow  durations  from  temperature  sensors,  Woodcutters  Wash  boreholes,  Sierra  Vista 
subwatershed,  Arizona.  A,  WC4a;  B ,  WC1. 
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error  of  estimate  bars 
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Figure  9.  Sediment-chloride  and  pore-water  tritium  data.  Woodcutters  Wash  boreholes,  Sierra  Vista  subwatershed,  Arizona. 
A,  WC4a;  B,  WC3;  C,  WC2;  0,  WC1 . 
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Figure  10.  Stable-isotope  data  for  pore  water  collected  from  Woodcutters 
Wash  boreholes  and  for  precipitation.  Sierra  Vista  subwatershed,  Arizona. 


The  tritium  activity  at  19.7  m  was  0.4  TU.  Five  flow  events 
were  detected  at  WC3  between  January  2002  and  December 
2002  (table  8).  There  were  no  sediment-temperature 
perturbations,  however,  at  1.5  m  and  below  during  or  after 
flow  events. 

The  temperature  data  indicate  that  the  infiltration  flux  at 
WC3  is  too  low  to  perturb  the  sediment-temperature  profile 
during  and  after  individual  flow  events.  Tritium  data  indicate 
an  infiltration  rate  of  about  4.1  cm/yr  (assuming  an  average 
volumetric  soil- water  content  of  0.15  cm3/cm3).  Stable-isotope 
values  in  pore  water  were  between  the  ranges  for  summer 
and  winter  precipitation,  indicating  that  infiltration  may  have 
occurred  during  both  summer  and  winter  streamflow  events  in 
the  past  (fig.  10). 

Borehole  WC2  is  about  12.3  km  downstream  from  the 
mountain  front.  The  site  is  in  a  22-m  wide  channel  reach  that 
is  incised  about  1.8  m  into  terrace  deposits  and  basin  fill.  The 
channel  bed  is  floored  with  only  a  veneer  of  stream  alluvium 
and  is  underlain  by  less  permeable,  poorly  sorted  sand,  silt, 
and  clay  of  the  basin  fill  (fig.  6C).  The  silt  and  clay  content 
of  the  basin  fill  is  generally  50  to  80  percent.  Geophysical 
borehole  logs  were  not  recorded  in  this  borehole. 

Physical  and  hydraulic  properties  were  determined  for 
five  core  samples  from  WC2  collected  at  depths  of  2.1,  3.7, 

6.7,  8.2,  and  9.8  m  (table  6).  The  measured  saturated  hydraulic 
conductivity  of  the  basin  fill  ranged  from  9.8-10'6  to 
2.2- 10'2  cm/s.  Porosity  of  the  basin  fill  ranged  from  0.34  to 
0.52  cm3/cm3.  The  soil- water  content  of  the  basin  fill  ranged 
from  0.15  to  0.47  cm3/cm3. 

Sediment-chloride  concentrations  to  12.5  m  were 
low  (less  than  7  mg/kg),  and  tritium  activity  to  1 1.1  m  was 
comparable  to  that  of  present-day  precipitation  (between 
3.9  and  5.7  TU),  indicating  that  recent  infiltration  has  occurred 
to  at  least  11.1m  (fig.  9C).  Pore- water  stable-isotope  values 
at  2.0  m  had  a  slight  evaporative  signature  and  indicate 
infiltration  of  both  summer  and  winter  streamflow  (fig.  10). 


This  borehole  was  not  completed  with  a  casing;  therefore, 
sediment-temperature  data  are  not  available  and  infiltration 
fluxes  could  not  be  determined. 

Borehole  WC1  is  at  the  SPRNCA  boundary  16.4  km 
downstream  from  the  Huachuca  Mountains.  The  site  is  in  a 
reach  10  m  wide  and  incised  several  meters  into  the  basin 
fill.  The  highly  permeable  stream  alluvium,  consisting  of 
sand  and  gravel,  extends  from  land  surface  to  3.7  m;  basin  fill 
consisting  of  silt  and  clay  extends  from  3.7  to  9.1  m;  and  basin 
fill  consisting  of  sand  extends  from  9.1  to  9.4  m  (fig.  6D). 

The  basin  fill  includes  a  large  amount  of  calcium  carbonate 
from  3.7  to  5.2  m,  and  a  slightly  sandy  zone  from  6.4  to  7.0  m. 
Electrical  conductivity  of  the  sediments  was  about  10  mS/m 
within  the  stream  alluvium  and  60  to  100  mS/m  within  the 
basin  fill  (fig.  6 D).  The  higher  electrical  conductivity  in  the 
basin  fill  is  likely  due  to  both  higher  soil- water  content  and  a 
higher  silt  and  clay  content  than  the  stream  alluvium. 

Drilling  of  borehole  WC1  in  January  2001  penetrated 
saturated  sediments  at  7.5  m.  Water  levels  were  highly  variable 
during  the  monitoring  period,  from  6.6  to  8.3  m  below  the  land 
surface,  and  rose  as  much  as  0.5  m  in  2  days  following  a  single 
flow  event  (fig.  8#).  Neutron  soil-moisture  counts  ranged 
from  5,000  to  7,000  in  the  stream  alluvium  and  7,000  to 
nearly  20,000  in  the  basin  fill  (fig.  6 D).  The  highest  counts 
correspond  to  saturated  or  near- saturated  sediments. 

Physical  and  hydraulic  properties  were  determined  for 
four  core  samples  from  WC1  collected  at  depths  of  0.3,  2.1, 
3.7,  and  6.7  m  (table  6).  The  measured  saturated  hydraulic 
conductivity  of  the  stream  alluvium  ranged  from  3.6-1  O'2  to 
7.2- 10'2  cm/s  and  that  of  the  basin  fill  ranged  from 
2.1-10'6  to  3.2- 10'3  cm/s.  Porosity  of  the  stream  alluvium 
ranged  from  0.37  to  0.46  cm3/cm3;  porosity  of  the  basin 
fill  was  not  determined.  Soil-water  contents  of  the  stream 
alluvium  ranged  from  0.07  to  0.13  cm3/cm3,  and  the  soil- 
water  content  at  one  depth  in  the  basin  fill  was  0.37  cm3/cm3. 
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Sediment-chloride  concentrations  to  8.8  m  were  low  (less  than 
6  mg/kg)  and  tritium  activity  values  to  8.1  m  were  between 
3.9  and  7.2  TU  (fig.  9 D),  indicating  that  recent  infiltration  has 
occurred  to  at  least  9.1  m.  Stable-isotope  values  in  pore  water 
were  between  the  isotope  ratio  ranges  for  summer  and  winter 
precipitation,  indicating  that  infiltration  prior  to  drilling  may 
have  occurred  during  both  summer  and  winter  streamflow 
events  (fig.  10).  An  evaporative  signature  is  not  evident  in  the 
stable-isotope  data. 

Eighteen  flow  events  occurred  at  WC1  between  January 
2001  and  December  2002  (tab;e  8).  During  four  events, 
temperature  perturbations  were  detected  at  1.5  and  3.0  m,  but 
not  below  3.0  m  (table  9).  The  water  level  in  the  saturated 


zone,  however,  rose  sharply  following  each  event,  indicating 
that  the  infiltrating  streamflow  is  recharging  the  saturated 
zone.  The  lack  of  temperature  variation  below  3.0  m,  coupled 
with  the  substantial  water-level  response,  supports  the 
occurrence  of  recharge  from  nearby  areas  but  not  directly 
through  the  silt  and  clay  at  the  site.  The  water-level  response 
may  represent  a  confined  response  in  sandy  stringers  within 
the  silt  and  clay  that  are  hydraulically  connected  to  upstream 
areas  where  infiltration  occurs. 

Infiltration  fluxes  to  3.0  m  below  the  land  surface  during 
flow  events  ranged  from  1.3  to  1.5  m/hr  (fig.  11  and  tables  9 
and  11).  Percolation  fluxes  to  3.0  m  after  flow  events  ranged 
from  0.4  to  1.5  m/hr  (table  9). 


JULY  7,  2001 
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Figure  11.  Streamflow  duration  and  sediment  temperatures  from  sensors  and  from  VS2DH  model  output,  Woodcutters 
Wash  borehole  WC1,  Sierra  Vista  subwatershed,  July  7, 2001. 
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Miller  Canyon  Wash 

Miller  Canyon  Wash  is  a  northeast-trending  channel 
with  headwaters  in  the  Huachuca  Mountains.  The  channel 
is  16.2  km  long  and  has  an  average  slope  of  20.5  m/km.  The 
channel  widens  from  less  than  1  m  along  the  upper  reaches 
to  more  than  4  m  along  the  lower  reaches.  It  is  incised  less 
than  1  m  along  the  upper  reaches  and  1  m  or  more  along 
the  lower  reaches.  The  channel  drains  a  32. 8 -km2  area  that 
ranges  in  elevation  from  1,252  to  2,859  m.  Desert  scrub, 
grassland,  and  mesquite  are  the  primary  land-cover  types  in 
the  drainage  basin;  oak  woodland  and  forest  are  dominant  in 
the  mountainous  areas  (fig.  2).  About  15  percent  of  the  basin, 
along  the  mountain  front,  is  urban.  The  urban  areas,  however, 
generally  have  less  paved  area  than  the  urban  areas  in  the 
Woodcutters  Wash  drainage  basin. 

Soils  in  the  drainage  basin  vary  from  loamy  soils  and  fine 
paleosoils  below  the  mountains  to  clayey  soils  downgradient. 
Along  the  stream  channel,  soils  are  sandy  (fig.  3).  The 
permeable  stream  alluvium  ranges  in  thickness  from  a  few 
centimeters  near  the  mountains  to  a  few  meters  in  the  lower 
reaches.  The  stream  alluvium  is  underlain  by  low  permeability 
basin  fill  that  has  variable  amounts  of  clay,  silt,  sand,  and 
gravel. 

Ephemeral  streamflows  recorded  in  2001-2002  ranged  in 
duration  from  about  30  min  to  19.8  hr  (table  8).  Flows  in  the 
upper  reaches  were  less  frequent  but  generally  were  of  longer 
duration  than  those  farther  downstream  (fig.  7).  The  farthest 
upstream  site,  temperature  sensor  377806,  detected  four 
flows,  which  had  an  average  duration  of  9.1  min.  The  farthest 
downstream  site,  temperature  sensor  374919,  detected  six 
flows,  which  had  an  average  duration  of  2.8  hr. 

Borehole  MC2  is  at  the  contact  between  the  crystalline 
rocks  of  the  Huachuca  Mountains  and  the  basin  fill.  The  site 
is  in  a  reach  about  2  m  wide  and  incised  less  than  1  m  into 
the  basin  fill.  Drilling  penetrated  a  layer  of  predominantly 
poorly  sorted  gravel  and  sand  of  the  basin  fill  from  the  land 
surface  to  15.8m  that  included  a  layer  of  silt  and  clay  from 
6.4  to  7.6  m  (fig.  6 E).  Electrical  conductivity  of  the  sediments 
at  MC2  ranged  from  about  20  to  40  mS/m  (fig.  6E).  Borehole 


MCI  is  at  the  SPRNCA  boundary,  7.1  km  downstream  from 
the  mountain  front.  The  site  is  in  a  reach  about  4  m  wide  and 
incised  several  meters  into  the  basin  fill.  Drilling  penetrated 
poorly  sorted  sand,  silt,  and  clay  of  the  basin  fill  from  the 
land  surface  to  25.0  m  (fig.  6 F).  Electrical  conductivity  of 
the  sediments  ranged  from  about  20  to  100  mS/m  (fig.  6 F). 
Neutron  soil-moisture  counts  ranged  from  about  7,000  to 
11,000  for  MC2  (fig.  6E);  a  neutron  soil-moisture  log  was  not 
completed  in  MCI. 

Physical  and  hydraulic  properties  of  the  basin  fill  were 
determined  for  three  core  samples  from  MC2  collected  at 
depths  of  1.8,  7.0,  and  10.1  m;  and  for  four  core  samples  from 
MCI  collected  at  depths  of  0.9,  4.0,  8.5,  and  19.2  m  (table  6). 
The  measured  saturated  hydraulic  conductivity  ranged  from 

1.7- 10'2  to  3.2-10  6  cm/s  at  MC2  and  from  9.3-10'3  to 

1.8- 10'6  cm/s  at  MCI.  Porosities  ranged  from  0.36  to 

0.42  cm3/cm3  at  MC2  and  from  0.36  to  0.39  cm3/cm3  at  MCI. 
Soil- water  contents  ranged  from  0.18  to  0.30  cm3/cm3  at  MC2 
and  from  0.28  to  0.31  cm3/cm3  at  MCI. 

Sediment-chloride  concentrations  were  lower  than 
5  mg/kg  in  the  upper  15.8  m  of  MC2  and  in  the  upper  21.9  m 
of  MCI.  Tritium  activity  in  the  upper  9.9  m  of  MC2  and 
in  the  upper  22.1  m  of  MCI  were  comparable  to  activity 
in  present-day  precipitation  (1.9  to  7.8  TU;  fig  12).  Recent 
infiltration  has  occurred  to  at  least  9.9  m  at  MC2  and  22.1  m 
at  MCI.  Stable-isotope  values  for  the  pore  water  at  both  sites 
indicated  that  infiltration  occurred  during  summer  and  winter; 
pore  water  below  0.2  m  did  not  show  evaporative  signatures 
(fig.  13). 

One  streamflow  event,  1 8  hours  in  duration,  was 
detected  between  January  2002  and  December  2002  at 
MC2;  streamflow  was  not  detected  between  April  2002  and 
December  2002  at  MCI  (table  8).  There  were  no  sediment- 
temperature  perturbations  at  or  below  1.5  m  during  or  after 
the  flow  event  at  MC2.  The  temperature  data  suggest  that  the 
infiltration  flux  at  MC2  is  too  low  to  perturb  the  sediment- 
temperature  profile  during  and  after  individual  events. 
Infiltration  fluxes  could  not  be  calculated  from  sediment- 
temperature  profiles  at  either  site  because  of  a  lack  of  data. 
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Figure  12.  Sediment-chloride  and  pore-water  tritium  data.  Miller  Canyon  Wash  boreholes.  Sierra  Vista  subwatershed,  Arizona. 
A,  MC2;  B,  MCI. 


Figure  13.  Stable-isotope  data  for  pore  water  collected  from  Miller  Canyon 
Wash  boreholes  and  for  precipitation.  Sierra  Vista  subwatershed,  Arizona. 
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Walnut  Gulch 


4.  WG3 


TRITIUM,  IN  TRITIUM  UNITS 


Walnut  Gulch  is  a  west-trending  channel  with  headwaters 
in  the  Tombstone  Hills  and  the  Mule  and  Dragoon  Mountains. 
The  main  channel  is  about  26.5  km  long  and  has  an  average 
slope  of  12.1  m/km.  The  channel  widens  from  less  than  2  m 
along  the  upper  reaches  to  more  than  40  m  along  the  lower 
reaches.  It  is  incised  less  than  1  m  along  the  upper  reaches  and 
4  m  or  more  along  the  lower  reaches.  The  channel  drains  a 
152.1-km2  area  that  ranges  in  elevation  from  1,164  to  1,939  m. 
Desert  scrub,  grassland,  and  mesquite  are  the  primary  land- 
cover  types  in  the  drainage  basin  (fig.  2).  Approximately 

3  percent  of  the  Walnut  Gulch  Basin,  near  the  northern  extent 
of  the  Tombstone  Hills,  is  urbanized. 

Soils  in  the  drainage  basin  vary  from  clayey  near  the 
mountains  to  secondarily  developed  calcium  carbonate 
downgradient;  soils  are  sandy  along  the  stream  channel 
(fig.  3).  The  permeable  stream  alluvium  ranges  in  thickness 
from  about  1  m  near  the  mountains  to  15  m  in  the  lower 
reaches  of  the  channel.  The  stream  alluvium  is  underlain 
by  low-permeability  consolidated  to  semiconsolidated 
conglomerate  of  the  lower  basin  fill  near  the  mountains  and 
by  upper  basin  fill  in  the  lower  reaches  of  the  channel. 

The  ARS  has  monitored  ephemeral  streamflows  in 
Walnut  Gulch  since  the  1950s  by  using  a  series  of  flumes. 
Ephemeral  streamflows  recorded  during  2001-2002  at 
flumes  1,  6,  and  1 1  ranged  in  duration  from  12  min  to  82.2  hr 
(table  8).  Flows  generally  were  more  frequent  and  of  longer 
duration  near  the  mountain  front  than  in  areas  downgradient 
(fig.  7).  The  farthest  upstream  site,  flume  11,  had  12  flows, 
which  had  an  average  duration  of  1 1 . 1  hr.  The  farthest 
downstream  site,  flume  1,  had  eight  flows,  which  had  an 
average  duration  of  3.2  hr. 

Borehole  WG3  is  just  upstream  from  flume  1 1  and  is 
7.2  km  downstream  from  the  mountain  front.  The  site  is  in  a 
reach  about  2  m  wide  and  incised  less  than  1  m  into  the  basin 
fill.  Drilling  penetrated  unconsolidated  stream  alluvium  of 
sand  and  gravel  to  1.8  m,  and  semiconsolidated  conglomerate 
of  the  lower  basin  fill  from  1.8  to  13.1  m  (fig.  6G).  The 
electrical  conductivity  of  the  conglomerate  was  about  5  to 
10  mS/m  (fig.  6G).  A  neutron  soil-moisture  log  was  not 
recorded  in  this  borehole.  The  high  degree  of  induration  of 
the  conglomerate  and  basalt  prevented  the  collection  of  cores 
below  1.8  m,  and  physical  and  hydraulic  properties  were  not 
determined  for  the  stream  alluvium  above  1.8  m. 

Sediment-chloride  concentrations  were  less  than 

4  mg/kg  in  the  upper  1.8  m,  and  tritium  activity  in  the  upper 
0.8  m  was  comparable  to  activity  for  present-day  precipitation 
(6.5  to  8.5  TU;  fig.  14A),  indicating  that  recent  infiltration 
has  occurred  to  at  least  0.8  m.  Although  the  primary  porosity 
of  the  semiconsolidated  conglomerate  is  small,  infiltration 
can  occur  through  secondary  fractures.  The  stable-isotope 
values  of  the  pore  water  at  0.2  and  0.8  m  had  a  strong 
evaporative  signature  (fig.  15).  Sediment  temperatures  were 
not  monitored  at  this  site,  and  therefore,  infiltration  fluxes 
were  not  calculated. 
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Figure  14.  Sediment-chloride  and  pore-water  tritium  data. 
Walnut  Gulch  boreholes.  Sierra  Vista  subwatershed,  Arizona. 
A,  WG3;  B,  WG2;  C,  WG 1 . 
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Figure  15.  Stable-isotope  data  for  pore  water  collected  from  Walnut  Gulch 
boreholes  and  for  precipitation,  Sierra  Vista  subwatershed,  Arizona. 


Borehole  WG2  is  just  upstream  from  flume  6  and  is 

13.2  km  downstream  from  the  mountain  front.  The  site  is  in 
a  reach  about  13  m  wide  and  incised  several  meters  into  the 
basin  fill.  Drilling  penetrated  unconsolidated  stream  alluvium 
of  sand  and  gravel  to  3.7  m,  semiconsolidated  conglomerate  of 
the  lower  basin  fill  from  3.7  to  1 1.6  m,  and  basalt  from  1 1.6  to 

12.2  m  (fig.  6 H).  Electrical  conductivity  of  the  sediments 
was  less  than  5  mS/m  in  the  stream  alluvium  and  from  10  to 
15  mS/m  in  the  conglomerate  (fig.  6 H).  Neutron  soil-moisture 
counts  in  the  borehole  ranged  from  about  6,000  to  10,000. 

As  at  WG3,  the  high  degree  of  induration  of  the  conglomerate 
and  basalt  prevented  the  collection  of  cores  below  3.7  m,  and 
physical  and  hydraulic  properties  were  not  determined  for  the 
stream  alluvium  above  3.7  m. 

Sediment-chloride  concentrations  in  the  upper  3.4  m 
were  less  than  the  detection  limit  of  1.80  mg/kg,  and  tritium 
activity  in  the  upper  2.3  m  was  comparable  to  values  for 
present-day  precipitation  (4.7  to  6.7  TU;  fig.  14 B),  indicating 
that  recent  infiltration  has  occurred  to  at  least  2.3  m.  The 
stable-isotope  values  for  the  pore  water  at  0.5  and  0.8  m  had 
a  strong  evaporative  signature;  however,  the  pore  water  at 
2.3  m  did  not  (fig.  15).  Although  the  primary  porosity  of  the 
consolidated  conglomerate  is  small,  infiltration  can  occur 
through  secondary  fractures. 

Seven  flow  events  occurred  at  flume  6  between  May  2002 
and  December  2002  (table  8).  Borehole  WG2  was  damaged 
sometime  after  mid-July  2002,  and  therefore  only  sediment 
temperatures  from  the  flow  on  July  19,  2002,  can  be  used 
to  determine  infiltration  fluxes  at  this  site.  During  this  flow, 
a  temperature  perturbation  was  detected  at  a  depth  of  3.0  m 
(table  9).  The  temperature  sensor  at  1.5  m  malfunctioned 
during  this  time  and  all  data  were  lost.  The  calculated 
infiltration  fluxes  to  3.0  m  during  flow  ranged  from  0.4  to 
1.0  m/hr  (fig.  16A  and  tables  9  and  11). 


Borehole  WG1  is  just  upstream  from  flume  1  and 
25.6  km  downstream  from  the  mountain  front.  The  site  is  in 
a  reach  about  31m  wide  and  incised  several  meters  into  the 
basin  fill.  Drilling  penetrated  stream  alluvium  of  sand  and 
gravel  in  the  upper  14.5  m  underlain  by  basin  fill  that  is  20  to 
55  percent  silt  and  clay  (fig.  61).  Electrical  conductivity  of  the 
sediments  ranged  from  less  than  20  mS/m  within  the  stream 
alluvium  to  25  to  65  mS/m  within  the  basin  fill  (fig.  61). 

A  neutron  soil-moisture  log  was  not  recorded  in  this  borehole. 

Physical  and  hydraulic  properties  were  determined  for 
four  core  samples  from  WG1  collected  at  depths  of  3.0,  3.7, 
13.4,  and  22.6  m  (table  6).  The  data  indicate  little  variation  of 
properties  between  the  stream  alluvium  and  the  upper  basin 
fill.  The  measured  saturated  hydraulic  conductivity  of  the 
cores  ranged  from  2.0- 10'2  to  4.6- 102  cm/s;  porosity  ranged 
from  0.35  to  0.47  cm3/cm3;  and  soil- water  content  ranged  from 
0.10  to  0.25  cm3/cm3. 

Sediment-chloride  concentrations  in  the  upper  22.3  m 
were  less  than  3  mg/kg,  and  tritium  activity  in  the  upper 
22.4  m  was  comparable  to  present-day  precipitation  (2.9  to 
4.5  TU;  fig.  14C),  indicating  that  recent  infiltration  has 
occurred  to  at  least  22.4  m.  Stable-isotope  values  at  2.9  m  had 
a  slight  evaporative  signature  (fig.  15)  and  also  indicated  that 
water  in  the  unsaturated  zone  infiltrated  during  both  summer 
and  winter  streamflow  events. 

Three  flow  events  occurred  at  flume  1  between  April 
2002  and  December  2002  (table  8).  All  were  during  the 
summer  and  had  durations  from  12  min  to  about  10  hr. 

During  each  event,  sediment-temperature  perturbations 
in  WG1  were  detected  at  1.5,  3.0,  and  4.6  m  (table  9). 
Temperature  perturbations  were  detected  at  6.1  m  after  the 
end  of  each  flow  event  (table  9),  but  no  perturbations  were 
detected  below  6.1  m. 
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Figure  16.  Streamflow  durations  and  sediment  temperatures  from  sensors  and  from  VS2DH  model  output.  Walnut  Gulch 
boreholes.  Sierra  Vista  subwatershed,  Arizona.  A,  WG2,  July  19, 2002;  B,  WG1,  July  26-27, 2002;  C,  WG1,  August  4, 2002. 
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Calculated  infiltration  fluxes  to  4.6  m  below  the  land 
surface  during  flow  events  ranged  from  0.4  to  2.1  m/hr 
(figs.  15 B  and  15C  and  tables  9  and  11).  Percolation  fluxes  to 
6.1  m  after  flow  ended  ranged  from  0.1  to  1.0  m/hr  (table  9). 

Banning  Creek 

Banning  Creek  is  a  west- trending  channel  with 
headwaters  in  the  Mule  Mountains.  The  channel  is  8.6  km 
long  and  has  an  average  slope  of  15.6  m/km.  The  creek 
drains  a  41.3-km2  area  that  ranges  in  elevation  from  1,257  to 
2,217  m.  Desert  scrub,  grassland,  and  mesquite  are  the 
primary  land-cover  types  in  the  drainage  basin;  oak  woodland 
dominates  the  mountainous  areas  (fig.  2).  Less  than  1  percent 
of  the  basin  is  urbanized  land. 

Soils  in  the  drainage  basin  vary  from  lithic  rock  near 
the  mountains  to  calcium  carbonate  downgradient;  soils  are 
sandy  along  the  stream  channel  (fig.  3).  The  permeable  stream 
alluvium  ranges  in  thickness  from  less  than  1  m  near  the 
mountains  to  more  than  7  m  in  the  lower  reaches  of  the  creek. 
The  stream  alluvium  is  underlain  by  limestone  in  the  upper 
reaches  and  unconsolidated  basin  fill  of  clay,  silt,  sand,  and 
gravel  in  the  lower  reaches. 

Ephemeral  streamflows  recorded  during  2001-2002 
ranged  in  duration  from  1  hr  to  48  days  (table  8).  Flows 
in  the  upper  reaches  were  not  more  frequent  but  were  of  a 
longer  duration  than  flows  farther  downstream  (fig.  7).  The 
farthest  upstream  streamflow-gaging  station,  09470700,  and 
temperature  sensor,  377822,  detected  two  flows  during 
2001-2002.  The  flows  had  an  average  duration  of  42.9  days. 
The  farthest  downstream  temperature  sensor,  at  site  BC1, 
detected  eight  flows,  which  had  with  an  average  duration 
of  4.6  hr. 

Borehole  BC1  is  3.9  km  downstream  from  the  crystalline 
rocks  of  the  Mule  Mountains.  The  channel  at  the  site  is  about 
4  m  wide  and  is  incised  several  meters  into  the  basin  fill. 


Drilling  penetrated  unconsolidated  stream  alluvium  of  clay, 
silt,  sand,  and  gravel  from  the  land  surface  to  6.4  m  (fig.  67). 
Basin  fill,  which  includes  greater  amounts  of  clay  and  silt 
than  the  stream  alluvium,  was  penetrated  from  6.4  to  25.9  m. 
Electrical  conductivity  ranged  from  10  to  50  mS/m  in  the 
stream  alluvium  and  from  10  to  70  mS/m  in  the  basin  fill 
(fig.  67).  Neutron  soil-moisture  counts  ranged  from  6,000  to 
10,000  in  the  borehole  (fig.  67). 

Physical  and  hydraulic  properties  were  determined  for 
six  core  samples  from  BC1  collected  at  depths  of  2.1,  5.2, 

8.2,  14.3,  15.8,  and  20.4  m  (table  6).  The  measured  saturated 
hydraulic  conductivity  ranged  from  9.6-10'3  to  1.2-101  cm/s 
for  the  stream  alluvium  and  from  4.1  TO'7  to  8. 9 TO'3  cm/s  for 
the  basin  fill.  Porosity  of  the  two  units  ranged  from  0.25  to 
0.43  cm3/cm3,  and  soil-water  content  ranged  from  0.08  to 
0.29  cm33/cm3. 

Sediment-chloride  concentrations  in  the  upper 
25.9  m  were  less  than  12  mg/kg,  and  tritium  activity  in  the 
upper  23.3  m  were  comparable  to  values  for  present-day 
precipitation  (3.2  to  6.7  TU;  fig.  17),  indicating  that  recent 
infiltration  has  occurred  to  at  least  23.3  m.  Stable-isotope 
values  for  pore  water  were  between  bulk  precipitation  ranges 
for  summer  and  winter  (fig.  18),  indicating  that  infiltration 
occurs  during  summer  and  winter  streamflow  events. 

Eight  flow  events  occurred  at  BC1  between  January  2001 
and  December  2002  (table  8).  During  seven  of  the  events, 
temperature  perturbations  were  detected  at  1.5  and  3.0  m; 
temperature  perturbations  were  also  detected  at  4.6  and  6.1  m 
after  the  end  of  each  event  (table  9). 

Calculated  infiltration  fluxes  to  3.0  m  during  flow  events 
ranged  from  0.4  to  2.8  m/hr  (fig.  19  and  tables  9  and  11). 
Percolation  fluxes  to  6.1  m  after  flow  ended  ranged  from 
0.7  to  2.8  m/hr  (tables  9  and  11). 
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Figure  17.  Sediment-chloride  and  pore-water  tritium 
data.  Banning  Creek  borehole  BC1,  Sierra  Vista 
subwatershed,  Arizona. 


Figure  18.  Stable-isotope  data  for  pore  water  collected  from  Banning  Creek 
borehole  BC1  and  for  precipitation.  Sierra  Vista  subwatershed,  Arizona. 
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Figure  19.  Streamflow  durations  and  sediment  temperatures  from  sensors  and  from  VS2DH  model  output.  Banning  Creek 
borehole  BC1,  Sierra  Vista  subwatershed,  Arizona.  A,  July  20, 2001;  B ,  July  25, 2001;  C,  August  1, 2001. 
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Greenbush  Draw 

Greenbush  Draw  is  a  northwest-trending  channel  with 
headwaters  in  the  Mule  Mountains  and  in  parts  of  Mexico. 

The  main  channel  is  19.5  km  long  and  has  an  average  slope 
of  8.8  m/km.  The  channel  widens  from  less  than  2  m  along 
the  upper  reaches  to  more  than  14  m  along  the  lower  reaches 
of  the  draw.  It  is  incised  2  m  or  more  along  the  lower  reaches. 
The  draw  drains  a  3 06. 3 -km2  area  that  ranges  in  elevation  from 
1,280  to  2,258  m.  Grassland,  mesquite,  and  desert  scrub  are 
the  primary  land-cover  types  in  the  drainage  basin  (fig.  2). 
About  4  percent  of  the  basin  is  urbanized  land. 

Soils  in  the  drainage  basin  are  primarily  sandy  along  the 
stream  channels  and  clayey  and  calcium  carbonate  elsewhere 
(fig.  3).  Secondary  calcium  carbonate  soils  dominate  the 
lower  downgradient  areas  to  the  north  of  Greenbush  Draw, 
and  clayey  soils  dominate  the  downgradient  areas  to  the 
south  of  Greenbush  Draw.  The  stream  alluvium  ranges  in 
thickness  from  2  m  or  less  in  the  upper  reaches  to  10  m  along 
the  lower  reaches  and  is  underlain  by  unconsolidated  to 
semiconsolidated  upper  basin  fill  that  includes  interbedded 
sand,  gravel,  silt,  and  clay. 

Ephemeral  streamflows  recorded  during  2001-2002 
ranged  in  duration  from  30  min  to  19  hr  (table  8).  Flows  in 
the  upper  reaches  of  Greenbush  Draw  were  not  more  frequent 
or  of  a  longer  duration  than  flows  at  sites  farther  downstream 
(fig.  7).  The  farthest  upstream  temperature  sensor,  377815, 
detected  five  flows  in  2001  and  2002,  which  had  an  average 
duration  of  7.5  hr.  The  farthest  downstream  temperature 
sensor,  377820,  detected  six  flows,  which  had  an  average 
duration  of  10.0  hr. 

Borehole  GD2  is  5.6  km  downstream  from  the  limestone 
and  crystalline  rocks  of  the  mountains.  The  channel  at  the 
site  is  about  2  m  wide  and  is  incised  less  than  1  m  into  the 
basin  fill.  Drilling  penetrated  unconsolidated  stream  alluvium 
of  clay,  silt,  sand,  and  gravel  from  the  surface  to  2.0  m  and 
basin  fill  from  2.0  to  25.0  m  (fig.  6 K).  The  basin  fill  has 
greater  percentages  of  clay  and  silt  than  the  stream  alluvium. 
Electrical  conductivity  of  the  sediments  ranged  from  10  to 
70  mS/m  in  the  stream  alluvium  and  from  10  to  100  mS/m 
in  the  basin  fill  (fig.  6 K).  Thin  saturated  intervals  (less  than 
1.2  m  in  thickness)  were  penetrated  at  depths  of  8.2,  12.8,  and 

24.1  m  during  drilling.  A  neutron  log  was  not  recorded  in 
this  borehole. 

Physical  and  hydraulic  properties  were  determined  for 
three  core  samples  from  GD2  collected  at  depths  of  2.4, 

5.3,  and  13.1  m  (table  6).  All  the  cores  were  from  basin  fill. 
The  saturated  hydraulic  conductivity  ranged  from  5.9-1 0'3  to 
5.3- 10'2  cm/s.  Porosity  ranged  from  0.41  to  0.53  cm3/cm3,  and 
soil- water  content  ranged  from  0.14  to  0.44  cm3/cm3. 

Sediment  chloride  concentrations  in  the  upper  25.0  m 
were  less  than  3  mg/kg,  and  tritium  activity  in  the  upper 

19.1  m  was  similar  to  values  for  present-day  precipitation 
(3.5  to  4.3  TU;fig.  20A),  indicating  that  recent  infiltration 
has  occurred  to  at  least  19.1  m.  Stable-isotope  values  for 
pore  water  were  between  ranges  for  summer  and  winter 
precipitation  (fig.  21),  indicating  that  infiltration  occurs  during 
summer  and  winter  streamflow  events. 
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TRITIUM,  IN  TRITIUM  UNITS 


NOTE:  Tritium  values  include  standard 
error  of  estimate  bars 

Figure  20.  Sediment-chloride  and  pore-water  tritium  data, 
Greenbush  Draw  boreholes,  Sierra  Vista  subwatershed,  Arizona. 
A,  GD2;  B,  GDI. 
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Figure  21.  Stable-isotope  data  for  pore  water  collected  from  Greenbush  Draw 
boreholes  and  for  precipitation. 


Three  flow  events  occurred  at  GD2  between  January 
2002  and  December  2002  (table  8).  Temperature  perturbations 
were  detected  at  1.5,  3.0,  6.1,  9.1,  and  12.2  m  after  flow 
ended  (table  9).  Calculated  percolation  fluxes  after  flow  ended 
ranged  from  0.3  to  1.8  m/hr  (table  9.) 

Borehole  GDI  is  12.4  km  downstream  from  the  limestone 
and  crystalline  rocks  of  the  mountains.  The  channel  at  the 
site  is  about  14  m  wide  and  is  incised  several  meters  into  the 
basin  fill.  Drilling  penetrated  unconsolidated  stream  alluvium 
of  clay,  silt,  sand,  and  gravel  from  the  land  surface  to  10.0  m 
and  basin  fill  from  10.0  to  19. 5m  (fig.  6L).  The  basin  fill  has 
a  greater  percentage  of  clay  and  silt  than  the  stream  alluvium 
and  includes  a  large  amount  of  calcium  carbonate  from  10.1  to 
10.4  m.  Electrical  conductivity  of  the  sediments  ranged  from 
10  to  20  mS/m  in  the  stream  alluvium  and  from  15  to 
100  mS/m  in  the  basin  fill  (fig.  6L).  Neutron  counts  did  not 
vary  greatly  (7,000  to  12,000)  in  the  borehole  (fig.  6L). 

Physical  and  hydraulic  properties  were  determined 
for  four  core  samples  from  GDI  collected  at  depths  of  4.0, 

7.0,  and  13.1  m  (table  6).  The  measured  saturated  hydraulic 
conductivity  ranged  from  4.6- 10'3  to  1.4-10'2  cm/s  for  the 
stream  alluvium  and  was  6.8*  10'8  cm/s  for  the  basin  fill. 
Porosity  was  determined  only  for  one  core,  from  the  stream 
alluvium,  and  was  0.27  cm3/cm3.  Soil-water  content  was 
determined  only  for  the  stream  alluvium,  and  ranged  from 
0.11  to  0.15  cm3/cm3. 

Sediment-chloride  concentrations  in  the  upper  18.9  m 
were  less  than  the  detection  limit  of  1.80  mg/kg,  and  tritium 
activity  in  the  upper  19.1  m  was  similar  to  values  for  present- 
day  precipitation  (3.3  to  4.6  TU;  fig.  20/?),  indicating 
that  recent  infiltration  has  occurred  to  at  least  19.1  m. 
Stable-isotope  values  for  pore  water  at  0.5  m  had  strong 
evaporative  signatures;  the  remaining  values  were  between 


bulk  precipitation  ranges  for  summer  and  winter  (fig.  21), 
indicating  that  infiltration  occurs  during  summer  and  winter 
streamflow  events. 

One  flow  event  was  detected  at  GDI  between  May  2002 
and  December  2002  (table  8).  Temperature  perturbations,  at 
depths  of  1.5,  3.0,  6.1,  and  9.1  m,  were  detected  only  after 
flow  ended  (table  9).  Percolation  fluxes  to  9.1  m  after  flow 
ended  ranged  from  0.1  to  0.8  m/hr  (table  9). 

Basin-Floor  Infiltration 

Infiltration  fluxes  were  determined  for  basin-floor  sites 
in  the  Walnut  Gulch,  Carr  Canyon  Wash,  and  Woodcutters 
Wash  drainage  basins.  Data  used  in  the  determination  of 
infiltration  fluxes  included  physical  and  hydraulic  properties 
of  sediments,  chloride  concentrations  of  sediments,  and 
pore-water  stable-isotope  values  and  tritium  activity.  Where 
possible,  infiltration  fluxes  through  the  upper  4.0  to  9.8  m  of 
sediments  were  estimated. 

Borehole  WG4  is  within  the  Walnut  Gulch  drainage 
basin  on  a  hilltop  in  the  ARS  Lucky  Hills  research  station. 
Mesquite  woodland  is  the  primary  land-cover  type  at  the 
site  (fig.  2).  Basin  fill  penetrated  by  the  borehole  consists  of 
semiconsolidated  to  unconsolidated  sand  and  gravel  from  the 
surface  to  3.4  m,  consolidated  clay  and  silt  from  3.4  to  5.5  m, 
semiconsolidated  to  unconsolidated  sand  and  gravel  from 
5.5  to  9.8  m,  and  semiconsolidated  to  consolidated  clay,  silt, 
and  sand  from  9.8  to  13.4  m  (fig.  6 M).  Electrical  conductivity 
of  the  sediments  ranged  from  0  to  30  mS/m  (fig.  6 M).  A 
neutron  soil-moisture  log  was  not  recorded  in  this  borehole. 
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Physical  and  hydraulic  properties  were  determined  for 
two  core  samples  collected  at  depths  of  4.3  and  13.4  m  in 
the  borehole  (table  6).  Saturated  hydraulic  conductivity  was 
1.3-10'2  and  1.7-10'2  cm/s.  Porosity  was  0.48  and  0.50  cm3/cm3, 
and  soil- water  content  was  0.10  and  0.19  cm3/cm3. 

Sediment-chloride  concentrations  increased  sharply  to 
226  mg/kg  from  1.5  to  2.1  m  as  a  result  of  evapotranspirative 
enrichment  (fig.  22A).  Concentrations  then  decreased 
sharply  to  34  mg/kg  from  2.4  to  3.4  m,  probably  as  a  result 
of  no  downward  water  movement  or  preferential  transport 
past  the  root  zone.  Below  3.4  m,  the  chloride  concentration 
ranged  from  3  to  90  mg/kg.  Chloride  below  2.1  m  likely  was 
emplaced  by  infiltration  during  a  wetter  period. 

Tritium  activity  of  the  pore  water  indicates  that  water 
has  infiltrated  to  at  least  7.5  m  since  1952  (fig.  22A).  No 
anthropogenic  tritium  was  detected  at  13.3  m.  The  location  of 
the  tritium  peak  at  2.0  m  indicates  an  infiltration  flux  of  about 
1.0  cm/yr  (assuming  an  average  volumetric  soil- water  content 
of  0.19  cm3/cm3)  to  2.0  m.  The  chloride  data  indicate  that 
about  694  years  is  required  for  water  to  infiltrate  to  this  depth, 
at  a  flux  of  about  0.2  cm/yr  (fig.  22 B).  Deeper  penetration  of 
tritium  relative  to  that  of  chloride  was  previously  identified 
in  desert  soils;  possible  explanations  have  included  the 
downward  movement  of  tritium  in  the  vapor  phase  (Phillips 
and  others,  1988;  Scanlon,  1991)  and  the  downward 
preferential  flow  of  tritium  through  cracks  and  macropores 
(Izbicki  and  others,  2000). 

Pore- water  stable-isotope  values  to  13.3  m  plot  along 
an  evaporative  trend  line  that  has  a  slope  of  2.5  (fig.  23). 

The  point  of  intersection  between  the  evaporative  trend  line 
and  the  meteoric  water  line  indicates  that  historical  infiltration 
at  WG4  was  similar  isotopically  to  recent  winter  precipitation. 


A.  Chloride  and  tritium 

TRITIUM,  IN  TRITIUM  UNITS 


The  chloride,  tritium,  and  stable-isotope  data  indicate  that 
little  to  no  recent  infiltration  has  occurred  past  shallow  depths 
at  site  WG4. 

Borehole  WG5  also  is  within  the  Walnut  Gulch  drainage 
basin.  The  borehole  is  on  a  hilltop  within  the  ARS  Kendall 
research  station.  Grassland  is  the  primary  land-cover  type 
at  the  site  (fig.  2).  Drilling  penetrated  consolidated  basin- 
fill  conglomerate  from  the  surface  to  12.2  m  (fig.  6N).  The 
electrical  conductivity  ranged  from  10  to  20  mS/m  (fig.  6N). 

A  neutron  soil-moisture  log  was  not  recorded  in  this  borehole. 
The  degree  of  consolidation  of  the  sediments  precluded 
collection  of  complete  cores;  therefore,  physical  and  hydraulic 
properties  of  the  conglomerate  were  not  determined. 

Sediment-chloride  concentrations  were  not  measured  at 
this  site.  Stable-isotope  values  of  pore  water  at  0.2  m  showed 
a  slight  evaporative  signature,  but  were  not  as  enriched  as  at 
similar  depths  at  other  basin-floor  sites  (fig.  23).  The  tritium 
activity  of  pore  water  at  0.2  m  was  0.4  TU.  These  isotopic 
dataindicate  that  recent  infiltration  at  this  site  has  not  extended 
past  shallow  depths. 

Borehole  BF1  is  within  the  Carr  Canyon  drainage  basin 
on  a  flat  area  of  the  basin  floor  between  the  Carr  Canyon 
and  Ramsey  Canyon  Washes.  Vegetation  at  the  site  is  sparse 
and  the  primary  land-cover  types  are  mesquite  woodland  and 
grassland  (fig.  2).  Drilling  penetrated  two  primary  lithologic 
layers  within  the  basin  fill:  unconsolidated  sand  and  gravel 
from  the  surface  to  about  8  m;  and  semiconsolidated  clay,  silt, 
sand,  and  gravel  from  8  m  to  15.8m  (fig.  60).  The  electrical 
conductivity  of  the  sediments  were  about  15  to  25  mS/m  from 
the  surface  to  1 1  m,  and  about  30  to  60  mS/m  from  1 1  to 
15.8  m.  The  higher  values  below  11m  are  indicative  of  a 
higher  soil-water  content  or  higher  silt  and  clay  content  below 
11m  (fig.  60).  A  neutron  soil-moisture  log  was  not  recorded 
in  this  borehole. 


B.  Residence  time  and  infiltration  flux 

INFILTRATION  FLUX,  IN  CENTIMETERS  PER  YEAR 


Figure  22.  Data  for  Walnut  Gulch  borehole  WG4,  Sierra  Vista  subwatershed,  Arizona.  A,  Sediment-chloride  and  pore-water 
tritium  data;  B,  Residence  time  and  infiltration  flux,  calculated  from  chloride  data. 
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Figure  23.  Stable-isotope  data  for  pore  water  collected  from  Walnut  Gulch 
basin-floor  boreholes  WG4  and  WG5  and  for  precipitation,  Sierra  Vista 
subwatershed,  Arizona. 


Physical  and  hydraulic  properties  were  determined  for 
seven  core  samples  collected  from  BF1  at  depths  of  0.3,  4.0, 
5.5,  7.0,  8.5,  10.1,  and  13.1  m  (table  6).  Saturated  hydraulic 
conductivity  at  0.3  and  7.0  m  was  2.1-10'2  and  1.4-1  O'2  cm/s, 
respectively.  Porosity  was  0.39  and  0.35  cm3/cm3,  and  soil- 
water  content  was  less  than  0.02  cm3/cm3  and  0.1 1  cm3/cm3, 
respectively.  The  silt  and  clay  content  from  0.3  to  7.0  m 
was  less  than  15  percent.  The  silt  and  clay  content  at 
8.5  m,  however,  was  30.7  percent.  The  saturated  hydraulic 
conductivity  of  this  core  was  5.1-10'3  cm/s,  the  porosity  was 
0.45  cm3/cm3,  and  the  soil-water  content  was  0.26  cm3/cm3. 

The  silt  and  clay  content  from  10.1  and  31.1m  was  less  than 
15  percent. 

Sediment-chloride  concentrations  were  less  than  5  mg/kg 
from  the  surface  to  5.8  m  (fig.  24A)  and  were  between  10  and 
32  mg/kg  from  6.1  to  14.3  m.  Below  14.3  m,  concentrations 
were  less  than  10  mg/kg.  This  chloride  profile  is  not  typical 
for  desert  soils  because  the  chloride  peak  is  not  near  the  land 
surface;  it  is  broad,  extending  almost  8  m;  and  it  has  a  low 
concentration. 

Tritium  activity  of  the  pore  water  indicates  that  water  has 
infiltrated  to  at  least  5.3  m,  and  possibly  to  10.2  m,  since  1952 
(fig.  24A).  An  infiltration  flux  could  not  be  calculated  using 
the  tritium  data  because  there  is  not  a  clear  depth  of  maximum 
tritium  activity.  The  presence  of  tritium  to  5.3  m,  however, 
corresponds  well  with  the  chloride  profile,  which  indicates  that 
about  75  years  is  required  for  water  to  infiltrate  to  this  depth. 

Pore- water  stable-isotope  values  show  strong  evaporative 
enrichment  at  0.5  m  and  a  slight  evaporative  enrichment  from 
3.8  to  13.0  m;  values  plot  along  an  evaporative  trend  line  that 
has  a  slope  of  2.0  (fig.  25).  The  point  of  intersection  between 
the  evaporative  trend  line  and  the  meteoric  water  line  indicates 
that  historical  infiltration  at  BF1  was  similar  isotopically  to 
recent  winter  precipitation. 


The  chloride,  tritium,  and  stable-isotope  data  indicate 
that  infiltration  occurs  at  a  rate  of  about  6  cm/yr  through  the 
more  permeable  upper  zone  of  sediments  to  a  depth  of  about 
8  m  (fig.  24 B).  Infiltration  is  much  slower  through  the  less 
permeable  sediments  at  depths  of  8  to  15.8  m.  Infiltrated 
water  maintains  a  low  chloride  content  and  high  tritium 
activity  through  the  upper,  permeable  sediments  and  transports 
chloride  to  the  deeper  zone  as  it  moves  slowly  downward 
through  the  less  permeable  sediments.  The  chloride  content 
in  the  deep  zone  is  enriched  through  transport  of  water  to  the 
atmosphere  by  plant  roots. 

Borehole  WC4c  is  within  the  Woodcutters  Wash  drainage 
basin.  The  borehole  is  about  150  m  south  of  borehole  WC4a 
on  an  area  of  the  basin  floor  that  slopes  northward  toward 
Woodcutters  Wash  (fig.  26).  Grassland  is  the  primary  land- 
cover  type  at  the  site  (fig.  2).  Drilling  penetrated  terrace 
deposits  of  mostly  unconsolidated  sand  and  gravel  from 
the  surface  to  5.8  m;  basin  fill  of  clay,  silt,  sand,  and  gravel 
from  5.8  to  7.6  m;  and  a  semiconsolidated  conglomerate,  the 
Pantano  Formation,  from  7.6  to  12.5  m  (fig.  6P).  Electrical 
conductivity  of  the  terrace  and  basin  fill  ranged  from  less  than 
10  to  50  mS/m,  and  the  electrical  conductivity  of  the  Pantano 
Formation  ranged  from  50  to  75  mS/m  (fig.  6 P).  A  neutron 
soil-moisture  log  was  not  recorded  in  this  borehole. 

Physical  and  hydraulic  properties  were  determined  for 
five  core  samples  from  WC4c  collected  from  depths  of  1.1, 

2.4,  3.8,  9.9,  and  11.3  m  (table  6).  The  saturated  hydraulic 
conductivity  ranged  from  5.1-10'3  to  1.5-1  O'2  cm/s  for  the 
terrace  deposits  and  from  7.2- 10'6  to  1 .9*  10 5  cm/s  for  the 
Pantano  Formation.  Porosity  ranged  from  0.33  to  0.40  cm3/cm3 
for  both  units.  Soil-water  content  ranged  from  0.07  to 
0.15  cm3/cm3  for  the  terrace  deposits  and  was  0.24  in  the 
Pantano  Formation. 
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A.  Chloride  and  tritium 


B.  Residence  time  and  infiltration  flux 


TRITIUM,  IN  TRITIUM  UNITS 


CHLORIDE,  IN  MILLIGRAMS  PER  KILOGRAM 


INFILTRATION  FLUX,  IN  CENTIMETERS  PER  YEAR 


Figure  24.  Data  for  Carr  Canyon  borehole  BF1,  Sierra  Vista  subwatershed,  Arizona.  A,  Sediment-chloride  and  pore-water  tritium  data; 
B,  Residence  time  and  infiltration  flux,  calculated  from  chloride  data. 


Figure  25.  Stable-isotope  data  for  pore  water  collected  from  basin- 
floor  boreholes  BF1  and  WC4c  and  for  precipitation.  Sierra  Vista 
subwatershed,  Arizona. 
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Figure  26.  Generalized  hydrogeologic  section  in  the  area  of  Woodcutters  Wash  boreholes  WC4a  and  WC4c,  Sierra  Vista 
subwatershed,  Arizona. 


During  drilling  in  January  2001,  the  borehole  intersected 
a  saturated  interval  within  the  terrace  deposits  and  basin  fill 
above  the  Pantano  Formation  at  4.7  m  (fig.  27).  Sediment- 
chloride  concentrations  were  17  mg/kg  at  1.5  m  and  generally 
were  less  than  6  mg/kg  below  2.7  m  (fig.  28A).  Tritium 
activity  of  2.6  to  7.0  TU  in  pore  water  to  3.7  m  indicated 
that  recent  water  is  present  in  the  unsaturated  zone  above 
the  saturated  interval  (fig.  28A).  The  chloride  data,  however, 
indicate  that  about  150  years  is  required  for  water  to  infiltrate 
to  the  saturated  interval,  at  a  flux  of  about  3  cm/yr  (fig.  28#). 
Similar  to  conditions  at  site  WC4a,  the  saturated  interval 
at  WC4c  is  probably  controlled  by  infiltration  from  an 


upgradient  retention  basin  and  subsequent  downgradient 
subflow  above  the  contact  of  the  Pantano  Formation  with 
terrace  and  basin  fill. 

Pore- water  tritium  activity  indicates  that  the  perched 
water  is  infiltrating  into  the  semiconsolidated  Pantano 
Formation  very  slowly  at  about  0.4  cm/yr  (assuming  an 
average  soil- water  content  of  0.07  cm3/cm3).  This  flux  is 
comparable  to  the  percolation  flux  for  the  Pantano  Formation 
at  WC4a  of  0.3  cm/yr.  Pore-water  stable-isotope  values  above 
the  water  table  and  within  the  Pantano  Formation  did  not 
have  an  evaporative  signature,  but  rather  fell  within  the  range 


Figure  27.  Water  level  in  Woodcutters  Wash  borehole  WC4c,  Sierra  Vista  subwatershed,  Arizona. 
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A.  Chloride  and  tritium  B.  Residence  time  and  infiltration  flux 

TRITIUM,  IN  TRITIUM  UNITS  INFILTRATION  FLUX,  IN  CENTIMETERS  PER  YEAR 


Figure  28.  Data  for  Woodcutters  Wash  borehole  WC4c,  Sierra  Vista  subwatershed,  Arizona.  A,  Sediment-chloride  and  pore-water 
tritium  data;  B,  Residence  time  and  infiltration  flux,  calculated  from  chloride  data. 


of  values  for  winter  precipitation  (fig.  25).  This  indicates 
that  infiltration  at  WC4c  could  be  dominated  by  winter 
precipitation. 

Subwatershed  Ephemeral-Stream 
Channel  and  Basin-Floor  Recharge 

Average  annual  recharge  during  2001-2002  was 
estimated  for  the  Sierra  Vista  subwatershed  by  upscaling 
the  infiltration  fluxes  determined  for  the  basin-floor  and 
ephemeral-stream  channel  sites.  This  recharge  estimate 
does  not  include  mountain-front  recharge;  this  study  did  not 
monitor  streamflow  or  infiltration  near  the  contact  between 
the  hard  rock  of  the  mountains  and  the  basin  fill.  Precipitation 
during  the  period  of  study  was  slightly  below  average  and 
primarily  occurred  from  July  through  September  (fig.  4). 

As  a  result,  the  calculated  annual  recharge  represents  recharge 
under  drier- than-normal  conditions. 

At  the  basin-floor  sites  in  the  Walnut  Gulch  drainage 
basin,  the  long-term  infiltration  fluxes  were  very  low  (less 
than  1  cm/yr).  Chloride,  tritium,  and  stable-isotope  data 
indicate  long  periods  of  no  net  deep  downward  percolation 
flux  beneath  the  basin  floor.  Previous  investigators  have  shown 
that  infiltration  of  precipitation  through  the  basin  floor  and 
past  the  root  zone  may  not  be  significant  in  arid  areas  of  the 
Southwestern  United  States  that  have  thick  unsaturated  zones 
(Phillips,  1994;  Walvoord  and  others,  2002).  Additionally, 
Scott  and  others  (2000)  reported  that  moisture  fronts  in  the 
unsaturated  zone  at  WG4  and  WG5  reach  a  maximum  depth 
of  only  0.5  m  in  the  summer  and  only  1.3  m  in  the  winter. 


The  chloride,  tritium,  and  stable-isotope  data,  coupled  with  the 
work  done  by  previous  investigators,  support  the  conclusion 
that  recharge  through  the  basin  floor  area  is  negligible. 
Infiltration  rates  could  also  be  negligible  in  other  areas  of  the 
basin  floor  that  have  soil  and  vegetation  that  are  similar  to 
those  at  WG4  and  WG5. 

Infiltration  fluxes  were  higher  at  the  basin-floor  site  in  the 
Carr  Canyon  Wash  drainage  basin  (2  to  6  cm/yr)  than  at  other 
basin-floor  sites.  Chloride,  tritium,  and  stable-isotope  data 
indicate  active  infiltration  to  8  m  and  a  decrease  in  infiltration 
below  8  m.  The  change  in  the  infiltration  rate  below  8  m  is 
controlled  by  an  increase  in  the  silt  and  clay  content  of  the 
soil.  Relatively  high  infiltration  rates  could  occur  in  other 
areas  of  the  basin  floor  that  have  soil  and  vegetation  that  are 
similar  to  those  at  BF1.  The  depth  that  infiltrating  water  moves 
through  the  unsaturated  zone,  and  whether  the  infiltrating 
water  will  ultimately  become  recharge,  however,  are  strongly 
dependant  on  the  silt  and  clay  content  of  limiting  soil  horizons 
in  the  unsaturated  zone. 

Recharge  through  ephemeral- stream  channels  is  strongly 
influenced  by  the  thickness  and  hydraulic  properties  of 
stream  alluvium  and  the  permeability  of  the  underlying 
basin  fill.  Streamflow  infiltration  readily  occurs  through  the 
highly  permeable  stream  alluvium  and  the  unconsolidated 
basin  fill  that  do  not  have  thick  intervals  of  silt  and  clay. 
Infiltration  is  limited  by  shallow  low-permeability  layers, 
such  as  bedrock,  silt  and  clay,  caliche,  or  consolidated  basin 
fill.  Where  the  contact  between  unconsolidated  sediments 
and  low-permeability  layers  is  deeper  than  the  depth  of  root 
penetration,  streamflow  infiltrates  to  depths  at  which  water 
is  probably  not  removed  by  evapotranspiration.  This  water 
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will  remain  in  storage  within  the  unsaturated  zone  until  it  is 
vertically  displaced  by  subsequent  infiltration  and  eventually 
recharges  the  regional  aquifer.  Where  the  contact  between 
unconsolidated  sediments  and  low-permeability  layers  is 
shallow,  however,  aquifer  recharge  may  be  minimal  because 
the  infiltrated  water  is  available  for  upward  transport  by 
evapotranspiration.  Water  not  lost  by  evapotranspiration  likely 
flows  laterally  above  the  contact  until  it  encounters  the  limit 
of  the  low-permeability  layer,  reaches  permeable  fractures,  or 
recharges  the  aquifer  where  the  water  table  intersects  the  top 
of  the  low-permeability  layer. 

Infiltration  fluxes  at  the  ephemeral- stream  channel 
sites  were  very  high  (0.4  to  2.8  m/hr)  during  streamflow 
events.  Percolation  fluxes  were  also  high  (0.1  to  2.8  m/hr) 
immediately  following  streamflow  events.  Chloride  data 
indicate  that  in  areas  where  low-permeability  layers  are  absent 
in  the  shallow  subsurface,  infiltrated  water  moves  past  depths 
where  it  is  available  for  evapotranspiration,  resulting  in  a  net 
downward  water  flux.  Tritium  data  indicate  that  water  that  has 
infiltrated  since  1952  is  reaching  depths  of  more  than  20  m. 
We  assume  that  the  net  water  flux  below  20  m  continues  to  be 
downward  and  that  infiltrated  streamflow  ultimately  recharges 
the  regional  aquifer.  Stable-isotope  data  indicate  that  both 
summer  and  winter  streamflow s  contribute  to  ephemeral- 
stream  channel  recharge. 

Sub  watershed  recharge  during  2001-2002  from 
ephemeral-stream  channel  infiltration  was  estimated  by 
upscaling  the  calculated  infiltration  fluxes  to  the  area  covered 
by  each  ephemeral- stream  channel  in  the  watershed  and 
weighting  the  fluxes  by  streamflow  duration,  evaporation,  and 
transpiration  for  the  2-year  period,  per  equation  11  (table  12 
and  fig.  29).  Recharge  estimates  were  not  derived  for  channels 
in  the  Babocomari  River  drainage  basin  (group  1)  because  the 
river  is  perennial  over  much  of  its  length  and,  therefore,  use  of 
the  infiltration  parameters  derived  for  the  ephemeral- stream 
channels  is  not  appropriate. 

The  infiltration  flux  used  in  equation  1 1  was  the  lowest 
flux  calculated  for  each  channel  by  sediment-temperature 
analysis.  The  calculated  infiltration  fluxes  represent  average 
values  of  the  initial  transient  infiltration  fluxes  during 
each  ephemeral-flow  event.  Consequently,  we  assume  that 
the  lowest  calculated  infiltration  flux  best  represents  the 
steady- state  infiltration  flux  when  determining  ephemeral- 
stream  channel  recharge.  This  assumption  was  validated 
by  comparing  ephemeral- stream  channel  recharge  rates 
determined  by  use  of  equation  1 1  with  rates  determined 
by  using  other  independent  methods  (channel-reach  water 
balance,  chloride  mass  balance,  ground- water  mounding 
model,  and  microgravity  measurements)  for  a  6.8-km  reach  of 
Walnut  Gulch  (Goodrich  and  others,  2004). 


Groups  3,  5,  6,  and  7  contain  ephemeral- stream  channels 
for  which  infiltration  fluxes  were  determined.  We  assume  that 
these  infiltration  fluxes  are  the  same  in  other  ephemeral- stream 
channels  within  each  group.  The  infiltration  flux  determined 
for  Woodcutters  Wash  was  also  used  for  channels  in 
groups  2  and  4,  which  did  not  contain  instrumented  borehole 
sites.  Infiltration  fluxes  for  periods  of  streamflow  were 
not  calculated  for  Greenbush  Draw;  therefore,  the  average 
percolation  flux  for  periods  after  cessation  of  streamflow  was 
used  in  the  Greenbush  Draw  calculation.  Channel  segments 
in  Woodcutters  Wash  and  Walnut  Gulch  underlain  at  shallow 
depths  by  the  low-permeability  sediments  of  the  lower  basin 
fill  were  not  included  in  the  recharge  calculations;  it  was 
assumed  that  similar  conditions  exist  in  Coyote  Wash,  because 
of  its  location  relative  to  Woodcutters  Wash,  Government 
Draw,  and  High  Knolls  Canyon  Wash,  because  of  their 
locations  relative  to  Walnut  Gulch  (fig.  29). 

The  streamflow  durations  measured  in  Woodcutters  Wash 
during  2001-2002  were  greater  than  those  measured  in  Miller 
Canyon  Wash  (table  8  and  fig.  7);  the  longer  durations  are 
assumed  to  be  related  to  the  corresponding  higher  percentage 
of  urban  land  cover  in  the  drainage  basin.  The  streamflow- 
duration  and  wetted-channel- width  regressions  determined 
for  Woodcutters  Wash  were  also  used  for  Coyote  and  Carr 
Canyon  Washes,  whose  drainage  basins  have  similar  amounts 
of  urban  land  cover  as  that  in  Woodcutters  Wash  and,  likely, 
have  similar  runoff  characteristics.  The  streamflow-duration 
and  wetted-channel-width  regressions  determined  for  Miller 
Canyon  Wash  were  used  for  the  remaining  channels  in 
groups  2,  3,  and  4. 

The  estimated  annual  recharge  rate  derived  from 
ephemeral-stream  channel  infiltration  during  2001-2002  is 
2.95  hm3/yr  (table  12).  This  rate  is  slightly  lower  than  average 
recharge  rates  of  3.4  to  9.7  hm3/yr  determined  by  Goodrich 
and  others  (2004)  for  1999-2000.  The  rate  calculated  during 
this  study  is  lower  probably  because  fewer  streamflow  events 
occurred  during  the  2001-2002  monitoring  period  than  during 
1999-2000.  At  the  Tombstone  and  Huachuca  Mountains  gages, 
precipitation  during  2001-2002  averaged  36.4  cm/yr,  whereas 
precipitation  during  1999-2000  averaged  47.4  cm/yr. 

The  ephemeral-stream  channel  recharge  rate  of 
2.95  hm3/yr  is  about  12  to  19  percent  of  previous  total 
recharge  estimates  for  the  Sierra  Vista  sub  watershed  of  15  to 
23  hm3/yr  (Freethey,  1982;  Corell  and  others,  1996;  Goode 
and  Maddock,  2000),  indicating  that  ephemeral- stream 
channel  recharge  is  a  significant  component  of  total  recharge. 
These  previous  investigators  have  assumed  that  recharge 
occurred  only  within  about  1.5  km  of  the  mountains;  however, 
the  borehole  data  clearly  indicates  that  when  streamflow  is 
available,  a  significant  portion  of  recharge  is  distributed  along 
the  entire  lengths  of  ephemeral-stream  channels  (fig.  29). 
During  this  study,  most  of  the  streamflow,  and  therefore 
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Table  12.  Estimated  annual  ephemeral-stream  channel  recharge,  2001  and  2002,  Sierra  Vista  subwatershed,  Arizona 

[Group  numbers  correspond  with  groups  shown  on  figure  5.  hm3/km/2yr,  cubic  hectometers  per  kilometer  per  2-year  period;  hr/2yr,  hours  per  2-year  period; 
m,  meters;  m/hr,  meters  per  hour;  km,  kilometer;  hm3/2yr,  cubic  hectometers  per  2-year  period] 


Tributary  name 

Channel 

evaporation 

(hm3/km/2yr) 

Near¬ 

channel 

transpiration 

(hm3/km/2yr) 

Streamflow 
duration, 
relation  to 
distance  from 
mountain  front 
(hr/2yr) 

Wetted 

channel  width, 
relation  to 
distance  from 
mountain  front 
(m) 

Infiltration 

flux 

(m/hr) 

Total  channel 
length, 
overlying 
alluvial 
sediments 
(km) 

Ephemeral- 
stream 
channel 
recharge, 
2001  and  2002 
(hm3/2yr) 

Group  2 

Unnamed  Wash  1 

3.0- 10-4 

1.4-102 

35.5e‘°'07x 

0.15(0.2x+2.0) 

1.3 

15.4 

0.01 

Graveyard  Gulch 

do. 

do. 

do. 

do. 

do. 

14.4 

.01 

Coyote  Wash 

do. 

do. 

89  9e-o.o3x 

0.15(0.9x+1.9) 

do. 

18.5 

1.16 

Group  3 

Woodcutters  Wash 

do. 

do. 

do. 

do. 

do. 

18.5 

1.16 

Garden  Canyon  Wash 

do. 

do. 

35.5e‘°'07x 

0.15(0.2x+2.0) 

do. 

21.7 

.01 

Ramsey  Canyon  Wash 

do. 

do. 

do. 

do. 

do. 

17.4 

.01 

Carr  Canyon  Wash 

do. 

do. 

89.9e'0'03x 

0.15(0.9x+1.9) 

do. 

16.8 

.88 

Miller  Canyon  Wash 

do. 

do. 

35.5e-°-07x 

0.15(0.2x+2.0) 

do. 

16.2 

.01 

Hunter  Canyon  Wash 

do. 

do. 

do. 

do. 

do. 

15.5 

.01 

Stump  Canyon  Wash 

do. 

do. 

do. 

do. 

do. 

14.6 

.01 

Ash  Canyon  Wash 

do. 

do. 

do. 

do. 

do. 

12.9 

.01 

Brown-Bob 

Thompson  Wash 

do. 

do. 

do. 

do. 

do. 

9.5 

.01 

Group  4 

Lewis  Springs  Wash 

do. 

do. 

do. 

do. 

do. 

12.1 

.01 

Bakarich-McCool  Wash 

do. 

do. 

do. 

do. 

do. 

11.8 

.01 

Three  Canyons  Wash 

do. 

do. 

do. 

do. 

do. 

6.6 

<.01 

Palominas  Wash 

do. 

do. 

do. 

do. 

do. 

11.1 

.01 

Group  5 

Greenbush  Draw 

do. 

do. 

40.2e001x 

0.15(1. 8x-8. 2) 

.8 

19.5 

.9 

Group  6 

Banning  Creek 

do. 

do. 

2060.8ellx 

0.15(0.8x+0.9) 

.4 

8.6 

.15 

Unnamed  Wash  2 

do. 

do. 

do. 

do. 

do. 

10.6 

.15 

Stagg  Ranch  Wash 

do. 

do. 

do. 

do. 

do. 

10.1 

.15 

Little  Dry  Creek 

do. 

do. 

do. 

do. 

do. 

10.5 

.15 

Spring  Creek 

do. 

do. 

do. 

do. 

do. 

11.4 

.15 

Group  7 

Willow  Creek 

do. 

do. 

126.6e  007x 

0.15(1.6x-8.5) 

do. 

7.7 

.23 

Walnut  Gulch 

do. 

do. 

do. 

do. 

do. 

7.7 

.23 

Government  Draw 

do. 

do. 

do. 

do. 

do. 

7.7 

.23 

High  Knolls  Canyon  Wash 

do. 

do. 

do. 

do. 

do. 

7.7 

.23 

Total  ephemeral-stream  channel  recharge,  2001  and  2002 

5.89 

Average  annual  ephemeral- stream  channel  recharge 

2.95 
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Figure  29.  Approximate  distribution  of  annual  ephemeral-stream  channel  recharge,  2001-2002,  Sierra  Vista  subwatershed,  Arizona. 
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most  of  the  recharge,  occurred  during  the  summer;  however, 
stable-isotope  data  from  most  of  the  ephemeral-stream  channel 
boreholes  indicate  infiltration  from  past  summer  and  winter 
streamflow  events. 

Sensitivity  Analysis 

A  sensitivity  analysis  was  done  to  quantify  bounds  on 
the  calculated  ephemeral- stream  channel  recharge  rates  by 
independently  varying  the  recharge  calculation  parameters 
in  equation  1 1  within  reasonable  limits  and  recalculating 
recharge  along  the  length  of  each  channel  (table  13).  Channel 
evaporation  and  near-channel  transpiration  were  increased  and 


decreased  by  100  percent.  Streamflow  duration  was  increased 
and  decreased  by  50  percent.  The  wetted  channel  width  was 
increased  to  100  percent,  and  decreased  to  10  percent,  of 
the  total  channel  width.  Calculated  infiltration  fluxes  were 
increased  by  100  percent  and  decreased  by  50  percent. 

Recharge  rates  were  most  sensitive  to  the  wetted 
channel  width  (fig.  30  and  table  13) — average  annual 
recharge  for  the  subwatershed  increased  by  937  percent  with 
the  assumption  that  100  percent  of  the  channel  widths  are 
wetted  during  streamflow  events.  Although  there  is  a  high 
degree  of  uncertainty  associated  with  the  wetted  channel 
width  value,  the  assumption  that  the  total  channel  width  is 
wetted  during  streamflow  is  most  likely  not  valid.  Recharge 
was  less  sensitive  to  other  parameters  (fig.  30  and  table  13). 


Table  13.  Sensitivity  of  estimated  annual  ephemeral-stream  channel  recharge  to  calculation  parameters,  2001  and  2002, 
Sierra  Vista  subwatershed,  Arizona 

[Group  numbers  correspond  with  groups  shown  on  figure  5.  hm3/2yr,  cubic  hectometers  per  2-year  period] 


Tributary  name 

Resulting  ephemeral-stream  channel  recharge,  2001  and  2002 
(hm72yr) 

Channel 
evaporation 
(percent  change) 

Near-channel 
transpiration 
(percent  change) 

Streamflow 

duration 

(percent 

change) 

Wetted  channel  width 
(percent  of 
channel  width) 

Infiltration  flux 
(percent  change) 

+100 

-100 

+100 

-100 

+50 

-50 

100 

10 

+100 

-50 

Group  2 

Unnamed  Wash  1 

<0.01 

0.01 

<0.01 

0.23 

0.12 

<0.01 

1.34 

<0.01 

0.24 

<0.01 

Graveyard  Gulch 

<.01 

.01 

<.01 

.22 

.12 

<.01 

1.27 

<.01 

.23 

<.01 

Coyote  Wash 

.15 

1.16 

1.06 

1.26 

1.79 

.53 

8.30 

.74 

2.42 

.53 

Group  3 

Woodcutters  Wash 

1.15 

1.16 

1.06 

1.26 

1.79 

.53 

8.30 

.74 

2.42 

.53 

Garden  Canyon  Wash 

<.01 

.01 

<.01 

.31 

.15 

<.01 

1.79 

<.01 

.31 

<.01 

Ramsey  Canyon  Wash 

<.01 

.01 

<.01 

.26 

.14 

<.01 

1.52 

<.01 

.27 

<.01 

Carr  Canyon  Wash 

.87 

.88 

.80 

.95 

1.35 

.40 

6.29 

.56 

1.83 

.40 

Miller  Canyon  Wash 

<.01 

.01 

<.01 

.24 

.13 

<.01 

1.40 

<.01 

.25 

<.01 

Hunter  Canyon  Wash 

<.01 

.01 

<.01 

.23 

.12 

<.01 

1.34 

<.01 

.24 

<.01 

Stump  Canyon  Wash 

<.01 

.01 

<.01 

.22 

.12 

<.01 

1.27 

<.01 

.23 

<.01 

Ash  Canyon  Wash 

<.01 

.01 

<.01 

.20 

.11 

<.01 

1.14 

<.01 

.21 

<.01 

Brown-Bob  Thompson  Wash 

<.01 

.01 

<.01 

.15 

.08 

<.01 

.87 

<.01 

.16 

<.01 

Group  4 

Lewis  Springs  Wash 

<.01 

.01 

<.01 

.19 

.10 

<.01 

1.08 

<.01 

.19 

<.01 

Bakarich-McCool  Wash 

<.01 

.01 

<.01 

.19 

.10 

<.01 

1.08 

<.01 

.19 

<.01 

Three  Canyons  Wash 

<.01 

<.01 

<.01 

.10 

.05 

<.01 

<.59 

<.01 

.11 

<.01 

Palominas  Wash 

<.01 

.01 

<.01 

.17 

.09 

<.01 

1.01 

<.01 

.18 

<.01 

Group  5 

Greenbush  Draw 

.90 

.90 

.72 

1.11 

1.45 

.36 

7.15 

.54 

2.00 

.36 
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Table  13.  Sensitivity  of  estimated  annual  ephemeral-stream  channel  recharge  to  calculation  parameters,  2001  and  2002, 
Sierra  Vista  subwatershed,  Arizona — Continued 


Resulting  ephemeral-stream  channel  recharge,  2001  and  2002  (hm3/2yr) 

Channel 

Near-channel 

Streamflow 

duration 

(percent 

change) 

Wetted  channel  width 

Infiltration  flux 
(percent  change) 

evaporation 
(percent  change) 

transpiration 
(percent  change) 

(percent  of  channel 
width) 

Tributary  name 

+100 

-100 

+100 

-100 

+50 

-50 

100 

10 

+100 

-50 

Group  6 

Banning  Creek 

.15 

.15 

.10 

.21 

.25 

.05 

1.32 

.08 

.35 

.05 

Unnamed  Wash  2 

.15 

.15 

.10 

.21 

.25 

.05 

1.32 

.08 

.35 

.05 

Stagg  Ranch  Wash 

.15 

.15 

.10 

.21 

.25 

.05 

1.32 

.08 

.35 

.05 

Little  Dry  Creek 

.15 

.15 

.10 

.21 

.25 

.05 

1.32 

.08 

.35 

.05 

Spring  Creek 

.15 

.15 

.10 

.21 

.25 

.05 

1.32 

.08 

.35 

.05 

Group  7 

Willow  Creek 

.23 

.24 

.12 

.35 

.41 

.06 

2.21 

.12 

.58 

.06 

Walnut  Gulch 

.23 

.24 

.12 

.35 

.41 

.06 

2.21 

.12 

.58 

.06 

Government  Draw 

.23 

.24 

.12 

.35 

.41 

.06 

2.21 

.12 

.58 

.06 

High  Knolls  Canyon  Wash 

.23 

.24 

.12 

.35 

.41 

.06 

2.21 

.12 

.58 

.06 

Total  ephemeral- stream 

channel  recharge, 
2001-2002 

5.74 

5.93 

4.62 

9.74 

10.70 

2.31 

61.18 

3.46 

15.55 

2.31 

Average  annual  ephemeral- 
stream  channel  recharge 

2.87 

2.97 

2.31 

4.87 

5.35 

1.16 

30.59 

1.73 

7.78 

1.16 

Percent  difference  to  previ¬ 

ously  determined  average 
annual  ephemeral- stream 
channel  recharge 

-2.7 

.7 

-21.7 

65.1 

81.4 

-60.7 

937 

-41.4 

164 

-60.7 

RECHARGE,  IN  CUBIC  RECHARGE,  IN  CUBIC 

HECTOMETERS  PER  2  YEARS  HECTOMETERS  PER  2  YEARS 
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Figure  30.  Examples  of  sensitivity  analysis  results,  Greenbush  Draw,  Sierra  Vista  subwatershed,  Arizona,  2001-2002. 
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The  highest  and  lowest  annual  recharge  rates  determined 
through  the  sensitivity  analysis  (excluding  the  rate  determined 
for  100  percent  wetted  channel  width)  produce  bounds  of 
+7.78  hm3/yr  and  -1.16  hm3/yr  on  annual  ephemeral- stream 
channel  recharge  for  2001-2002. 


Summary  and  Conclusions 

As  part  of  a  larger  project  to  better  understand  the 
hydrologic  system  in  the  Sierra  Vista  subwatershed,  infiltration 
and  recharge  fluxes  through  the  sub  water  shed’s  ephemeral- 
stream  channels  and  basin  floor  were  investigated.  Data  used 
in  this  study  include  information  resulting  from  the  drilling 
and  coring  of  16  boreholes:  12  within  ephemeral- stream 
channels  and  4  on  the  basin  floor.  The  collected  data  includes 
physical,  hydraulic,  and  thermal  properties  of  sediment; 
chloride  concentrations  of  sediment;  pore-water  stable- 
isotope  values  and  tritium  activity;  neutron  soil-moisture  and 
EMI  logs;  continuous  sediment- temperature  profiles;  and 
continuous  water  levels. 

Infiltration  fluxes  through  the  basin  floor  were  0.2  to 
6  cm/yr.  At  the  basin-floor  sites  in  the  Walnut  Gulch 
drainage  basin,  long-term  infiltration  fluxes  were  very  low 
(less  than  1  cm/yr).  Chloride,  tritium,  and  stable-isotope  data 
indicate  long  periods  of  no  net  deep  downward  percolation 
flux  beneath  the  basin  floor.  The  sediment-chloride  profiles 
coupled  with  previous  studies  in  the  arid  Southwestern  United 
States  indicate  that  water  infiltrated  through  the  basin  floor 
does  not  move  downward  past  the  root  zone  for  long  periods 
of  time.  Infiltration  rates  could  also  be  negligible  in  other  areas 
of  the  basin  floor  that  have  soil  and  vegetation  similar  to  those 
at  WG4  and  WG5. 

Infiltration  fluxes  in  the  Carr  Canyon  drainage  basin  (2  to 
6  cm/yr)  were  higher  than  at  other  basin-floor  sites.  Chloride, 
tritium,  and  stable-isotope  data  indicate  active  infiltration  to 
8  m,  and  a  decrease  in  infiltration  below  8  m.  The  change  in 
the  infiltration  rate  below  8  m  is  controlled  by  an  increase 
in  the  silt  and  clay  content  of  the  sediment.  At  this  site,  the 
depth  that  infiltrating  water  moves  through  the  unsaturated 
zone  and  whether  the  infiltrating  water  will  ultimately  become 
recharge  are  strongly  dependant  on  the  silt  and  clay  content 
of  limiting  soil  horizons  in  the  unsaturated  zone.  Relatively 
high  infiltration  rates  could  occur  in  other  areas  of  the  basin- 
floor  that  have  soil  and  vegetation  similar  to  those  at  BF1, 
but  detailed  subsurface  soil  data  are  necessary  to  determine  if 
infiltrated  water  will  become  recharge. 

Infiltration  through  ephemeral-stream  channels  in  the 
Sierra  Vista  subwatershed  is  rapid,  occurs  to  depths  of  at  least 
20  m,  and  occurs  from  the  edge  of  the  mountain  front  to  as  far 
as  25.6  km  downstream.  During  streamflow  events,  infiltration 
fluxes  ranged  from  0.4  to  2.8  m/hr.  After  streamflow  ended, 
percolation  fluxes  ranged  from  0.1  to  2.8  m/hr.  In  areas  where 
low-permeability  layers  are  absent  in  the  shallow  subsurface, 
infiltrated  water  moves  past  depths  where  it  is  available  for 


evapotranspiration,  resulting  in  a  net  downward  flux.  This 
deep  downward  water  flux  is  assumed  to  ultimately  recharge 
the  regional  aquifer. 

Recharge  from  ephemeral- stream  channel  infiltration 
is  highly  dependant  on  the  thickness  and  hydraulic 
properties  of  the  stream  alluvium,  the  permeability  of  the 
sediments  underlying  the  stream  alluvium,  the  streambed 
geomorphology,  the  runoff  characteristics  of  the  contributing 
drainage  basin,  and  the  availability  of  streamflow.  Ephemeral- 
stream  channel  recharge  for  the  entire  Sierra  Vista 
subwatershed  was  estimated  by  upscaling  the  calculated 
infiltration  fluxes  to  the  area  covered  by  each  ephemeral- 
stream  channel  and  weighting  the  fluxes  by  streamflow 
duration,  evaporation,  and  transpiration.  The  resulting  average 
annual  recharge  during  2001-2002  is  2.95  hm3/yr.  Analysis 
indicates  that  the  calculated  recharge  rates  are  most  sensitive 
to  wetted  channel  width  and  are  less  sensitive  to  the  other 
parameters  of  evaporation,  transpiration,  streamflow  duration, 
and  infiltration  flux.  The  upper  and  lower  bounds  for  the 
annual  recharge  estimate  determined  from  the  sensitivity 
analysis  are  +7.78  hm3/yr  and  -1.16  hm3/yr,  respectively. 

The  calculated  annual  ephemeral- stream  channel 
recharge  rate  indicates  that  ephemeral- stream  channel 
infiltration  is  a  significant  contributor  to  total  recharge  in 
the  subwatershed.  When  streamflow  is  available,  recharge  is 
significant  along  the  entire  length  of  the  ephemeral- stream 
channels.  In  addition,  ephemeral- stream  channel  recharge 
occurs  during  both  summer  and  winter  streamflow  events. 
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